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SECTION 1 . 

INTRODUCTION AND SUMMARY 

Electric propulsion (EP) mission studies are usually conducted on 
the basis that thruster characteristics are given. The objective is to 
define the best mission using these characteristics. Over the past 
decade, the 30-cm electron bombardment thruster using mercury propellant 
has been the primary candidate for EP mission studies. The bulk of near- 

term mission possibilitip could be accomplished relatively well with the 
30-cm thruster. . 

Within the last few years, several new high energy mission possibil- 
ities have beau suggested. These include transportation of large apace 
systems from low Earth orbit (LEO) to geosynchronous Earth orbit (GEO), 

stationkeeping and attitude control of large apace systems, ’ 
planetary sample returns. For these missions, the 30-cm thruster 
operating at baseline conditions (i.e., 3,000-sec specific Impulse 
(Igp), 3 kW into the power processor) may not result in optimum mission 
performance. Considerations such as propellant type, number of modules, 
power level, and 1^^ suggest that a different thruster operating point 
would benefit several missions. 

This study of advanced electrostatic ion thrusters for space propul- 
sion was initiated to determine the suitability of the baseline 30-cm 
thru3ter for future missions and to identify other thruster concepts 
that would better satisfy mission requirements. In developing advanced 
thruster concepts, scaling and performance assumptions must Include 
technological realities. Thus, the general scope of the study was to 
review mission requiL-.-icnts , select thruster designs to meet these 
requirements, assess the associated thruster technology requirements, 
and recommend short- and long-term technology directions that would 
support future thruster needs. Preliminary design concepts for several 
advanced thrusters were developed to assess the potential practical 
difficulties of a new design. 
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A. OBJECTIVES 

The overall objective of this program was to develop recommendations 
for future ion thruster technology directions. Three specific objectives 
guided the study: (1) identification of the probable useful mission 
range of the existing 30-cm thruster, (2) investigation of methods and 
benefits of expanding the 30-cm thruster mission capabilities, and (3) 
investigation of methods and benefits of developing new thruster concepts 
(e.g., another size thruster) for those missions for which the 30-cm 
thruster is not well suited. 

To carry out the investigation implied by these guiding objectives, 
several detailed objectives were established. Including (1) selection 
and analysis of a representative set of potential ion propulsion missions, 
(2) assessment of mission performance as a function of thruster capabil- 
ity* (3) assessment of thruster technology boundaries, (4) identifica- 
tion of thruster or system parameters that scrongly influence mission 
design, and (5) development of a cost modeling technique to investigate 
the relationship between thruster parameters and overall mission cost. 
Relative to the last objective, only Earth orbit missions were considered, 
although, with minor modification, planetary missions could also be 
accommodated by the computer model. 

B. STUDY PLAN 

The study plan follows rather directly from the stated objectives. 
However, at the beginning of the study, the "working level" objectives 
had not yet been so clearly defined; they were, in fact, partially devel- 
oped during the study. The thruster technology assessment task was 
clear from the start. Mission set selection was reasonably clear but 
of rather broad scope. Since the analysis of the missions was deperdent 
on the mlsFlon set selected, the exact analytical approach was Initially 
difficult to define. Similar to the mission analysis task, the process 
of defining Important parameters was not straightforward because of the 
wide range of mission possibilities and mission objectives. The coat 
modeling task was developed during the study as the result of attempts 
to develop a generalized approach to Earth orbit missions. A logic 
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diagram of the various study tasks and their Interrelationships Is 
shown In Figure I**!; included In the diagram are the numbered tasks 
called for under the contract. Although the contract tasks ';ere not 
necessarily performed In the sequence originally planned, all contract 
objectives were satisfied. 

After a period that included data collection, selection of a pre- 
liminary mission set, analysis of several missions, and the start of the 
technology assessments, two significant study approaches became apparent. 
Tlie first, a generalized mission analysis approach* allowed a wide 
variety of available mission results to be displayed In a common format 
and provided a straightforward method for evaluating the sensitivity 
of "mission performance" to thruster characteristics. With the adoption 
of this technique, the mission set final selection process was greatly 
simplified since many missions could be handled easily and the problem 
of arbitrarily selecting or eliminating missions was essentially elim- 
inated. The output of the generalized mission analysis approach Is an 
estimate of the degree to which a given thruster technology could support 
the selected mission. 

The second approach that provided significant direction to the 
study evolved from an Earth orbit mission cost model technique. A com- 
puter program was developed to include cost models, mass models, thrust 
system performance models, and approximate relations for mission per- 
formance (l.e., "rocket equation"). For a given set of mission param- 
eters (e.g. , velocity requirement (&V), vehicle total mass, flight time), 
costs are computed as a function of specific impulse. This cost model 
approach provides a simple method fnr evaluating the relative sensitivity 
of total mission cost to a large number of thruster and system parameters 
(■*100). Most efforts during the last half of the study concentrated on 
the development and use of the cost model program. 

As Figure 1-1 Indicates, both the generalized analysis and the cost 
model approaches were supported by the technology assessment task. 
Performance scaling relations, physical scaling experience, and operat- 
ing limits (e.g., perveance, thermal) were Incorporated Into the 
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Figure 1-1. Study logic diagram 













analytical approaches and Into advanced thruster design concepts. 

Several design concepts for thrusters with larger beam areas were pre- 
pared and analysed in terms of technology advancement requirements. 

C. MISSION SET SELECTION AND ANALYSIS 

Initial efforts under the study were directed toward preparing a 
mission set to be used as the basis for the remainder of the work. The 
selection process, described in detail in Section 2. resulted in the 
following set of missions! 

• Earth orbit transportation 

• Earth orbit sta.t ionkeeping 

• Earth observatory 

• Outer planets (composite of seven missions) 

• Solar system escape 

• Mercury orblter 

• Close solar probe 

• Comet Encke rendezvous/ flyby 

• Comet Halley rendezvous/flyby 

• Asteroid rendezvous 

• Asteroid sample return 

• Mara sample return 

• Out-of-ecllptic. 

A generalized analytical approach was developed to assess the 
capabilities of several levels of thruster technology to perform these 
missions. The analytical approach utilizes mission trajectory calcula- 
tion results obtained from the literature and from work done under this 
contract. For each mission, characteristics such as AV, initial accel- 
eration (A^), initial mass, and payload mass were determined. These 
characteristics and an analytical formulation based on the rocket equation 
were then used to develop sets of curves in a AV vs A^ coordinate system 
using payload mass fraction (MP, payload mass divided by vehicle initial 
mass) as a parameter. Part of the analytical formulation required models 
of thrust system performance and mass. Thruster and thrust system 
technology enter the analysis through these models. An example of the 
AV vs A^ curve set format la shown in Figure 1-2. 

3y entering a given curve set (l.e., a set based on a given tech- 
nology) with mission data for AV and A^, the payload mass fraction required 
to perform a specific mission (MFR) can be compared with the mass fraction 
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PAYLOAD MASS FRACTION 



Figure 1“2 
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Le of AV vs A format and results, 
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capability of the assumed technology. The required mass fraction is 

Inpoeed on the curve map as indicated. Thus* if the AV va A coordinate 

o 

point occurs at a mass fraction higher than that required, the mission 
can be accomplished with the assumed technology. If the point falls at 
a lower mass fraction, the assumed technology level is not adequate. 
Since this approach Involves several approximations, the conclusions 
are not necessarily quantitatively accurate. VIhen the results show a 
definite capability margin or a definite inability to perform the 
mission, a relatively firm conclusion is possible. However, when the 
results are close, slight variations in assumptions could change the 
conclusion. Such precise Interpretation of the results is not warranted. 
Analysis details and results are presented in Section 3. 

D. EARTH ORBIT COST MODELING 

The generalized analysis approach provides a relatively simple 

method of gauging the ability of a given technology to perform a mission. 

However, several technologies potentially could satisfy the mission 
« 

technical requirements and remain within the boundaries of technology 
projections. The selection of one technology direction over another 
will certainly involve cost (technology development cost and recurring 
mission cost). 

Although this study does not attempt to predict technology develop- 
ment cost, significant effort was devoted to modeling the relationship 
of "transportation” cost to numerous technology-related thruster and 
system parameters. Since this study was primarily directed toward ion 
thruster technology, the scope of the cost modeling was limited to the 
propulsion function of a given mission. 

The cost model work was applied to several earth orbit missions. 
Earth orbit missions were selected because of the future potential of 
electric propulsion in the transportation and on-orbit support of large 
space systems. Selection of a propulsion system type (i.e., electric 
or chemical) will probably be strongly influenced by total system coat. 
Therefore, this cost model work was directed toward determining those 


7 


factors that elgnlftcancly affect ayste* coat aod toaaM detetnlBlBg the 
aeoaltlvlt, of ayatea coat to variations io throater, ayatem, and alaaloo 

parameters'. ' ' " 

The coat modeling analysis approach Is Illustrated qualitatively 
n Figure 1-3. The oMective of the calculations la to obtain total 
propulsion coat (transportation or on-orblt) as a function of various 
parameters. With total coat and varloua mas, hreakdowna. coat per unit 
of payload mesa (trenaportatlon) or coat per unit of net satellite mass 

otr vontrol) ere program out- 

puts. A conve..lent format for displaying the results la ahoUn In- 

Figure 1-4. Specific coat (?/k8 of payload) la plotted as a function of 
specific impulse ( 1 ^). „,th module power as a pncemeter. The nu*er „f 
modules needed to perform a given mlaalon varies with I and module 

power. Details of the analysis method and extensive results are pro- 
seated In Section 4. 


E. SUMMARY OF CffllCLUSIONS 

This study produced useful general methodologies for aasesalng both 
snetary and Earth orbit missions. For planetary missions, the assess- 
ment is in terms of payload performance as a function of propulsion 
system technology level. For Earth orbit missions, the asseasmsnt Is 

made on the basis of coat (coat sensitivity to propulsion system tech- 
nology level). 

The selection of the 30-c. thruster for near-term missions (for 

power levels below 100 hw) was reinforced by this work. Except for hlgh- 

power missions or those requiring a propellant other than mercury, the 

30-cm engineering-model thruster :EMr) is a good choice. However, for 

systems larger than about 100 kW (Earth orbit or planetary) . a larger 

t roster sire would be advantageous In terms of cost and payload 
performance. ^ 

Based on the ■Avalon studies, cost modeling, and technology 
ments. a 50-cm-dIameter thruster Is Suggested as the next step In 
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Figure 1-3. Earth orbit mission cost modeling, simplified 
block diagram. 
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transportation costs per kG OF PAYLOAD, $/kG 
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Figure 1-4. Example of cost modeling results. 
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thruster technology development. By proper selection of propellants 
and beam currents, a 50'-cm thruster designed to operate with a beam 
voltage of about 2A00 V would satisfy moat of the requirements of 
future missions. Thruster concepts. Including an oval cross section, 
were developed to illustrate design techniques for thrusters with larger 
beam areas. 
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SECTION 2 

MISSION SET SELECTION 

The goal of the mlsalon set selection process was to establish a 
representative set of 10 missions that would fairly reflect the range of 
future electrostatic Ion thruster requirements. These requirements pro- 
vided direction and motivation for future thruster technology develop- 
ment. The general selection process followed Is Illustrated in 
Figure 2-1. 

A. CONTRACT GUIDELINES 

Guidelines for selecting 6 of the 10 missions were provided In the 
contract. Specifically, the mission set was to Include: 

• Transportation of large space systems from LEO to GEO 

• Stationkeeping and attitude control of large space systems 
in GEO 

• A mission to less than 0.7 AU 

• A solar system escape mission 

• At least two missions between 0.7 and 5.0 AU, 

These mission categories Indicated the desired scope of the study without 
being overly restrictive. 

To aid in the mission selection process, NASA provided several 
reports on previous electric propulsion studies and a list of documents 
to be considered in establishing the mission set. The document list 
used is presented in Appendix A. This background information provided 
results of trajectory calculations for various missions, system concepts, 
and designs and the general sensitivities between thruster requirements 
and mission objectives. Although a certain amount of mission/trajectory 
analysis was anticipated, it was assumed that the bulk of the -nlBslon 
data needed in the study would be obtained from existing literature. 



HRL 

SCREENING 


NASA 

INPUTS 


PROCESS 


Figure 2-1. Mission set selectlori process flow diagram. 
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B. MISSION CANDIDATES 


The mlaalons initially cocaldered for selection are shown in 
Table 2-1 and are characterized In terms of payload mass range and AV 
range. The relatively wide spread In the parameters of Table 2-1 
reflects the flexibility of mission requirements, mission design, and 
utilization of electric propulsion. Clearly, a specific mission design 
would require defining many other factors. However, AV and pay loqd mass 
tend to allow a first-order examination of propulsion system requirements. 

Candidate missions are simmarized briefly here to provide a basis for • 
evaluation. 

Transportation of Large Masses 

The types of missions that could be considered In this category 
include transportation of full Shuttle-size payloads, portions of large 
sr«ce power generation satellites from LEO to GEO, transportation of 
space stations, or recovery of large satellites for repair or 
disposal. Propulsion requirements might include single one-way- 

trip-to-orblt, single round trips with refueling, or multiple round trips 
with and without refueling. In some cases, power might be provided by the 
payload, while in others (e.g., a return to LEO) the EP transportation 
vehicle would require its own power source. EP Is a prime candidate for 

this category of missions and will probably benefit from advanced 
technology. 

^.atlonkeeplng and Attitude Control of Large Systems 

Massive and/or large-surface-area systems in Earth orbit will 
experience large torques and disturbance forces. Correcting these 
attitude and orbit perturbations could require a significant propulsion 
capability on-orbit. The specific mission of the satellite will dictate 
the detailed requirements in terms of allowable eccentricity change or 
pointing variations. However, propulsion system hardware and propellant 
for on-orblt operations represent a penalty to the satellite since this 
material must be transported to and stored on orbit. In addition. 


Table 2-1. Candidate Missions and Characteristics 


Mission 


Earth orbit 

Transportation of 
large masses 


Stationkeeping and 
attitude control of 
large systems 

Multimission modular 
spacecraft missions 

Earth observatory 

Planetary 

Out-of-ecliptic 

Close solar probe 


Solar system escape 
Mercury orbiter 


Mercury lander 


Mercury sample return 


Venus orbiter 
Veniis lander 
Venus sample return 
Mars sample return 
Jupiter orbiter 
Jupiter satellite orbiter 
Saturn flyby 


Approximate Payload 
AV,® km/ sec Range, kg 


10 ^- 10 ^ 


10^-10^ 


Comments 




5x10^-10^ 


10^-2x10^ 


3-5x10' 

10^-2x10^ 

10^-2x10^ 


Thermal 

Constraints 

Jupiter swingby 

Thermal 

constraints 

Thermal 

constraints 

Thermal 

constraints 


EP has little 

advantage 

■over 

ballistic 

mission 





tabu 2-1. candidate Cbar.curUtlc. (ConUnued) 
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Mission Approxl 

AV,® km 

Planetary (Continued) 

Saturn orbiter 

Titan orbiter - * 

Uranus flyby ^^3 

Neptune flyby ^_3 

Asteroid rendezvous 6 - 1 ' 

Asteroid sample return 10-15 


Approximate Payload 
AV, km/sec Rangei kg I 


Comments 


5x10-10^ 

10^-3x10^ 

SxlO^-10^ 



5x10^-10^ 


Flyby velocity 
dependent 

Nuclear, using 
payload as 
power source 




pHuosop.,. J^LTLlmT 

i«l« th.t should mlnisass the ou-orh,! 

vehicles «y ^ tech°T°^°” ‘^p*' 

operation. technology for coat-effective 

Modular Cp . — ^rift 

The multtoisslon nodular spacecraft fMMS^ 

satellite bus that, with certain mla 9 < ^ ^ Provide a basic 

Potentlsl tor algnificeutlv reducl "»0"lea. hes the 

Althoua, EP could be hene£lclelly appUeT"'* ppp‘“-*‘ 

syateu re,ulremsnts can ha adecuat.l “ ‘ *'’® P">P“l»lon 

8- thrustets. Bacaua. o/th T' ’“-P- 

Phv HHE is unllhaly t. ta,ulte 

^-arth O bservaroT-y 
8aalrea.21 

=vyatal dynanlcs. ocean dynaMca and t a 

this type night also Include •’spleCl A satellite 

»°VW Ptobebly ba synchronous for atotlor' ‘T^ *’“'tt™«>ta. The orbit 
satellite capability for orbit measurements; a 

-..■a.. 

bla of operating for nany years with “’’“"story satellite, capa- 

ba substantially larger than ordinary slnT ‘““""“ts. would 

should be well suited for orbit ral.J “stellites. BP 

““““"■'“P‘»8. .SB positioning. 

Transfers out of the ecliotlr 

spsce scientists because of the n to 

such an orbit would provide **W-28 “"PSttlve on the solar system chat 

bscause of the high by. the'unconstralned^lTn^" 

IncUnatlon change (as opposed to a , 8radual 

Ihe power level and satellite size d transfer), 

complement. ““”8 largely on the Instrument 
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6. nose Solar Probe 


The close solar probe is a good EP mission candidate because of the 
relatively high AV and th. effective reduction in propulsion system 
specific mass resulting from increased solar intensity. Such a 

satellite is expected to be relatively small, requiring a propulsion 
power level in the range of 15 to 25 kW (1 AU). However, "advanced- 
thruster requirements might result from thermal constraints. 

7, Solar System Escape 

The degree to which EP is needed for this mission depends on pay- 
load mass. Using a Jupiter swing-by trajectory and relatively long 
flight times, a reasonable payload could be transported ballistlcally. 
Electric systems in the 15- to 25-kW range improve payload capability 
and/or reduce flight time. However, since solar power is drastically 
reduced beyond 5 AU, larger EP systems probably would not be required 
unless a nuclear power source was used. 

8. Mercury Orbiter 

This mission Is similar to the close solar probe mission, with 
somewhat relaxed thermil constraints. * Typical payloads of a tew 
hundred kilograms will require propulsion power levels of 20 to 25 kW. 
Thus, except for possible thruster thermal limitations, existing tech- 
nology is probably adequate. 

9, Mercury Lander 

The payload (transported by an EP system) for a Mercury lander 
mission probably would be about an order of magnitude greater than that 
for an orbiter mission. An electric propulsion system power level of 
100 to 150 kW would be needed to transport such a vehicle (initial mass 
'X. 4 X 10^ kg) . Existing 30-cm thruster technology could be applied to 
such a large system, but it would probably be less effective than a 
system using a higher power thruster. 


10. Mt> r curv Sample Return 

This mission voul'i need about two orders oT magnitude more pay loud 
than an orbiter mission. Thus, the KP system would be propelling a 
vehicle having an initial mass in the range of i to 4 x 10^ kg and would 
requite a power level on the order of 10^ kW, Such a system, which Is 
in the range of "large space systems'' being considered under Mission I, 
would probably benefit significantly from the use of an advanced 
thruster. 

1 1 . Venus Orbiter 

The relatively short transfer time to Venus reduces the usual poteti- 
tlal advantages of EP, If the mission were performed with EP, the exist- 
ing 30-cm thruster would probably be sufficient. 

12. Venus Lange r 

The Increased payload required for a lander would make EP more 
advantageous, but probably not mandatory. 


13. Venus Sample Return 


With the Increased flight time demanded by a return flight and 
possible maneuvers in orbit at Venus, an EP system would probably be 
beneficial. The large mass and '.ii.gh power level might require an 
advanced thruster. 


• M ars Sample Return 

Compared with a ballistic approach to this mission, EP can increase 

the payload returned to Earth from a few hundred grams to several kllo- 
22 

grams. Typical missions would use a shuttle launch, EP for Earth/Mars 
transfeis, and chemical systems for Mars descent and ascent maneuvers. 
Mission times of three to four years and power levels of 50 to 100 kW 
would be needed. 


15. Jupiter Orbiter 

If performed using EP, this mission would require power levels of 
15 to 25 kW for typical payloads. Larger payloads (for instance. 
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if surface probes were included) might require higher power levels. 

Thus, conventional mission approaches of the type considered in many EP 
studies would require only existing technology. More aggressive missions 
might benefit from advanced technology. 

16. Jupiter Satellite Orblter 

Limited reconnaissance of Jovian satellites will probably be per- 
formed as part of Jupiter orblter missions. However, more extensive 
satellite investigations could require larger payloads. For such vehi- 
cles, advanced EP systems would be beneficial. 

17. Saturn Fly-By 

EP improves the payload capability or reduces flight time 
missions of this type with only modest power levels (15 to 25 kW). 
However, with solar-powered systems, propulsion is limited to distances 
of 4 to 5. AU. Significant payload improvements will occur when space 
nuclear power sources become available. Such systems in the multi- 
100-kW range would benefit from advanced technology thrusters with long 

life at high power levels. 

18. Saturn Orbiter 

See comments on Mission 17 . 

19. Titan Orblter 

See comments on Mission 17 . 

20. Uranus Fly-By 

See comments on Mission 17. 

21. Neptune Fly-By 

See comments on Mission 17. 

22. Asteriod Rendezvous 

Several asteriod missions have been studied 
be well within existing EP technology capabilities. » * ’ Because 

of the high AVs associated with asteriod missions and the reasonable AU 
distances, EP is ideally suited to such missions. 
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23. Asterlod Sample Rettirn 


The higher &V associated with such a mission improves the advantage 
of an EP system. The larger masses and higher power levels 
required, compared with a rendezvous, would increase the size of the 
propulsion system, but would probably not require advanced thruster 
technology. 


24. Comet Encke Rendezvous/Fly^By 

This mission is typical of several cometary missions, but is of 
particular interest because of its three year orbital period and near 
earth passages. Encke is a good EP mission 

because of the high AV. However, various navigation errors might make 
a true rendezvous rather difficult. With typical science payloads, pro- 
pulsion power levels of 20 to 25 kW would be adequate for this mission. 
Advanced thruster technology should not be needed. 


25. Comet Halley Rendezvous/Fly-By 

13-16 

Comet Halley is of both scientific and popular interest. 

Because the comet's orbit is retrograde, the mission is essentially 
Impossible chemically. It would, however, be quite possible with EP. 
Significant factors in evaluating EP system requirements are flight time 
and fly-by velocity, -iven some flexibility in these parameters, exist- 
ing technology would be adequate for the comet Halley mission. Unfor- 
tunately, the time available to develop a spacecraft for the 1986 passage 


is very short. 

26. Nuclear Waste Disposal 


Strictly from an EP point of view, the nuclear waste mission is 
technologically feasible. Since power would be derived from the payload 
and fli^t time would be ijnconstralned, the existing 30-cm thruster 
would probably be adequate. 


C. MISSION SELECTION 

Based on considerations of the type presented in Section 2.B, dis 
cusslons with the NASA LeRC study manager, and contract guidelines, a 


Table 2 2. generally eatisfled all re,„lre«„ta. Hewever. Inprevenenra 

were recommended as the study proaressed ai#.k u u ^ ” 

set ahn™ 4 n. U, . Progressed. Although the final mission 

. n able 2-3, includes more than the 10 missions required it 

la eere representative el tbe study goals than „as the Initial set 

In the final nlsslon set, the various outer planet missions uere 
aken as a group because their requlrefents are slrtlar. The uultl- 

oTthTpTlT 1-a bulk 

payloads represented by this category require low dVa. have low 

mass, and would most likely use exist^no t-h 

sample return mission and the comet Halley MllIL'wIrHZd b'*"" 
bath were of current Interest and both wight require T alee T' 

thruster technology. The close solar probe was added because It was . 
extension of the mercury orbiter mission. " 

Each mission selected was analyzed from the standpoint of its 
man s on thruster and propulsion system technology. These demands 

Id requirements, indicate the de;ree 

c on o technological change needed. At various points in the 

extrapolations were not unreasonable. 
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Table 2-2. Initial Mission Set 


Mission 

Comment 

Earth orbit 


EO transportation 

Contract requirement 

EO stationkeeping 

Multlmlsslon modular spacecraft missions 
Earth observatory 

Contract requirement 

Planetary • 


Out-of-ecllptlc 

Contract requirement 

Solar system escape 

Contract requirement 

Mercury orblter 

Contract requirement 

Neptune flyby or orblter 


Asteroid rendezvous 

Contract requirement 

Comet Encke rendezvous 




Table 2-3. Final Mission Set Uaed in Study 


Mission 

Comment 

Earth orbit 


transportation 

Contract requirement® 

EO stationkeeping 

Earth observatory 

Planetary 

Contract requirement 

* 

Outer planets 

Jupiter orblter 

Jupiter satellite orblter 


Saturn flyby 

Saturn orblter 

Titan orblter 

Uranus orbiter 

Neptune flyby 

1 

^Considered as a single type of mission 

Solar system escape 

Contract requirement 

Mercury orbiter 

Close solar probe 

Comet Encke rendezvous 

Comet Halley rendezvous 

Cont ract requirement 

Asteroid rendezvous 

Asteroid sample return 

Mars sample return 

Contract requirement 

Out-of-ecliptic 

Contract requirement 

a 

See Section 2. A for discussion of requirements. 
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SECTION 3 

MISSION SET ANALYSIS 

To use the information gathered from the various reference sources 
and the trajectory data generated by this study, a generalized analysis 
technique was developed. All the pertinent data was converted to a com- 
mon format to provide a uniform basis for comparison. The technique is 
discussed in Section 3, A. Section 3.B presents propulsion system perfor- 
mance and mass models used in the analysis, including estimates for 
several technology levels. Section 3.C summarizes the mission character- 
istics required for the generalized approach, including AV, Initial 
acceleration, payload, and initial mass. Section 3.D discusses the 
generalized analysis. In addition, results for specific missions ana- 
lyzed in detail (e.g. , trajectovy calculations) are summarized. Detailed 
information that does not fall in the "generalized analysis" category is 
Included in the Appendix B. Conclusions from this portion of the study 
are presented in Section 3.E, 

The generalized mission analysis technique was applied to the total 
mission set, including the Earth orbit missions. However, the transpor- 
tation and stationkeeping of large systems. Missions I and 2, were not 
clearly defined. That is, since wide ranges of mass, size, and power 
level were to be considered, with only broad mission objectives, the 
generalized analysis could not provide all the desired results. In 
particular, large-scale operation in Earth orbit will require cost- 
effective propulsion systems. The economic tradeoffs related to 
Missions 1 and 2 were developed through cost modeling techniques 
(discussed in Section 4). 

A. GENERALIZED ANALYSIS APPROACH 

1. Results Format Basis 

Mission analysis results obtained from various studies are generally 
sojaewhat difficult to compare because of the multitude of assumptions 
required in the basic calculations. However, by distilling the pertinent 
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information and displaying it in a common format, as illustrated in 
Figuto 3-1, approximate comparisons can be made. The basis of this for- 
mat is derived from the rocket equation with an appropriate use of EP 
system definitions. 

The rocket equation has the form 

• ( 3 -« 

where 

AV - velocity change characterizing the mission 
“ gravitational constant 

Igp true specific Impulse (corrected for multiply charged ions, 
beam divergence, and propellant utilization efficiency) 

M • Initial mass of the vehicle 

O' 

“final mass of the vehicle. 

The initial mass can be written as 


- M- + M 
o f p 


(3-2) 


or 


where 


M 


M 


pi 


M + 
ps . 


M 


( 3 - 3 ) 


M 


M 


pi 


propellant mass 
payload mass 


**ps “ P^oP“i8i-t>n system mass. 


A simple reartangenent of Eqs. 3-2 and 3-3 yields 


M 


H 


H 


M 




( 3 - 4 ) 


H. 

^ = Mp + -JES. 
o o 


(3-5) 


or 


AV, km/aac 


60 


72S4-22iKa 



Figure 3 - 1 . Generalized mission analysis results format. 
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where MF payload mass fraction. 

EF system parameters are related to propulsion system mass as 
follows: 


ps 


M 

_££ 

P • 
ps 


(3-6) 


where P__ “ propulsion system Input power, or 
ps 


M 

o 


P a 

M 

o 


(3-7) 


Power Is related to other parameters by means of the energy equation: 


ps 


TI g 
sp *^o 

2 n * 
ps 


(3-8) 


where 


T 


True thrust (corrected for multiply charged Ions and 
beam divergence) 


n “ True propulsion system total efficiency (corrected for 
^ multiply charged Ions and beam divergence) . 


This allows us to write 


M TI g 

-BS. „ 3P o PS 

o o ps 


(3-9) 


In terms of vehicle Initial acceleration, 

A T 

A “ 77 — 

° Vo 


(3-M) 


For convenience In using mission analysis literature, has the units 
of "g's." Eqs. 3-8 and 3-9 become 


A I g 

p „ „ I P . PB ^ 

ps o\ 2 npg 


(3-11) 


M A I ot g 

ps o SP PS ^o 

2n 

o 'ps 


(3-12) 


If is expressed in conventional units of kg/kW, then 


M 

-28 - 
M 
o 


Ala 
O sp PS 


20.79 n 


(3-13) 


pa 


Substituting Eq. 3-5 into Eq. 3-1 and using Eq. 3-*13 yields 

" - So V^“ (” + 

As illustrated in Figure 3-1 for a given set of propulsion system 

parameters that are technology- level related, lines of constant pfyload 

mass fraction can be defined in a AVvsA coordinate system. To use this 

o 

format, propulsion system models were developed for the desired technol- 
ogy level and were used iii displaying Eq. 3-lA graphically. The result 
is a graph having lines of constant mass fraction that represent a given 
propulsion technology level. Mission data is then used to obtain 

typical values of AV, A , and MFR. 

c 

The mission study value of MFR is then imposed as a line on the 
appropriate graph. This line represents the payload mass fraction typi- 
cally required to satisfactorily perform the mission. Then, the study 

values of AV and A are located on the graph. If the AVvsA point is at 
o o 

higher mass fraction th 2 in that required, the mission probably can be 
performed with the assumed technology. If the AVvsA^ point is at a 
lower mass fraction than required, the mission probably cannot be accom- 
plished with the assumed technology. 


2. Number of Thrusters 

An estimate of the number of thrusters required for a given mission 
and an assumed technology level can be obtained as follows. The maximum 
number of operating modules (one thruster and power processor combination) 
is approximately the ratio of propulsion system input power to the total 
power per module (i.e., power into a power processing unit (PPU)): 


N 


mod 



(3-15) 
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where power preceseor end 


■pp 


Vc ''eVc 


where 


" thruster input power 
Ij, - beam current per thruster 
\ " beam voltage 

^pp " power processor efficiency 
Hg - thruster electrical ef^ciency 
\ - cabling efficiency. 


The been perenerers ere deterwined by the technology choice 
I„ depende on thruster design, end depends on I fron 


V, « 

D 


(») 


sp 


where 


*■ ion mass 

e = electronic charge 


Oy = measured propellant utilization efficiency 
Y " thrust loss parameter. 

Combining Eqs. 3-11, 3-15, end 3-16 yields 


N 


mod 


( 4 ) (^0 


A I g 
o sp ®o / 


because 


"ps " VuHppV . 


it follows that 


N 


mod 


/M A \ /n \ 
1 - 0 o t /e\ iW 

\ K } vn»/ ll ) 

' b / ' \ sp/ 


32 


(3-16) 


since 


(3-17) 


(3-18) 


f 


This Is the approximate number of operating modules needed to utilize 
the propulsion system power. 

Two additional conditions must be satisfied: (1) there must be 
enough modules provided to account for wearout, and (2) standby redun- 
dancy must be provided independently of wearout. The total number of 
inodules required for Hearout (now) oen be eetlmated from total propellant 
requlrenenta. being Eqe. 3-2. 3-5.and 3-13, the propellant mane la 

Mp - M„ Jl - (mf + 

and 

Mp -y^Mdt N__^(t)dt . (3-20) 

where M =* total mass flowrate. 


The value of the Integral In Eq. 3-20 la almply the total nunber of 
anpere-houra of operation required to uae the propellant. If a thruater 
is characterlted by a wearout life In ampere-hours . Independent of 
throttling, the total propellant mass can be simply divided by the mass 
equivalent of the ampere-hour life. Thus, 

”00 ’ (m) ’ (3-21) 

where AH = thruster life in ampere-hours. 


To establish the total number of thrusters required, N and N 
must he scaled up to account for redundancy/rellablllty. To'first ™der. 
a simple redundancy factor. E^, can be used. Ihen, the approximate 
number required is the larger of 


N 


ot 


No. of thrusters 
required 



j(l + R^) 
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(3-22) 


B. PROPULSION SYSTEM MODELS 


The generalized analysis technique discussed above reqOitres models 
for propulsion system specific mass and totaJ efficiency as 

Indicated by Eq. 3-14. Models for these parameters are presented below. 

1. Performance Models 

Although the desired final result of this section Is a model for 
Hpg, several other relationships must also be defined. Total efficiency, 
as Illustrated in Eq. 3-8, originates from the energy equation and Is a 
product of ail thruster, power processor, and cabling efficiencies: 


ps 


•^t ’^pp ^c 


where 


(3-23) 


Hj. = thruster total efficiency 
Hpp = power processor efficiency 
rip = cabll’ag t’fflclency. 

Thruster total efficiency includes electrical efficiency, propellant 
utilization efficiency, and factors to account for ion beam divergence 
and ir,u-.{;iply charged ions; 

2 


= measured propellant utilization efficiency 

= n, + 2n,, 

rij mass .fraction flowrate of singly charged ions 
’i2 = mass fractional flowrate of doubly charged ions 

Up = thruster electrical efficiency 



Y = experimentally determined correction factor for beam 
divergence and multiply charged ions 
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“ beam voltage 

Ijj » beam current per thruster 
" thruster Input power. 

Beam voltage is related to other paraii 


energy: 


parameters through the definition of ion 


since 


2 "'i \ 


' ■■ (*y 


( 3 - 24 ) 


( 3 - 25 ) 


where 


I - -!L 
SP 




“ ion mass 

e = electronic charge 


= ion velocity due to acceleration through V 
— b 

V - average exhaust velocity 

■ W- 

Beam current cen be related to flowrate through 


where 


'b ■ "1 (^) = K) . 


( 3 - 26 ) 


( 3 - 27 ) 


( 3 - 28 ) 


ion mass flowrate 
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Calculations made for the generalized analysis utilized the following 

empirical relationships associated with the 30-cm-diameter Hg ion thruster 

\ 43,59,60 
(in mks units) : ’ 



I = neutral atom flowrate (equivalent amperes) 
n 


6 


Y 


0.24 + 0.032 


1 - 0.08 




+ bx 


2.2 + I 


bx 


bx 


bx 


0.942 - 0.005 Ii, 

b 


‘ 0 for I. i 1 
b 


(3-29) 

(3-30) 

(3-31) 

(3-32) 

(3-33) 





1 - 


0.13 

,0.5 


(3-34) 


Substituting these terms into Eq. 3-23 is ouitted since the equation 

would be somewhat awltward. But, to summarize, n can now be evaluated 

P“ 

for a given as follows : 

• from Eqs. 3-30, 3-31, and 3-32 

• n from the definition and Eq. 3-29 

• Y from Eq. 3-33 

• n n from Eq. 3-34. 

PP c 

Propulsion system total efficiency is plotted in Figures 3-2 and 3-3. 

The apparent slight decrease in efficiency that occurs at high beam 
current is a consequence of the Increased influence of multiply charged 
ions as modelled in equation 3-33. It should be noted that equation 
3-33 has been formulated to fit data in the 1- to 4-A beam current range 
and the accuracy of the expression for lOA has not been verified. 




V 


I 

I 


I 
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BEAM CURRENT. A 


Figure 3-2. Propuisifi.i system total efficiency for 30-cjn- 
dlameter mercury thruster as a function of , 
beam current. 
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EFFICIENCY 
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Figure 3-3, 


Propulsion system total efficiency as a fui ction of 

I . 

sp 
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Although th. omclency ooloulotton. «. basod oo data obtaload 
fto. oparatlon ot lO-om thruatara on Hg.*’*“-“, ualng thaaa oquationa 
for doactlhtog thruator parfomanca ualng othar propaUanta ot for othar 
aiaa thtuatata can ha axpactad to produca ralatlvaly accurat, result,. 
More apaclflcally, thtuatar electrical efficiency »1U acala exactly; 
utUlaatlon efficiency may be aonaubat In error but not oorlously; 

Y should acala roughly aa defined by gq. y-M for other propellanta. but 
for other alaea an equivalent 30-cm bean current should be uaed; and 

pouer ptoceaaor efflcl«,cy (Eq. 3-M) la probably reasonable lor beam 
voltages greater than 1 kV. 

The general dependence on propellMit type la through atomic maaa. 

y notmaUtlng the atomic mass to mercury, the V. and L relatlonahlpa 
can be written In the form ® 



(3-35) 

(3-36) 

(3-37) 


y ■ atomic mass unit • 1.66056SS x 

W - atomic mass ' 
yW » 

X » Atomic m ass of propellant 
atomic mass of siercury 

■ I (mercury) 

■ 0.19915 (argon) 

• 0.41777 (krypton) 

*“ 0.65457 (xenon). 


X - 39.948 

A 200,59 

X . -83.8 

Kr 200.59 

X - ^31.3 

Xe 200.59 
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Additional relationships that are of use In various calculations 
Include thrust T and vehicle acceleration A; 



(3-38) 


■ mV ' - 

o o \ o sp / 

where M^ le Che vehicle Initial mans and all values are In mka unite 
except A^, which is expressed in "g's#" 


2. Mass Models 


Propulsion system specific mass a is assumed to include all the 
elements required to produce thrust, including a power source (e.g., 
solar array). Historically, o's of 30 kg/kW were used in mission studies 
based on optimistic technology goals. Generally, these a's included only 
the element masses and not system considerations such as thermal control 
and vehicle structure required for support. Depending on the particular 

study, such system factors may or may not have been Included in the 
overall mass accounting. 

The mass models discussed in this section include allowances for 

the major system factors assuming certain design approaches. Models for 

present technology are based on hardware currently under development^^' 

(specifically, the 30-cm EMT and its corresponding power processor). For 

the generalized analysis, advanced technology was modeled by projecting 

the capability of the 30-cm mercury thruster to higher beam current and 
voltage. 


mass models used to define present and advanced thruster technol- 
ogy are shown in Table 3-1. Although many system designs exist, these 

Tl ““e based on typical values used In the various referenee docu- 
ments and should be reasonably representative. 
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Table 3-1. propulsion System Hass Model Used In Generalized 
Mission Analysis 


Sub-Assembly 

Present 
Technology 
Mass, kg 


Solar array* 

18 Pq 


Power conditioning panel*^ 

(Including structure and thermal control) 

+ 3 Ppp 


Thruster 

(Ineludes gimbals and structure) 

17 


Tankage, structure, and miscellaneous^ 

5^0 


JKBQuXlOBIICy XttCifOlCr . ^ 

Power conditioning panels 

1.2 


Thrusters, 

1.4 


Lifetime 

c 

a. 

Pq Is solar panel power at 1 AU after degradation. 

b. 

Pppls Input power to power conditioning panel (maximum). 

c. 

‘Constant A-hr wearout* 
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The. sour arroy apeolfio mass la baaed on the aaaomptlon that all 
depredation occurs earl. In the nlaalon. The IS hg/hW aaaunea approa - 
nately a 15 to 20T degradation of the original Inatalled power (15 kg/ )• 
Since aolar array atudlea Indicate that apeclflc^muea on the order o 
5 kg/kw aay he feasible with present technology, the res a o 

work discussed here arc probably conaeryatlve. . 

Power processor mass was modeled using historical data for PPlIs alone. 

and current design data for the series resonant Inverter (Ski) 
developed by NASA leKC, which Includes thermal control and structure. 

The basis for the equaVlon used In Table 3-1 la Udlcated U 
The “hlstorlcalP curve was obtataed as shown In Figure 3-5 using data for 
various PPNs without structure and thermal control. Since the mass of 
the packaged FFU Is dependent on many vehicle design factors, the straight 
line shown In Figure 3-4 was assumed as a reasonable estimate. 

Thruster mass, as used In the model. Includes an allowance 
gimbals and miscellaneous structure associated with the thruster. Typi- 
cally, the total "thruster array” Is about twice the weight of the 
thrusters.^**’ All remaining propulsion system elements (such as tanks, 
valves, structure. miscellaneous hardware) are assumed to be proportions 

to power at the rate of 5 kg/kW. 

As Table 3-1 indicates, thruster mass is constant in this formula- 
non. However, since It Is assumed that higher power can be obtained, 
thruster lifetime and the required number of modules plays a strong ro e 
in modeling propulsion system specific mass. Figure 3-6 Illustrates the 
influence of lifetime on mission performance. In the first example, the 
on operates at 2 A. The total mass of two thrusters and their 

power conditioners Is about 136 kg. Modifying the “ “ “ 

a 4-A beam current, but keeping specific Impulse at 3000 sec (modified EMI), 
yields a total mass of 77 kg without considering thruster lifetime (l.e. . 

If no wearout of the thruster occurs). This Urge decres=» In mass la 
«lnl, due to the lower specific mass of the FFO at the higher power 
levels. The trajectory dynamic performance U unchanged for this second 
example. sUce running fewer thrusters at higher current can produce the 
same acceleration time profile with no Increase In system propellan o 
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g POWER PROCESSOR INPUT POWER, kW 

Figure 3-4. Power processor mass model. 
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Figure 3-5. Power processor 


historical mass data. 
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EX. 

EMT 



POWER 



PROCESSOR 

THRUSTER 


LjlJ — 

-*|aj 


Mpp '»1O2k0 


m — 


- 34 kg 


EX. I40WEAR0UT 

MODIFIED EMT (1^ « 3000 SEC) 
POWER 

PROCESSOR THRUSTER 
4« I - . » f 4i 


Mpp *60kg 
Mj -17 kg 


EX. CONSTANT A-HR WEAROUT 

MODIFIED EMT (l,p - 3000 SEC) 
POWER 

PROCESSOR THRUSTER 



Mpp -60 kg 
M| - 34 kg 


NOTE: UTILIZATION EFFICIENCY AFFECTS WEAROUT 


Figure 3-6. Exan^le of mass modeling philosophy applied 
In generalized analysis. 
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power requlremente. The only chenge le a lighter propulelon ayseem with 
potentially Improved payload. 

The mass modeling of thruster wearout lifetime is assumed to occur at 
constant ampere-hours. That is. if the thruster lasts for a given duration 
it 2 A, it will wear out in half that time when operated at 4 A. As in 
♦•.le third example of Figure 3-6, a spare thruster must then be carried. 
Switching occurs for wearout in a manner analogous to switching for 
thruster failures. The total number of thrusters is the same as in the 
first example, but a significant mass savings accrues from raising 
module power. Total mass would then be 94 kg, a savings of 42 kg. 
Furthermore, even if a third thruster were necessary to complete the 
mission (as a result of a variable thruster lifetime in ampere-hours), it 
would still only add 17 kg, and net savings would then be 25 kg. 

Savings in propulsion system mass thus can apparently be achieved 
by increasing beam current, permitting wearout to occur, and saving 
power processing mass. Even if the constant ampere-hour wearout assump- 
tion is not completely accurate the effect tends to be balanced by 
reliability considerations. A favorabla reliability interaction results 
at higher beam currents because the total thruster hours are reduced and, 
consequently, there will be fewer failures and fewer thrusters required 
for reliability purposes. Based on typical system studies, redundancies 
of 20% and 40% for the PPU and thrusters. respectively. were selected. 

Using the models in Table 3-1 and the previous philosphy for 
thruster wearout, propulsion system specific mass curves were developed. 

For the simplified example shown in Figure 3-7, mission thrust time is 
less than thruster lifetime (i.e., no thruster wearout occurs). The 
trends caused by increasing beam current and specific Impulse are then 
divorced from lifetime considerations. As is seen for present technology 
modifying the EMT from its nominal value of 2 A to around 6 A will generate 
significant decreases in specific Tnass. On the other hand, holding beam 
current constant and increasing the specific impulse from 3000 to 
4500 sec will accomplish approximately the same decrease in specific 
mass (10 kg/kW). Simultaneously increasing both beam current and specific 
impulse would only further improve specific mass by 5 kg/kW. The 
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Change le celatLely snail hecanse the hnlh of the benefU cones fron Che 

eichec vacurrsi^a:"!;.""’ 

“ 

^ho 4 . ® include wearout on 

the basis of a 20,000 A-hr thruster life A hr«.ie 4 

,, ire. A break point occurs at 2 A 

because additional thrusters are required ar - at 2 A 

maintain a . . required at higher current levels to 

maintain a constant mission life. The a j ^ 

nKfo4 j £ ‘ %s needed in Eq. 3-14 was 

obtained from Figure 3-8. ^ 

C. MISSION CHARACTERISTICS 

Referring to Eq. 3-u, the generalized analysis approach reonlrea 

att“ ^n»lal acceLraclon I ) 

nllr --h of the nlsslons Included In ^b’ 

laclor °*d ftO" the literature and fron trajectory calcu- 

a one node under this study, as Indicated In Tables 3-2 id 3-3 Pay- 

Thls dT “<=•■ "foofon are sho™ m Table 3-2 

This ^a es used CO select a naalnun payload nass for conpuclng 

dlsclni f Tn f«=onced In Table 3-3 L 

«ra A and°dv ^ '''^®'^®”‘^od literature. The paraneters of Interest, 

MFR la seed i tb° ! tcapectlvely. 

ased on the naxlnun payload values Indicated In Table 3-2. 

D. RESULTS OF THE GENERALIZED ANALYSIS 

and sonerallzed analysis nere derived fron Eq. 3-14 

e ase on the perfomance nodels, nass nodels, and nlsslon cb a t - 

istlcs discussed above. A general ezanple of the relatlonsbi t , 
is shown in Figure 3-9 Tho i j relationships involved 

n rigure j-9. The legend on top of Figure 3-9 dpsgir>r4u 

technology of the nartloiil at- 1 .u 8 re 3 9 describes the 

sldered ("Pri propulsion systen being con- 

in fabU liT ptariT-TT' 

»nre or present power)' la ehoin C-lstlng bard- 

Chosen. Then, for a modified EMT within the 


48 










Table >2. Mission Pay load Mass Summary and References (Continued) 
Mission 


Outer planets 
(Continued) 

Solar system 
escape 

Mercury orblter 


Close solar probe 
(0.1 AU) 


Earth observatory 

Comet Encke 
Slow flyby 


Payload 
Mass, kg 

References 

Remarks 

320 

8 

Hellos class 

1335 

30 

HAC study; Pioneer 
Venus Module (793 kg) 
Includes retro; drop 
propulsion system 
before Injection 

1480 

31 

Drops propulsion sys- 
tem before Injection 

824 

8, 29 

350 kg Helios 

1317 

51 

Mariner class 
(775 kg) plus retro 

465 

55 

Solar probe uses 

Venus swlngby 

1000 

55 

Low cost (heavy) 
payload 

8000-10,000 

11, 18, 21 

Shuttle class 
maximum payloads 

113 

8 

Small science 
package 

508 

56 

JPL Mariner Encke 
study; 4 km/sec 

358-608 

33 

Pioneer 10/11 class 
and Advanced Pioneer 

450 

13, 15 

Representative pay- 
load for flyby 
opportunity and back- 
ground for comet 
missions 
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Table 3-2. Mission Payload Mass Summary and Reference (Continued) 


Mission 


Comet Encke 
(Continued) 


Payload 
Mass, kg 


References 


Remarks 


Slow flyby 
(Continued) 

400 

38 


533 

25 


350 

32 

Rendezvous 

635 

41 


545 

51 


500 

38 

Asteroid rendezvous 

350 

8 


500 

34 


500 

35 


635 

34 


1 500 

38 


350 

32 

Asteroid sample 
return (Eros) 

940 

37 


1000 

38 

Out-of-ecllptic 
(1 AU) 

600 

11 


391 

51 


200-400 

2 



Mariner class 

HAC /Metis (Pioneer 
class) 

Eros 

Eros 

Ceres 

Various 

Various 

Return capsule, 
leave stage at Eros 


Includes science 
package of 31 kg 

Inclination dependeint 
payload 





I' »• •• 1 4^ 111* I **! ! i m ■! ^ 

ti . r* : iti,,, I. f ‘ ^ 4 Tiwuf 


ixansi^ SeaS 


Table 3-3. Mission Characteristics 


Mission 

C3 

km^/s’ec^ 

«o» kg 

Mpl 

Conservative 
Payload, kg 

MFR 

Aq. 

10-5g 

AV, 

km/ sec 

Refer- 

ence 

Remarks 

Outer planets 









Jupiter Orbiter 

37 

3050 

1360 

0.45 

3.4 

4.3 

8,29 


Jupiter 

satellite 

37 

3050 

1265 

0.41 

3.4 

4.3 

8,29 


arbiter 









Saturn flyby 

64 

2320 

840 

0.36 

3.4 

3.6 

8,29 


Saturn orbiter 

64 

2320 

1031 

0.44 

3.4 

3.6 

8,29 


Titan orbiter 

64 

2320 

1096 

0.47 

3.4 

3.6 

8,29 


Uranus flyby 

42 

2920 

825 

0.28 

3.8 

4.5 

8 

Jupiter swlngby 


72 

2120 

825 

0.39 

3.8 

3.8 

8 

Uranus flyby to 
Neptune 

Neptune flyby 

72 

2120 

750 

0.35 

3.8 

3.8 

8 

Uranus swlngby 

Solar system 
escape 

37 

3050 

450 

0.15 

3.4 

4.3 

8,29 

Jupiter swlngby 

Mercury orbiter 

4 

5720 

1480 

0.26 

2.1 

17.4 


HAC trajectory 
studies; thermal 
constraints, pay- 









load is approach 









mass 
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Table 3-3. Mission Characteristics (Continued) 


Mission 

^3, 

Ian2/sec2 

“o* ^8 

Mpl 

Conservative 
Payload, kg 

MPR 

Aq, 

10-5g 

AV. 

kia/sec 

1 Refer- 
ence 

Remarks 

Mercu^ orbiter 
(Continued) 

25 

3827 

1480 

0.39 

1.6 

10.9 


HAG trajectonr 
studies, thermal 
constraints, pay- 
load Is approach 



















mass 


25 

3827 

1480 

0.39 

1.7 

13.1 

31 


Close solar probe 
(0.1 AD) 

20 

4260 

1000 

0.23 

2.9 

14.1 


HAG trajectory 
studies; thermal 
constraints, (Sun 
angle > 90*) 

Earth observatory 

- 

27,215 

10,000 

0.37 


6.0 


Comet Encke 







HAG trajectory 
studies to sync, 
eq. 



Slow flyby 

55 

2560 

608 

0.24 

4. 1-4. 7 

8.8 


HAG studies; 
approach * 

4.5 km/sec 


54 

2600 

608 

0.23 

4-6 

9.5 

8 

4 Icm/sec 


77 

1970 

608 

0.31 

3.1 

11.2 

39 

4 km/sec 

— L 

49 

. 

2720 

608 

0.22 



2.7 


39 

Venus gravity 
assist 


6U9 










«^ iK ai tf'- r |, 


dttHifeff 

ir ^ 


’’’'-.'ijle 3-3. Mission Characteristics (Continued) 


Mission 

?3. 

km2/sec2 

«o» ^'-8 

“pi 

Conservative 
Payload, kg 

MFR 

lo-h 

AV, 

km/ sec 

Refer- 

ence 

Remarks 

Comet Encke 

55 

2560 

608 

0.24 

3.4 

11.3 

29 


(Continued) 

54 

2600 

608 

0.23 

2.5 

10.4 

2 


Slow flyby 
(Continued) 









Rendezvous 

49 

2720 

635 

0.23 

4. 4-4. 5 

9. 1-9. 8 

41 

Different 

opportunities 


54 

2600 

635 

0.24 

4.0 

11.1 

41 



42 

2920 

635 

0.22 

4.1 

9.5 

8 

Long flight time 


44 

2860 

635 

0.22 

4.2 

14.3 

8 

Short flight time 


51 

2680 

635 

0.24 

4.2 

10.4 

2 


Asteroid 

25 

3827 

635 

0.17 

3.3 

8.5 

8,29 

Metis 

rendezvous 

3 

5840 

635 

0.11 

3.6 

8.1 

34 

Eros 


4 

5720 

635 

0.11 

4.2 

5.9 

32 

Eros 


56 

2530 

635 

0.25 

4.6 

00 

34 

Ceres 

Asteroid sample 

4 

5720 

1000 

0.17 

4.2 

12.0 

36 

Return capsule. 

return (Eros) 








leave stage at 










Eros 
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Figure 3-9. Example of the Influence of I on 
AV vsA^ format results. 
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selected power technology, one must choose thruster operating conditions 
in terms of beam current and specific impulse. These choices will serve 
to describe the propulsion system according to the modeling in Sec- 
tions 3. A and 3.B. Mission requirements are defined by the payload 
(MFR) for mission success. 

1 . Specific Impulse Conslderatlgns 

The dependence of performance on is illustrated in Figure 3-9. 

Any mission can be categorized in terms of an initial acceleration A^ 
and a velocity AV required to be delivered by the propulsion system. 

The spectrum of A^ and AV combinations for all missions is shown in 
Figure 3-9 for a constant (MFR). The solid line on the curve corresponds 
to a particular propulsion system (system technology level, thruster 
operating conditions) that can deliver the specified MFR. To the left 
of these curves, the propulsion system performs the mission with an 
MF greater than that required; to the right, it performs the mission 
with the inadequate MF (or payload). The effect of increasing 
is shown by arrows. In the region of higher AV and lower A^, where cor- 
respondingly less thrust is needed, the preference is to keep the pro- 
pellant mass as small as possible, and increasing helps to do this. 

On the other hand, when higher thrust is needed at low AV, increasing 
I causes a decrease in performance. The reason is that the power 
required, in terms of mass, when increasing I^p is greater than the 
propellant savings that would result from the higher I^p- 

Figure 3-10 outlines the tradeoffs between I^p and various thruster 
modifications and MFRs. These charts also show equivalent thrust time, 
which is the mission duration parameter used to characterize the required 
thruster lifetime for any mission. This equivalent thrust time is 
actually the time It would take to expel all of the propellant with the 
power available at 1 AU. For orientation, the Halley’s comet mission, 
which is a relatively long thrusting mission, requires only about one 
year of equivalent thrust time. 

The trends, as far as I^p is concerned, show a very flat optimal. 
Considering Figure 3- 10 (a), for example, for very high accelerations 
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AV. km/sec 


7264-12111 


SPECIFIC IMPULSE, sec 



Ao.1tf®9 


Figure 3- 10(b). AV vs format for present power technology with 
modified EMT; MF - 0.2. 
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Figure 3- 10(c). 


AVvsAq format for 
modified EMT; MP • 


present power technology with 
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^V, km/sec 
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figure 3- 


AVvsAo format for 

modified EMT; MP a power technology with 
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(6 6 X 10 ^ g) and low AV (10 to 15 knt/sec), the performance resulting 

from an 1^^ up to 5000 sec Is relatively constant. The fact that there 
is not much difference indicates the very fiat optimal. At the extreme 
right (the envelope of maximum performance), little variation occurs In 
^sp ® great range of and AV combinations. The message of these 
curves is that, since increasing 1^^ may not significantly increase mass 
performance, the higher cost of higher power may not be justified. 

The other parts of Figure 3—10 have the same propulsion system with 
Increased MFs. Since the particular propulsion system now has all of 
the curves moved over to the left and down, fewer potential missions can 
be captured by the particular technology level and thruster modifications 
because of the higher MFR. The trends, as far as I is concerned, are 
basically the same* Since these missions are characterized by rather flat 
optlmals, the basic 1^^ considerations become a cost tradeoff: Increasing 
^sp improve performance, although sometimes possible, is expensive. 

In Bummsxy, the above considerations suggest that I should be 

sp 

raised only after other alternatives have been examined. It is generally 
better to first attempt to lower the specific mass of the propulsion sys- 
tem by raising beam current. 

2. Mission Performance Results 

The space shuttle was assumed available as a launch vehicle for 
these advanced missions (Table 3-3). The major Implication of the shuttle 
is that larger injected mass and larger propulsion systems are possible 
as compared to earlier launch vehicles (e.g.. Atlas Centaur). Fig- 
ure 3-11 shows the assumed capability of the Shuttle. The two-stage lUS 
is assumed for the lower launch energy missions and a three-stage con- 
cept for the higher values of 

The mission performance comparison curves used in the remainder of 

this analysis are shown in Figures 3-12 through 3-16. These figures 

display lines of constant mass fraction in the AVvsA format for 

o 

various technology assumptions. Given a particular mission, specified 
in terms of A^ and AV, any propulsion system (power technology level. 
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Figure 3-11. Interim upper stage 



6 ' 



Figure 3-12. AVvaA^ format; EMT, existing power subsystems. 


65 




8 13. AV vs Aq format; modified EMT (4 A) oreaeni- 

power technology, 3000 sec i ^ ’ Present 

sp* 
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aV, km/iec 



Figure J-16. AV vs format; modified EMT (4 A), present 

power technology, 6000 sec I . 

ap 
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Thia MF must be groeter than HFR it the pro? 

an adequate paylba^i* _ j „of<fo™ance of the existing 

»lg«. 3 - 12 ^«.e«ts ae^B»er. lot preooot 

rrtr^ - --- - 

"ZTtlTlr^iTr r^«» l^rrmote -isslons (») lot low A„ and 
but baa a low capability lor high and low AV. 

a, ntiter Planets 

itar.ai-« are depicted in Figure 3-17. 
mssiona to the outer planet^ „,„l,ementa superimpoaed. 

which la Figure 3-12 “ „f acceletatlona and velocitiea 

The email elliptical area ^ ahown ia a line repre- 

in Table 3-3 lor the outer ,,„n orbiter). 

aenting Che largeat MFR in a e ^ perlorm 

Obviously, since the BMT can deliver t e a information in Fig- 

the remaining rthruat and redundancy. 

„Pe 3-17 ia lur'becauac ol the abort miaaion eguivalent 

Thruster weatout does n n 

thrust time (EMT life was assumed to be 20,0 

b. Mercury Orbiter 

The results of the detailed trajectory studies of the Mercury 
. T.J 1 R The areas indicated on the 

orbiter mission are shown in Figure - • ^^w launch 

f-vim different trajectory classes investig ..... 

curves are the differ J ..^niilrements (MF - 0.26) and high 

energy with high propulsion system power ^ 0.39). The EMT 

Taunch energy with low propulaion pc,»r 

With existing hardware ^ of both classes of 

actual trajectory runs easily meets tne 

""^^\he generalized plots in Figures 3-12 and 

weeront oeeura lor the -^ 0 "/ Ition longer 

orbiter is one ol the l« mlaaiona that baa a thr 
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Figure 3-17. 


Outer planets mission results; EMT, existins 
power subsystems. * 
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Fl(;ure 3-18. Mercury orbiter mission results; EMT, existing 
power subsystems. 
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than 20.000 A-hr. Its equivalent thrust time la 16,000 hr (32,000 A-hr). 
Thus, four thrusters were Included In the maximum of 13 required for the 
high-power case to compensate for wearout (assuming a 20,000 A-hr life- 
time). The HAG trajectory studies of the mercury orbiter included 
thermal constraints which kept the thrust attitude 90» from the aunline. 
c* Close Solar Prnb<a 

The close solar probe (0.1 AU) was also analyzed using HAG 
trajectory programs with conservative MFs (Table 3-3). As. Figure 3-19 
shows, the EMT adequately performs this mission with a moderate number of 
rtirusters using existing hardware. Thlr. mission has a relatively long 
equivalent thrust time compared to other missions, but it is still 
only 10,400 hr, which is Just slightly more than the assumed 
20,000 A-hr lifetime of the EKT. The HAG trajectory studies again 

included thermal constraints that kept thrust attitude 90» from the 
sunline. 


Earth Observatory 

/n conservative payload mass fraction 

(0.37) required for the Earth observatory mission could be easily accom- 
plished with existing hardware, but the number of thrusters might be too 
arge for the Shuttle bay. To decrease the number of thrusters, present 
power technology could Le used with a modification of the beam current. 
For example, if beam current were Increased to 4 A, a reasonable number 
of thrusters for the very large payload results (as shown in Figure 3-21) 
Increasing the beam current also Increases the already adequate payload, 
ore performance and a further decrease in the number of thrusters could 
e ac leved by raising but this is not required to perform the Earth 

servatory mission. In addition, advanced power system technology is 
not required to perform this mission. 


Comet Encke Missions 

Th. comet Encke elom flyby mleslone are depicted by the ehaded 
eteas Indicated In Figure 3-22 (eee Table 3-3). Enlatlng herdware la 
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AV, km/sec 



Figure 3-19. Close solar probe (0.1 AU) mission results; EMT, 
existing power subsystems. 
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Figure 3-20. Earth observatory mission results; EMT, existing 
power subsystems. 
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Figure 3-22. Comet Encke slow flyby mission results; EMT 
existing power subsystems. 
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not adequate to perfom all plsalous but, aa Figure 3-23 sbous, all 
propoaed Couet Eucke alow flyby ulaaloua can be performed with adequate 
payload, using a small number of thrusters and a modified EWr. 

Similar results hold for the comet Encke rendezvous missions, as 
a own in Figures 3-24 and 3-25. Modifying the EMT to operate with a 4-A 
earn current with present power technology yields an adequate payload MF. 

Asteroid Rendezvnma 

The asteroid rendezvous mission requirements for Eros, Metis, 
and Ceres are plotted in Figure 3-26 from the data in Table 3-3 using the 
conservative payload of 635 kg. The existing hardware (EMT) easily meets 
the payload requirements. The 17 thrusters required is a reasonable 
ound for the moat conservative payload and acceleration but can be 
decreased, if desired, by increasing beam current. 

g* Asteroid Sample Ret-tim 

The asteroid sample return (Eros), with its somewhat higher 
payload and velocity requirements may need a modified thruster to obtain 
additional performance depending on launch opportunity (as Indicated in 

gures 3-27 and 3-28). Changing the beam current to 4 A provides this 
extra performance with present power technology. 

h. Out of Ecliptic 

Figure 3-29 depicts the performance for the existing hardware 
torn the HAC trajectory runs as a function of the out-of-ecllptlc Incll- 
pstion angle. Existing hardware can provide up to 50- with a payload mass 
ractlon of 0.1. Raising the beam current to 4 A Increases maximum Incll- 
^tion to about 60-. In addition to using a 4 A thruster, it is necessary 
o Increase 1^^ to obtain a 90* inclination because the very high velocity 
requirements for 90” inclination can only be met by reducing propellant 
^ss ^ Trajectory results for a 6000-sec are given in Figure 3-30. 

90 Inclination is required, then I and cost will be higher. 

o ecliptic missions can also have very long equivalent thrust 
imes Acceleration increases as propellant is consumed, but the average 
time is about 11,000 hr for each 30" of inclination. Thus, the 
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Figure 3-23. Comet Encke slow flyby mission results; modified 
EMT (4 A, 3000 sec), present power technology. 
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Figure 3-24. Comet Encke rendezvous mission results; EMT, 
existing power subsystems. 
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Igure 3-25. ^^J/ncke rendezvous mission results; modified 
EMT (4 A, 3000 sec), present power technology. 
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gure 3-27. ERM sample return mission results; EMT, existing 
power subsystems. ** 
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Figure 3-28. EROS sample return mission results; modified 
EMT, present power technology. 
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Figure 3-29 


1 AU out-of— the-ecllptlc mission results; EMT 
existing power subsystems. 
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Figure 3-30. 1 AU out-of-the-ecliptlc mission results; modified 

EMT, (4 A, 6000 sec) present pover technology. 
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thrusters required (in Figures 3-29 and 3-30) are mainly for wearout 
replacement. In fact, the modified EMT (4 A, 3000 sec) requires only 
four operating thrusters to provide the highest acceleration. Thus, for 
the highest inclination case (90°), eight thrusters are required for 
wearout and another eight for failure for the conservative payload and 
maximum acceleration. The equivalent thrust time compensation provided 
in the propulsion system assumed in Figure 3-30 results in an increased 
mass fraction of less than 0.05 for the 6040 kg injected mass in this 
worst case. 

3- Conclusions from the Generalized 'Vnalysis 

Tn summary, the EMT can perform a large fraction of the missions 
analyzed. Development of a thruster (modified EMT) capable of 4-A or 
higher beam current at 3000-sec specific impulse \ 7 ill provide the 
capability to perform all the missions considered, with the possible 
exception of the out-of-ecliptic mission. Increasing I is generally 
good from a performance viewpoint because it improves the payload and 
reduces the number of thrusters required, but it could increase mission 
cost because of higher power. Similarly, the 90° out of ecliptic mission 
can be accomplished with present power technology (solar array and PPU) 
by increasing specific impulse. Advanced lightweight power technology 
would improve performance, but it is not required to perform any of these 
missions and may not be economically justifiable. 

This analysis did not consider Earth orbit missions Involving large 
space systems. This subject is discussed in Section 4. In addition, the 
generalized analysis did not consider other thruster sizes since the 
existing EMT or a modified EMT could adequately accomplish all the 
missions considered. 
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SECTION 4 

COST MODELING FOR EARTH ORBIT MISSIONS 


The generalized analysis approach, discussed in Section 2, is an 
efficient means for assessing the first-order impact of various missions 
on thruster technology. However, the generalized approach omits cost, 
which, under certain circumstances, may be decisive in the use of EP. 
Although cost is an important factor in planetary and small body 
missions, for these the use of EP will probably rest on the tradeoff of 
mission risk against technical advantage, with cost as a secondary 
factor. 

With new large space systems, for which most elements will have 
limited bases for assessing reliability/risk, EP will be judged largely 
on its merits as a cost-effective propulsion system. This section 
describes cost models developed for several Earth orbit missions to 
indicate where ion thruster technology would best satisfy selected cost- 
effectiveness criteria. In general, the goal is the technology that 
would produce the lowest mission cost consistent with the reality of 
developing such technology. In other words, from a mission standpoint, 
the goal is to determine the thruster (and possibly other) parameters 
that would minimize mission cost, and then, through an understanding of 
technology limits, recommend a direction that has a reasonable proba- 
bility of success. 

Two general types of Earth orbit missions were considered: 

(1) orbit raising (transportation) of large masses from low (altitude) 
Earth orbit (LEO) to higher altitude orbit and (2) stationkeeping and 
attitude control of large masses in Earth orbit. Details of the assumed 
missions are discussed in "mission description" sections. The computer 
mjdels developed are presented in the "general model" sections, which 
include descriptions of the mission analysis portions of the models. 
Propulsion system performance models, which apply to both types of 
missions, are discussed in the orbit raising mission section. Similarly, 
power source degradation and redundancy models, which also apply to 
both types of missions, are discussed in a single section. Mass and 
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cost models are discussed separately for each type of mission because 
the models and terminology vary somewhat. Results of each model are 
grouped into two categories: (1) baseline results using model parameter 
values selected as being reasonable compromises between 1977 technology 
and future technology and (2) sensitivity results obtained by varying 
each major parameter over a wide range of values. Certain combination 
sensitivities, obtained by varying two or more parameters, are also 
presented in the latter category. 

The basic graphical format presents cost as a function of specific 
impulse. This format was selected because I is one of the most 
fundamental parameters in thruster technology. For orbit raising mis- 
sions, the cost scale is in terms of cost per unit mass of net payload 
delivered to the destination (i.e., $/kg of payload). For on-orbit 
missions, the cost scale is in terms of cost per unit mass of net satel- 
lite maintained on orbit per year ($/kg/yr). Other plots are also pre- 
sented to aid in understanding the results and the various relationships. 


A. ORBIT RAISING (TRANSPORTATION) MISSIONS 

To carry out the economics analysis with the cost modeling program 
described in this section, a relatively large number of assumptions and 
detailed models were required. The orbit raising mission work was 
generally oriented toward the transportation of masses from LEO to GEO. 
However, this particular transfer is only representative, and other 
transfers can be simply handled by the model through the specification 
of AV. To further clarify assumptions, two baselines were considered: 

(1) matured Shuttle era and (2) large space systems era. Such a dis- 
tinction was made to identify differences in thruster technology require- 
ments that might result from different launch vehicle sizes and levels of 
activity in space. For instance, the matured Shuttle era baseline assumes 
a payload size of 25,000 kg, a launch cost to LEO of $300/kg, solar power 
cost of $100/W, and a solar power specific mass of 6 kg/kW. The magni- 
tude of these factors significantly influences the results of the cost 
analysis. 
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For the Iprge space systems baseline, it is assumed that the 
applications being considered are in the area o£ highly expanded use of 
space, the masses of interest are larger than those of current-day satel- 
lites, launch costs are lower, and power sources are lighter and cheaper. 
Specifically, masses in the range of 10^ to 10^ kg are considered. Such 

masses would be delivered to LEO by new launch vehicles (e.g. , a heavy 

18—21 

lift launch vehicle) at relatively low cost ($10 to 50/kg). Power 

sources are projected to weigh 3 kg/kW at a cost of $0.5/W. Other dif- 
ferences in these two baselines are discussed in conjunction with the 
results. 

Once in LEO, the payload must be placed in another orbit. An 
electrostatic ion propulsion system, using a solar source, transports 
the mass to the specified orbit. Except for characterizing this trans- 
fer by AV, the details of the transfer were not considered in this study. 
Clearly, many factors must be considered in designing such a transfer 
from LEO (e.g., occultation, attitude control, solar array orientation 
and vehicle maneuvers). However, the impact of most of these factors 
on propulsion system design is through AV. The mission options considered 
(such as one-way or round trips) are discussed in Section 4.A.I. 

For the orbit raising mission, the propulsion system is part of an 
orbit raising vehicle (ORV) , shown schematically in Figure 4-1, which 
contains the subsystems (e.g., computers, data handling, communications) 
necessary to operate with only limited ground interfacing. Models of 
the various elements of the ORV used in this study are described in 
Sections 4. A. 6 and 4. A. 7. As mentioned earlier, the analysis of the 
sensitivity of the cost results to the various models was a significant 
part of the work. 

1. Mission Descriptions 

The three types of orbit raising missions considered are shown 
schematically in Figure 4-2; they include; 

• Self-powered round trip — ORV contains the power source to 
supply the necessary propulsion power. On the "down trip," 
the full power is used. Power source mass and cost are 
charged to the ORV. 
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Figure 4-2. Orbit raising mission scenarios. 
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Pavload nowered up; salf-pouared douu - payload Is, or contains, 
I pow7r sourL cap;i,la of supplying tha propulsion poyar dur- 
ing tha ur^rip. On tha down trip, the ORV provides Its o^ 
power. Up trip power source mass and/or cost may e - ai g 

to payload. 


• One-way trip (up) “ ORV contains the power source, but stays 
with the payload on-orbit. 

Although many other mission options can be defined, theoe three should 
be representative of a wide range of mission requirements. 


a. Self-Powered Mode 

In this mode, the ORV may make a given number of round trips 
from LEO to a destination orbit (defined by AV) consistent with the 
specified ORV life. Propulsion power for both the up and down trips is 
provided by the ORV, and the power source mass and cost are charged against 
the ORV. The number of thrusters required per trip is based on the number 
operating, wearout rate, and redundancy for random failures. 

Assuming that the use of the ORV for a given transfer is paid 
for by a "user," the fraction of the worn out thrusters must be taken 
into account. This is done by subtracting from the cost to the user 
the fraction of unused thruster life and adding the cost for refurbish- 
ment. In this way, each user effectively on]y pays for the fraction of 
the ORV worn out. This approach also eliminates the need to transport 
a large number of spare thrusters through many round trips before they 
are needed. In addition, propellant for only one round trip is carried 
for the same reason. More specific details of this mode are discussed 

in the "general model" section. 


b, Pavload Powered Up-Trip : Self-Powered Down-Trip Mo^ 

Since the ORV mass is generally only a fraction of the payload 
mass, and down-trip time may be less critical than up-trip time, signifi- 
cantly less power may be needed for the down-trip. In addition, many 
future payloads may contain large power sources that could be used on the 
up-trip. For such missions, the cost charged to a user might be reduceo 
since a smaller power source would be paid for as part of the transporta- 
tion cost, and the net payload delivered on-orbit is increased. 
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c. One~Wav Trip (Up) 

In certain cases, it might be advantageous to leave the ^rv 
with the payload on orbit. For this assumed mode, the power source is 
charged to the ORV. The cost results will be conservative if the pay- 
load can make use of this power on orbit. Propellant is sized only for 
the up-trip, and operation of the ORV on orbit would require using part 
of the payload for propellant. 

2. General Model; Orbit Raising 

A logic diagram for the orbit-raising cost model program is shown 
in Figure 4-3. This diagram correctly indicates the functional form of 
the program, but does not necessarily correspond identically with the 
computer program listing. Efficient coding and the many details not 
shown in this general diagram introduce format variations. A typical 
program listing and a discussion of the listing are included in 
Appendix C. 

Mission mode selection indicated in Figure 4-3 is actually accom- 
plished by selecting a few parameters within the program (see the listing 
discussion in Appendix C) . For a given set of mission variables, the 
mission computations are started. Propellant mass is obtained from the 
rocket equation (see Eqs. 3-1 and 3-2). Although this equation is not 
strictly accurate for low-thrust propulsion, conservatism in selecting 
AV tends to compensate. 

Propellant mass can be related to ion beam "jet power" as follows. 
The ion beam energy is 


ev^ . Ej . i , 


or, in terms of beam power. 


P _ 1 * 2 

- 2 Vi 


(4-1) 


(4-2) 
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Figure 4-3. Orbit raising cost model logic diagram. 
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Figure 4-3 


Continued 













where 




ion beam kinetic energy 



rate derivative of E 

j 


= 

total jet power 

''i 


ion velocity 

”l 

=: 

ion mass 

”1 

= 

total ion flowrate. 


Total propellant mass flowrate Is related to total ion flowrate 


by 


m. = M q 

1 P u 

(4-3) 

= M t- , 

P P f 

(4-4) 

where t^ is the flight time appropriate to the 
(self-powered, one-way, etc). 

mission being considered 

_ ^8p®0 

V. = — » 

1 n Y 

u' 

(4-5) 

where r\^ and y are as defined in Section 3, A. 2. 
4-4, and 4-5 into 4-2 yields 

Substituting Eqs, 4-3, 

M (I g 

P = - P sp^o^' 
j 2t^ (q^y2) 

(4-6) 


To simplify the jet power computations, the quantity n was assumed 
equal to one. This approximation, at this point in the calculation, 
introduces an error of 5 to 10 percent in Pj . However, since the rocket 
equation for computing Mp is approximate, AV can be adjusted for con- 
servatism, and thruster and PPU redundancy can be adjusted, the approxi- 
mation will no*- significantly affect the results. (P^ is used only in 

computing the number of operating modules; P is divided by P where 

P = T V ^ 

mod -^b b*'' 


98 



The next step in Figure 4-3 is to calculate total propulsion system 
input power according to 


P 

ps 


M (g I )' 
p o sp 

2t- n 
f ps 


(4-7) 


where r| is as defined in Eq. 3-8. With the use of n in Eq. 4-7, the 
ps ps 

approximation made in computing P. is not included in P since n 

23 ps ps 

includes n^Y as well as other efficiencies. Total propulsion system 

power is used for sizing the solar power source, as explained in 

Section 4, A, 5. 

The total number of modules needed for the ORV for one trip Includes 
those operating, redundant modules for random failures, and spares for 


wearout if required: 




4J 

o 

= N , 
mod 

/totrip\ \ 

1 ^ ) r^) 

(4-8) 

N 

op 

= N , 
mod 

(l + Rp) , 



where 


N 


ot 


N 


mod 
N 

op 

^raod 

TOTRIP 

R 


P./P , 
j mod 


total number of thrusters installed in ORV 

number of operating thrusters 

total number of power processors 

selected module power (beam power) 

selected up-trip time plus down-trip time 
thruster redundancy 

PPU redundancy 

thruster life. 


If TOTRIP is < then TOTRIP /L^. = 1. PPU life is assumed to be long 
compared with trip time, and wearout can be neglected. Only enough 
thrusters are carried to support one round trip. Those used during the 
mission are assumed lo be replaced. 
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The round-trip mode' s require an iterative procedure such as 
illustrated in Figure 4-3. Initially, the down-trip time is assumed 
to be half the up-trip time. Since the down-trip time depends on the 
down propellant, which is not known until is found, the calculations 

proceed until the correct down-trip times are found (typically three 
iterations). The net payload is obtained by subtracting the down pro- 
pellant as indicated in Figure 4-3. M-ss models used to define the ORV 
are discussed in Section 4. A. 7. 

Once a solution to the trip time and payload mass calculations is 

obtained, cost models are applied to determine the total cost. Details 

of the cost models are presented in Section 4. A. 6. Total cost divided 

by net payload mass is the basic answer sought for each data point. 

Generally, this calculation cycle is repeated for a range of I ’s and 

for a series of values of a selected parameter (e.g., P M .^flleht 
^ mod o ® 

time) . 

Since down propellant is not used for a one-way up-trip, the itera- 
tive procedure is not required, and net payload is calculated directly. 

The other difference for one-way trips is the lack of amortization 
of the OR/. For multiple trips, a user is charged for only a fraction 
of the ORV life. However, for a one-way trip, the total cost is charged 
for the one trip. The specifics of the cost allocation are given in 
Section 4. A. 7. 

^ • Propulsion System Performance 

The various general equations defining the performance of an EP 
system are presented in Section 3.2. The purpose of this section la to 
define the assumptions used in simplifying and applying the general 
equations to the Earth orbit cost modeling analysis. Since the cost 
modeling includes several approximations, simplifications of the per- 
formance models are believed to be consistent with the overall modeling 
accuracy. 

Propulsion system total efficiency (Eq. 3-23) can be written in 
the form 


ps 


= n n n n y 

u e pp 'c ' 


(4-9) 
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In the cost model analysis, the following assumptions were made to 

59 50 

simplif> the computations: ’ 


n n 

PP c 


(4-10) 


'e V. + 220 
b 

These values are representative of highly developed components (e.g., the 
30-cm Hg thruster and FPU) as might be expected from a development pro- 
gram. These values are assumed to apply to all propellants considered 

(Hg, Xe, Kr, Ar), (This assumption will be shown to be justified by the 
lack of sensitivity of cost to minor variations in performance.) Using 
the values defined by Eq. 4-10, Eq. 4-9 becomes 


= 0.80 


V, + 220 
□ 


(4-11) 


Total propulsion system input power, Eq. 4-7, is determined using this 
relationship. The relationship between and 1^^, Eq. 3-35, can now 
be simplified to 


1 > 
T, I ^P 

\ 90 


(4-12) 


. (Hg) 


(k)' 


I 

140/’ 


\200/’ 


(4-13) 


These relationships (solid lines) are plotted in Figure 4-4. Also shown 
are lines of constant thruster diameter for a module power of 100 kW. 
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Figure 4-4 
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^sp \ for various propellants. 





To determine thruster size for given beam current and beam voltage 
requirements, the considerations illustrated in Figure 4-5 must be 
addressed. Perveance, the proportionality between beam current and beam 
voltage, is discussed in Section 5, A. The result of that technology 
assessment was used in the cost modeling analysis in the form 

I = 8.19 X 10"® V , (4-14) 

b b 


where 

D « thruster beam diameter, cm 
X = ratio of atomic mass of propellant to Hg, 


or, normalized to the 30-cm thruster. 


I 


b 


2.33 

2.33 

2.88 

3.60 

5.22 



Vioooy 

3/2 

1 

x-1/2 


\1000^ 

3/2 
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ilof 1 

f 

uooo; 

3/2 
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, (Xe) 
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3/2 
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\iooo; 
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, (Kr) 
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/ V. N 

3/2 


i^o) 

Vioooj 

1 . 

, (A) . 


(4-15) 


As shown in Section 5. A, Eq. 4-14 assumes a net-to-total accelerating 
voltage ratio, "R", of 0.7. 

Thermal limits are established by thruster maximum operating tem- 
peratures. Results of the thermal ‘analysis, presented in Section 5.D, 
show that the thermal limit can be expressed in the form 


db) - 15 ( 20 ) 

max 


(4-16) 
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Figure 4-5. Thruster operating limits. 
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Thus, in Che process of calculating the size of modules required to 

use the selected power, Eq. 4-15 is used if the current is less than the 

thermal-limited value. For example, the cost model calculations select 

a module power and an I for each data point. Thus, I, and V, are 

sp b b 

established. If the beam current found from Eq. 4-15 satisfies the 
inequality 



then the perveance equation is used to determine beam diameter. If the 

inequality is not satisfied, thruster diameter is determined from the 

thermal limit equation. Thruster diameter as a function of I and 

sp 

module power is shown for argon in Figure 4-6. 

The maximum voltage limit indicated in Figure 4-5 is of practical 
significance in developing a thruster. As shown in Figure 4-4, the 
maximum beam voltages are below 7 kV for thruster diameters larger than 
30 cm. In terms of breakdown, this voltage level may be slightly 
inconsistent with the close grid spacing assumptions used to develop 
Eq. 4-14. However, moat of the calculations for the cost model Involve 
lower voltages and should not produce breakdown concerns. 

The lower limit shown in Figure 4-5 is of little practical impor- 
tance here since maximum power per module is of interest. At very low 
beam currents, thruster stability must be considered. 

Thruster life is affected by several operating conditions (e.g., 
beam current density, discharge voltage, propellant utilization efficiency). 
However, in developing a new thruster for a specific application, a given 

life could probably be achieved by selecting operating conditions. 

61 6A * 

materials, and designs. ’ Conversely, present knowledge of thruster 

wearout life could be extrapolated to larger diameters by normalizing 

to the 30-cm beam current density. For cost model program purposes, 

two thruster life models that assume either constant life or a life 
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proportional to beam current density have been formulated and described 
in the relationship 


where 


10 


(i} 


(4-17) 


w = (jj “ 4 or 5 for constant life 
c 

or 

u) = « 4.4 - 1.27 log (I^^) 

where 

I = equivalent 30-cm thruster beam current 
eq 



The model for was based on 30-cm thruster life test results. A third 
life model that assumes constant ampere-hour life (as used in Section 3) 
would predict lives intermediate between the constant life and the beam 
current dependent models. 


4. Redundancy 

In a typical set of cost model calculations (e.g. , cost versus I 

sp 

as a function of module power), the number of thruster and PPU modules 
varies from point to point. To be consistent, each point in this map 
should have the same reliability. Since reliability calculations are 
somewhat more complicated than simply adding redundancy, an analysis was 
performed to assess the error introduced by using a fixed redundancy. 
Details of that analysis are included in Appendix D. The basic result 
of the analysis is illustrated in Figure 4-7. 

For reasonable assumptions on failure rate, standby redundancy, and 
operating time, reliability becomes relatively high for large numbers of 
modules. Thus, if a redundancy is selected that makes the reliability 
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Figure 4-7. Reliability/redundancy results. 
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high for any number of modules, little error will be Introduced by 
assuming fixed redundancy. Based on these results, a redundancy factor 
of 20 percent was used for both PPUs and thrusters. Other results shown 
in Appendix D illustrate that, even for small numbers of modules (e.g., 
25) and lower redundancy for which the reliability would be reduced, the 
final transportation cost is not significantly affected. 


5. Degradation 

The mission approach used in the orbit-raising mission cost modeling 
assumes low thrust propulsion through the Van Allen Belt. The degrada- 
tion model used in sizing the solar power source is shown in Figure 4-8. 
The model was based on published degradation factors used for silicon 
P/N solar cells (curve A), Typically, low-thrust missions are expected 
to incur degradations in solar array outj.ut on the order of 25 percent to 
50 percent depending on the trajectory. 

Since future solar cell or design Improvements may reduce degrada- 
tion, two additional curves (B and C) were added to Figure 4-8 to reflect 

10 percent degradation in array output with equivalent 1 MeV electron 
15 16 2 

fluxes of 10 and 10 cm , respectively. This model is also included 
in the on-orbit calculations since the long operating times may accumu- 
late large dosages even though the flux rate is low. 

The power degradation factor used in the model is viewed as an 
efficiency; 


_ final power 
pd initial power 

Thus, the power used for sizing the power source is 


’’tot ■ 


(4-18) 


6. Mass Models 


Mass models for the propulsion system and ORV, as defined in Fig- 
ure 4-1, are summarized in Table 4-1. In the cost model program (see 
Figure 4-3) , net payload is obtained by subtracting the ORV mass from 
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Table 4-1. Ilass Models Used in Orbit- Raising Cost Model 



the initial mass. Except for determining total propellant mass, the 
calculations do not require iteration. 


a. Thruster Mass 

The mass assigned to each thruster in the system model makes 
an allowance for gimbals and miscellaneous structures by d'/ubling the 
thruster individual mass. The model for the individual thrueter mass 
as a function of diameter is shown in Figure 4— 9« A resi^i't able fit to 
the data points is provided by 

M^ = 0.078 (D in cji) . 


Ill 




THRUSTER BEAM DtAMETER, CM 


7264-41 HI 



THRUSTER MASS. KG 


Figure 4-9. Thruster mass versus diameter. 
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Thruster total mass Is estimated by assuming Items such as gimbals and 
thruster "array" structure double the thruster effective mass. 


M ^ = 2(0.078 N ^ . 

rt ot 


(4-19) 


b. PPU Maas 

The historical curve for total PPU mass plotted In Figure 3-4 
was converted to specific mass, as shown in Figure 4-10. This curve is 
represented by the equation 

ttpp = 10^“°’^°^ (kg/kW) . (4-20) 

where 

Ppp = power processor input power (kW) . 

Since the model is only approximate, the computer program assumes 
module power, is equal to Ppp for convenience. These powers are 

actually related by 


^raod " ’^PP^e^PP » 


(4-21) 


l^ut Ppp is not calculated in the program. This results in a more con- 
servative estimate for ttpp than would be obtained if Ppp were actually 
used in Eq. 4-20. Power processor total mass is then obtained from 


10"^ a__ P^^ (1 + R ) , 

PP ps p 


(4-22) 


where a„_ is in kilograms per kilowatt, P is in watts, and R is the 
PP ps p 

PPU redundancy factor. 
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c . Propellant Tankage 

Propellant tankage, including valves, lines, and miscellaneous 


hardware, is modeled as a fixed fraction of the propellant mass; 


M 


rpt 


= a 


M 

rpc orp 


(4-23) 


where 


M = M 
orp o 


1 - exp \- 


(-rt) 


- tankage fraction 

The baseline cases use = 0.1. Sensitivity studies show that cost 


results are not sensitive to the reasonable variations in a 
d. Power Source 


rpt* 


Power source mass is typically considered in terras of specific 
mass (kilogram? per kilowatt). Thus, 


M = lO"^ a P 
tw rw tot 


(4-24) 


where 


= total installed solar power (W) 

= power source specific mass (kg/kW) . 

For the matured Shuttle era and large space systems era baselines, a's 
of 6 kg/kW and 3 kg/kW, respectively, were assumed. 

e. Subsystems 

The model used for ORV subsystem mass is 


M 


rss 


100 


(1 + 10"^ P^ J + 
tot 


fe-‘) 


(4-25) 


where the first term in brackets represents the dependence on power 
level, and the second term accounts for complexity added by PPU/thruster 
switching. In the baseline cases, equals 1.0 because total trip 

time is less than the thruster lifetime, and no switching is assumed. 
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f. 


Structure and Meehan 1 stns 


The mass of ORV structure and mechanisms is modeled as a 
fraction of the sum of the other ORV component masses, including 
propellant: 


M = a (m 
rsm sm \ rt 

The baseline value of 0.2 

7 . Cost Models 


\pp “orp ^rpt \ss + ■ (4-26) 

was assigned to throughout this study. 


Cost models for the ORV are summarized i: Table 4-2. The overall 
objective of these cost equations is to obtain total orbit-raising 
cost, C^^. This total cost divided by the net delivered payload H 
is the basic final product of each calculation. 

The total cost to launch and transport a given payload mass to 
orbit is expressed in the form 


or 


%ef ^1 + ^rt ^opf %t S + C^rp > 


where 


C4-Z7) 


^oef ~ utilization cost factor 
= ORV hardware cost 
~ cost of thrusters 

^opf “ flight operations 

‘^opt " time-dependent cost of flight operations 
C^ = launch cost to LEO 
^orp " propellant • 

Each of these terms is defined quantitatively in Table 4-2, and each is 
discussed below. Many of the subscripts are defined In Figure 4-1. 


r 

! 


Table 4-2. Cost Models Used in Orbit-Raising Cost Model 
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QRV Utilization Factor 


The ORV is considered to be a Leasable vehicle that is rented 
to a "user" for payload transport. A user utilization factor has been 

included in the form 


C 


oef 



+ r TOTRIP 


(A-28) 


where 

TOTRIP 

L 

orv 

r 


flight time per trip 
ORV useful life 
interest rate 


Part of C j represents the fraction of ORV useful life (Lgry) charged 
to the user. The second part represents interest on the cost of the 
ORV during the flight. An ORV life of 2000 days and an interest rate 
of 10 percent per year were used throughout this study. 


b. ORV Hardware Cost (C^^) 

This factor includes recurring costs and development costs 


amortized over a given fleet of ORVs. 


— r 4-r +C +C +C (A 29) 

^rp rpt rsm rss rw ori ord 


where 


C = cost of PPUs 

rp 

C = cost of propellant tankage 

rpt 

C = cost of structure and mechanisms 
rsm 

C = cost of subsystems 

rss 

C = cost of power source 

rw 

C = cost of ORV integration and testing 
ori 

C = cost of design,' development, test, and evaluation (DDT&E) 

ord 
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Thrusters, propellant, md operations are not included in because they 
are the components that are worn or fully used each trip. These costs 
are calculated separately. 


c. Thruster Cost 

The relatively large systems being considered require a large 
number of thruster and PPU modules. The thruster cost model assumes 
that a mass-production type fabrication approach will be used to reduce 
costs. Furthermore, thruster costs are modeled by a learning curve, in 
which the unit cost decreases as quantity increases in the form 



where 



/ T0TRIP\ 

\ / 


first unit cost 

learning curve parameter 
thruster life (hours) 



The last term in Eq. 4-30 is the fraction of thrusters worn out although 
this wear is distributed among many thrusters. A user only pays for 
those used. The factor of 2 in Eq. 4-30 is included to account for gim- 
bals, structure, and other miscellaneous costs associated with the thrus- 
ters. All results were obtained using an 80 percent learning curve, 
giving a = 0.34. This means that the unit cost is reduced by 20 percent 
each time the quantity doubles. For instance, if the 100th unit costs 
$100, the 200th unit will cost $80. Since the summation in Eq. 4-29 is 
inconvenient, an approximate equation \jas derived in the form 


N 


rt 


= 2C’ 


2a / 

ot (TOTRIPj 


rt ,a \" L 


(4-31) 


119 



Thruster first unit cost was obtained from Figure 4—11, which displays 
current thruster costs and a projection for future costs. The rationale 
used for reducing cosLs is th^t future systems for Earth orbit missions 
will allow reduced fabrication tolerances and greater performance dis- 
persions. The close tolerances of present thrusters contribute signifi- 
cantly to cost. Using the equation that fits the projected thruster 
cost in Figure 4-11, Eq. 4-31 becomes 


rt 


= 2 (4.2 X 10^ l >g D) 


(1.47 


( T0TRIP \ 

^ ) 


(4-32) 


where D is in cm, 

d. Power Processing Unit Cost (C^p) 

PPU cost is modeled similarly to the thruster using an 80 per- 
cent learning curve. The first unit cost was derived by fitting a curve 
to cost data for 8-cm and 30-cm PPUs. The PPU cost is then expressed as 

= 1.88 X 10^ 

where 

P = power processor input power. 

PP 


e. 


mass: 


Propellant Cost (C^^p) 

Propellant costs are assumed to be proportional to propellant 


C 


orp 


C’ M 
orp orp 


(4-34) 


where 


C' = cost per unit mass 
orp 

= $0.4/kg for baseline • 


f. Propellant Tankage Cost (Cj.pj.) 

Tankage costs are assumed to be proportional to tankage mass 
using Eq. 4-23 
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where 


(4-35) 


C . = C' a M 

rpt rpt pt orp 


^rpt * tankage cost per unit mass 

= $100/kg for baseline 

= tankage factor 

= 0.1 for baseline • 

S' Structure and Mechanisms C ost (c ) 

^ ~ rsm'^ 

The structure and mechanisms cost model was derived from 
consideration of past vehicle costs and has the form 

^rsm ” ^rsm ^\sm^ , (4-36) 

where 

*"rsm baselines. 

Since the details and complexity of these elements would affect their 
costs, any general model necessarily must be rather arbitrary. However, 
using the model, the sensitivity of overall transportation cost to 
variations in structure and mechanisms cost can be evaluated. 

h. Subsystems Coat (C ^ 

^ rss 

As with structure and mechanisms, the ORV subsystem cost will 
be highly dependent on design details. Again, through discussions of 
various programs, the following model was assumed: 


where 

C 

rss 


C 

rss 


C' (M ) 
rss ras*^ 


0.9 


10 $/kg for both baselines. 


(4-37) 


compared wUh E,. 4 - 35 . the difference In the exponents reflecte the 
difference. In complexity. If the sire of . etructur. is Increaeed. the 
eost Is probably not too strongly affected. However, with electronice, 
increased mass certainly means Increased cost. 
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i. Power Source Cost fC > 

— rw 


Power costs are modeled in the conventional form 


C = C P 
rw rw tot 


(4-38) 


where 


- cost per unit of power 
= $0.5/W (large system baseline) 

= $100/W (Shuttle era baseline) 

Tne $100/W value is only a slight reduction from present technology and 
is consistent with goals for solar power source development. The $0.5/W 
is about three orders of magnitude less than current technology. How- 
ever, such low power costs will probably be necessary to ma'ke a large 
space system economically viabl-? (independent of SEP). 

j. Integration and Test Cost fC I 

^ 

This factor is intended to account for ORV recurring costs for 
assembly and testing. It is assumed that a fleet of ORVs would be 
produced in a time scrle consistent with production-line assembly. The 
cost of such assembly is modeled in the form 


ori ^ori (^^orv ^orp "^used) ’ (4-39) 

where 


C , 
ori 


used 


cost of integration per unit of ORV dry mass 


mass of thrusters worn out during one round trip 


M. 



By including factors for thruster wearout, the user pays only for the 
thrusters used (see Eq. 4-30), but must pay for the future Installation 
of new thrusters. 
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k. 


D esign, Development, Test, and Evaluation 


Design, development, test, and evaluation (DDT&E) costs are 
Included In the model for completeness In the sensitivity analysis. The 
form of the relationship, 

C’ 

ord Nr \ orv orp/ 

Is similar to that used for structure and mechanisms. Using the rela- 
tively small exponent seems appropriate since should depend on mass 

and complexity, but significant non-engineering costs should be largely 
independent of mass. The other parameters in Eq. 4-40 are defined as 
follows: 

C' , = DDT&E cost coefficient 

ord 

= $10^ $/kg for both baselines 

N = number of ORVs over which DDT&E costs should be amortized 
r 

= 10 for both baselines 


1. Flight Operations Costs (C , C . ) 

^ opp opt 

These cost factors arise from the philosophy that an SEP ORV 
flight would incur both fixed costs (similar to those for a ballistic 
transfer) and time-dependent costs due to the multi-month duration of 
the flight. For both factors, it is assumed that only the efforts 
directly related to the ORV (i.e., not the payload or launch vehicle) 
should be charged to ORV cost. Thus, the following values were used 
throughout; 


^opf 

^opt 




TOTRIP • 


(4-41) 


m. Launch Cost (Cj^) 

This element accounts for the coat of launching the initial 
mass from Earth to LEO and thruster refurbishment. The model for has 
the form 

■ ‘=1 »o * ^used) • 
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where 


C 


I 


cost per unit mass 

$300/kg for Shuttle era baseline 


= $15/kg for large system baseline 

I 

For the multiple round trip mode, the results are conservative because 
the full launch cost was charged to each launch. The $300/kg value is 
based on a possible reduction in current Shuttle projections of $500 to 
$800/kg. The $15/kg is based on launch cost estimates developed for 
heavy lift launch vehicles. Overall cost sensitivity studies for each 
of these quantities are presented in Section A. A. 8, 

8. Orbit Raising Laree System Baseline Missi on Results 

Cost model program results are structured around the two baselines 
selected to represent near future applications (matured Shuttle era) and 
late 20th century applications (large space systems era). For each 
baseline, a set of parameters required by the computer model was 
selected. Computer runs using variations around this baseline set of 
parameters were then performed to study the sensitivity of the results 
to various assumptions. 

Recall that the goal of this part of the work is to identify 
thruster characteristics and tljruster technology options that have sig- 
nificant impact on transportation system overall cost. Although the 
model incorporates many anticipated technology limits (e.g., thermal, 
ion optics) , the sensitivity studies also include parameter variations 
that might require violation of presently expected limits. Thus, inter- 
pretation of the final results in terms of new directions in technology 
will require certain qualitative judgements to temper the "black or 
white" computer answers. 

Initially, the cost model analysis was directed only toward large 
space systems. Later, the Shuttle baseline was included to provide a 
link between the present potential applications of the 30-cm thruster to 
relatively small satellites and to far future systems of the solar power 
satellite type. Many of the parameters studied in the large system base' 
line sensitivity analysis proved to be relatively unimportant. Thus, 
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Che number of runs needed for the Shuttle baseline analysis was signif- 
icantly reduced. For convenience in presenting the results, the large- 
system baseline is discussed first; the Shuttle baseline follows, draw- 
ing on similarities with the large system baseline to reduce duplication. 

To simplify the designation of the baselines and to distinguish 
between various orbit raising and on-orbit missions, the following 
abbreviations are used: 

OR = orbit-raising (transfer) mission 

ORS = self-powered orbit- raising mission 

ORP = payload powered orbit 

ORO = one-way orbit-raising mission 

00 = on-orbit mission 

LS = large space system era baseline 

MS = matured Shuttle era baseline. 

Within the OR/LS category, there are three modes (see Section 4.A.1): 
self-powered, payload powered, and one-way. For these, an additional 
descriptor is used, as indicated above. The self-powered mode is used as 
the baseline, and most sensitivities were investigated using this mode. 

Large Space Systems Era Baseline 

A summary of the large system baseline parameter values is 
shown in Table A-3 for use with the mass and cost equations listed in 
Tables 4-1 and 4-2, respectively. All of these parameters were pres- 
ented in previous sections. The selection of these particular values 
was accomplished through discussions with NASA LeRC personnel and 
reviews of preliminary results. Argon propellant was assumed because 
of the large quantities required. 

The results obtained with the computer cost model are catego- 
rized as indicated in Table 4-4. The general category presents the 
baseline results and several breakdowns of mass and cost to show the 
contributions of the various model elements. The other categories were 
selected for the sensitivity studies. In each case, the parameter of 
interest is varied about the baseline value with all other parameters 
fixed. The results are presented below by category. ^ 


126 



J able 4-3. OR/L S Ba seline Cost Model Paramete rs 
Argon Propellant 

= 100 kW 1 

Module power (beam) 


10^ hr 


6 km/ sec 
150 days 
10^ kg 


Curve 


2000 days 
$10/kg 

$0.5/W 

$4.2 X 10^ log (0) 

$4. 3 X 10^ (P 1^ 
PP 


Thruster life 

Velocity change 
Op trip flight time 
Initial mass in LEO 
Radiation flux 
Radiation model 
Life of ORV 

ORV Integration and testing 
Cost of ORV associated power 
First unit thruster cost 
First unit PPU cost 


"ot 

= 0.2 

PPU redundancy 

R 

op 

= 0.2 

Thruster redundancy 

C 

rpt 

= $100/kg 

Cost per kg of tankage 

pt 

= 0.10 

Tankage factor 

C’ 

orp 

“ $. 40/kg 

Propellant cost per kg 

C 

opt 

= $10^/day 

Time-dependent cost 

C 

rsm 

= $10^ /m \ 

\ rsm ) 

Cost of ORV structure and 

c . 

ord 

10 ^ orv orp/ 

mechanisms 

DDT&E cost amortized over 
10 ORVs 

C 

rss 

= lx 10^ /m \ 

\ rss f 

Cost of ORV subsystems 




Table 4-3. OR/LS Baseline Cost Model Parameters (Continued) 


^opf 

= $10^ 

Cost not associated with 
electric ORV 

“ps 

= 3 kg/kW 

Specific mass for power 


= $15/kg 

Launch cost per kg to LEO 

r 

= mjyx 

Interest rate 
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b. 

OR/LS General Results 



The typical data format and results for the baseline condi- 
tions are presented in Figure 4-12. The 100-kW curve represents the 
baseline case, while the other module power curves show the general 
cost trends with module power. Lines of constant thruster diameter 
imposed on the power curves indicate the general relationship of size 
to cost and 1^^. Unless otherwise marked, all data was developed using 
the baseline parameters. 

The shape of the curves in Figure 4-12 is the result of 
changes in many parameters. Figures 4-13 through 4-16 are provided to 
aid in the Interpretation. Breakdowns of the cost model and mass model 
elements as a function of are shown in Figures 4-13 and 4-14, respec- 
tively. In addition, plots of total cost and of cost per kilogram as a 
function of payload mass are presented in Figure 4-15. An additional 
aid in Interpreting Figure 4-12 is the number of thruster modules (N ) 
required, as shown in Figure 4-16. 

At low Igp» the steep slope of all the curves in Figure 4-12 
is primarily due to the rapidly increasing payload. Since, as shown in 
Figure 4-15, total cost is relatively constant, the reduction in payload, 
due to Increases in other masses as decreases dominate the low I 
region. (Payload is derived by subtracting the ORV mass and propellant 
from the initial mass in LEO.) The mass elements that Increase wIlIi 
decreasing are, as shown in Figure 4-14, thrusters, propellant (up 
and down), and structure and mechanisms. Structure increases because 
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Table 4-4. Data Run and Sensitivity Study Categories, OR/LS Baseline 


1- General 

Cost/kg versus Igp; 

Cost breakdown 
Mass breakdown 
Cost/kg versus payload mass 
Number of modules 

2. Mission design sensitivites 
AV 

Trip time (up-trip) 

Initial mass 
Launch cost 

3. Thruster sensitivities 
Life 

Efficiency 

Mass 

First unit cost 
"R" ratio 
Redundancy 

4. Propellant sensitivities 

Type 

Cost 


5. PPU sensitivities 
Mass 

First unit cost 
Redundancy 

6. Power source sensitivities 

Mass 

Cost 

Degradation 

7. ORV system design sensitivities 

Structure and mechanisms 
Subsystems 

Integration and testing 
DFT&E 

Flight operations 

8. Payload powered options (ORP) 
Power fraction 

9. One-way trip options (ORO) 

AV 

Trip time 
Initial time 
Launch cost 
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the model assumes a proportionality to total ORV mass which increases 

rapidly as I decreases. Another significant point is Illustrated in 

Figure 4-16. In the low I area, large thruster diameters are 

sp 

required, even for the 20-kW and 40-kW moduljs. Since a thruster diam- 
eter of less than 100 cm would be desirable, and since the number of 
modules is high (500 to 1000) for low module power, a 100-kW module for 
the baseline is reasonable. 

Moving to the high end of Figure 4-12, total cost (Fig- 
ure 4-13) is increasing while payload mass (Figure 4-14) is decreasing. 

The basic cause of both of these changes is the Increased power associ- 
ated with higher I ♦ Even with the low cost of power assumed for the 

baseline, raw power and PPU costs are important at high I • 

sp 

The most cost-effective range of I for the ORS/LS baseline 

sp 

is from about 4,000 to 10,000 sec (Figure 4-15), The transportation 

cost per kg of payload decreases rapidly as I increases to 6000 sec, 

sp 

As explained above, this occurs because payload mass is Increasing rapidly. 

At higher values of I , the specific cost increases again. This occurs 

sp 

because the payload mass is nearly constant while the costs, mainly for 
power and PPU, are increasing. Similar comments apply to the total cost 
as a function of payload mass. At the low end of this range, relatively 
large diameter thrusters or low-power modules would be needed. However, 
substantial cost penalties are incurred at low module power. The region 
around 8,000 to 9,000 sec and a module power of 100 kW provides a 
good compromise between low cost and thruster diameter. Lower module 
power would allow smaller diameters, but at a higher cost. Higher 
mocnle power would reduce costs, but would require larger diameters; 
modules of even 100 kW will present a challenge for test facilities. 

These preliminary conclusions will be discussed further in light of 
the sensitivity results presented below, 

c- Mission Design Sensitivities 

The sensitivities reported in this section are for variations 
around the baseline self-powered rojnd-trip mission. Other mission 
options are presented in Section 4,A.8.i. The sensitivities to mission 
AV, trip time initial mass in LEO, and launch cost can be summarized 
as follows: 
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As an example of the sensitivity of cost to lives shorter than 
10^, consider the baseline region discussed in Section A.A.S.b. 
At 9000 sec, the lOd-kW module has a diameter of about 55 cm and 
beam conditions of 50 A at 2000 V (see Figure 4-4) . With the 
beam current dependent model, a life of about 815 hr is pre- 
dicted. Since the cost difference is about $8/kg, the ratio 
of cost change to life change is approximately 9 x 10“^ 

$/kg/hr of life. Thus, for a typical payload of 8 x 10^ kg, 
the absolute ratio would be 'v $700/hr. Ir other words, in 
developing a new thruster for the ORL/LS aj Mcation, if the 
cost to achieve thruster lives longer than 'v 1000 hr is less 
than $700/hr, the increase would be cost effective. Since 
this ratio is nonlinear and is probably substantially less 
near 10^ hr, a life goal of a few thousand hours would prob- 
ably be reasonable. The main point to be made is that 
thruster lives of a few thousand hours are quite acceptable 
for the assumed application. 

• Thruster Efficiency . Cost sensitivity to thruster (or system) 
efficiency is shown in Figure 4-22. The baseline efficiency, 
as discussed in Section 4. A. 3, is only moderately optimistic. 

It assumes that electrical losses, utilization efficiency, 
and correction factors can be maintained at the same values 
as for the present 30- cm thruster with a 2 A beam. At higher 
beam current densities, measured utilization increases while 
the correction factors decrease, tending to produce a constant 
total efficiency. The probability is low for achieving effi- 
ciencies higher than the model predicts. The probability is 
moderate for slight reductions in the model predictions. But 
the probability is low for efficiencies significantly lower 
than predicted. 

Thus, the results indicated in Figure 4-22 show that system 
cost should not be too sensitive to reasonable changes in 
efficiency. More specifically, efforts to Increase efficiency 
are not warranted, and minor reductions from the predicted 
values will not substantially affect cost. 

• Thruster Mass . Cost sensitivity to thruster mass is shown in 
Figure 4-23. Except at low Igp, thruster mass variations 
from 0 to 200 percent of the baseline value have little impact 
on total cost. Such an effect is explained by the small frac- 
tion thrusters contribute to ORV total mass. In developing a 
new thruster, mass could be traced for other factors such as 
larger fabrication tolerances, lower cost, or increased life. 

• Thruster Cost . The sensitivity to thruster first unit cost is 
shown in Figure 4-24. Small variations around the baseline 
have little impact on transportation cost, although signifi- 
cant increases in cost would begin to be important. 
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Dean, voltagl'to'toui vIlLgerplays'TJoU ‘n deflT^" 

with the R ratio to the 3/2 ! thruster varies inversely 

even „lth ratloa aa Jew aa 0 T^'t f n T 
costs are not strongly affected^ \ht *'^^"®P°rtatlon 

diameter, since the pLer (thrust^ Ilf effect is in thruster 

invanaaly wUh the /"uo 3/2"p„“a“. 

tulll onia?? 

given reliability requirement with 

(particularly since 0 ? fr L ^ ^ ‘^ost 

arry since 0.2 provides good reliability). 

Appellant Sensitivities 

Two propellant parameters were Investigated: propellant type 
.. propellant eoat. Propellant tankage mass and tankage coat aenalLv- 

V mass (see Figure 4-U) that even large variations would not aftect 
transportation cost. 

xenon and'"’ ~ <«B0n. krypton, 

xenon, and mercury) were considered, as shown in Figure 4-27 Althouah 

the same cost was used for all propellants In Figure 4-27. the lack of 

sensitivity to type la valid. The W.4/kg cost la repreaentatlve of 

He^°X propellants. With higher costs for 

. Xe. and Kr. the cost curves will compact even more. The primary 
lerence etween the various propellants la thruster diameter at a 

100 kw , r' “>‘>“1^ p— of 

0 kw. the for a go-cm thruster would range from about 16.000 sec 
argon to about 7.000 sec with Hg as shown In Figure 4-4 Prooell . 

Ill™ !n° r'^ transportation “ 

DOS (e.g.. environment, availability, thruster development, or 
packaging) . 
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Figure 4-27. ORS/LS effect of propellant type 
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(2) Propellant Cost — The transportation cost sensitivity to 

propellant cost per kilogram is shown in Figure 4-“28 for argon. For 

I._’s in the range of 3,000 to 9,000 sec, propellant cost can affect 
sp 

transportation cost. However, for I 's greater than 9,000 sec, the 

sp 

dependence on propellant cost is rather small and should not influence 
ORV or propulsion system design. 

f . Power Processing Unit Sensitivities 

The variables considered in the PPU study Include specific 
mass, cost, and redundancy. PPU life was assumed to be long compared 
with ORV life, and wearout is negligible. In addition, PPU efficiency 
was omitted since the sensitivity result would be similar to that shown 
in Figure 4-22. 

(1) PPU Specific Mass — The results presented in Figure 4-29 

Indicate that mass decreases from the baseline would not significantly 

reduce transportation costs. However, large Increases in mass from the 

baseline would have a major impact. In the region of 9,000 sec, the 

transportation cost sensitivity is about $3/kg/unit change in a (base- 

PP c 

line value at 100 kW is 1.4 kg/kW) . Thus, with a payload of 8 x 10 kg, 

the sensitivity is about $2.4 x 10^ per unit change in a . Such sen- 

PP 

sitivity would probably provide Incentive to maintain low specific mass 
or to reduce specific mass if initial developments are significantly 
higher than the baseline. 

(2) PPU Cost — The Importance of PPU first unit cost is 
shown in Figure. 4-30. At 9,000 oec the sensitivity of transportation 
cost to PPU cost is about $0. 07/kg/percent change in PPU unit cost. 

In absolute terms with 8 x 10 kg of payload the sensitivity is about 

4 

$5.6 X 10 per 1 percent change in PPU cost. Since the PPU cost model 
is probably accurate only within a factor of two (i.e. , from 1/2 to 2 
times the model estimate) , a 100 percent increase in PPU cost would be 
quite significant. 
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(3) PPU Redundancy — The effect of PPU redundancy is more 
dramatic than for the thruster, as shown in Figure 4-31, since PPU cost 
represents a larger fraction. However, with reasonable redundancy fac- 
tors (i.e., less than 1.0), transportation cost will not be affected 
significantly. 

g. Power Source Sensitivity 

Power source mass and cost are among the most sensitive param- 
eters in the present cost model as well as in essentially all SEP appli- 
cations. These factors and degradation are considered here. However, 
under category 6 in Table 4-4, other power sensitivities are presented. 

(1) Power Source Specific Mass — As indicated in Eq. 4-24, 
solar power source mass is the product of specific mass (kg/kW) and 
total installed power (for propulsion). Transportation cost sensitivity 
to specific mass is shown in Figure 4-32. Current technology is about 
15 kg/kW, but projected advances from the literature were assumed in 
selecting the 3 kW/kg baseline. Using the 9,000 sec point again as a 
reference point, the transportation cost sensitivity is about 4 to 5 $/kg 
per unit increase is specific mass. With a payload of 8 x 10^ kg, the 

g 

sensitivity is $30 to 40 x 10 per kg/kW. Clearly, there is a strong 
incentive to minimize power source specific mass. 

(2) Power Source Specific Cost — Sensitivity to power source 
specific cost ($/kW) is shown in Figure 4-33. Although the baseline of 

3 

$0.5/W is about a factor of 10 less than current costs, studies of large 
power systems project such substantial reductions in the next few decades. 
The sensitivity at 9,000 sec is about 5.5 $/kg per unit change in power 
cost, or $44 X 10^ per $/W of power. This extreme sensitivity places a 
large premium on low-cost power for this application. 
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Figure 4-31. ORS/LS effect of power processor redundancy. 
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Figure 4-32, ORS/LS effect of power source specific mass 






(3) Power Source Degradation — It is assumed that the power 
source utilizes some form of photovoltaic conversion. In raising the 
orbit through the Van Allen belt, degradation may occur depending on the 
array design and on projected improvements in cell technology. The 
impact of degradation for two equivalent flux levels is shown in Fig- 
ures 4-3A and 4-35. As defined in Section 4. A. 5. curve A is representa- 
tive of current technology, while curves B and C are assumed parametric 
improvements. At the 9.000 sec point, the sensitivities are ;,bout 0.3 
and 0.5 $/kg per 1 percent change in degradation at 10^^ and 10 equiv- 
alent flux, respectively. These convert to absolute sensitivities of 
$2.4 X 10^ and $4 x 10^, respectively, per 1 percent change in degrada- 
tion. Although improvements in degradation are not quite as beneficial 
as cost or mass improvements, costs in the range of $2 to 8 x 10 would 
be incurred if degradation followed curve A or B rather than the assumed 

baseline. 


h. ORV System Design Sensitivities 

This category includes some of the system factors that woui^ 
go into developing and building an ORV, as indicated In Table 4-4. Since 
most of the values selected for these factors are rather arbitrary, the 
main Importance of this particular Bet of results Is in assessing which 
factors might have the most influence. With such sensitivity Informa- 
tion, a more logical or balanced direction of thruster technology can be 

chosen. 


(1) Structure and Mechanisms Sensitivity — As Figures 4-36 
and 4-37 show, within reasonable bounds, structure and mechanisms param- 
eters (C and a ) have little effect on transportation cost. Combin- 
T stn stu 

ing high mass fraction and high specific cost would have a greater 
effect but would not greatly influence the curve shape or thruster 

selection. 


(2) Sbsys terns Sensitivity — Because of the large exponent 
used in the subsystems cost model (see Eq. 4-37), subsystem cost can 
strongly affect transportation cost, as shown in Figure 4-38. Again, 
the shape of the curves is not greatly affected over a rather large 
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range of specific costs (C|.gg), and thruster selection will not be 
Influenced. However, this sensitivity shows that ORV electronics com- 
plexity should be minimized. 

(3) ORV Integration Sensitivity — In the range of C’ from 
0 to $100/kg, transportation cost is not strongly affected, as^Fig- 

ure 4-39 shows. As a figure of merit for assembly- line- type operations, 
automobiles are probably assembled for about $l/kg. The baseline 
selected was an order of magnitude greater, but even costs up to $lOO/kg 
would have only a minor effect on total cost. However, a more significant 
impact would be introduced if integration and testing exceeded $100/kg. 

(4) DDT&E Sensitivity — Design, development, test, and evalua- 

tion are always a large program cost factor. However, as shown in Fig- 
ure 4-40, unless the coefficient is raised, and the number of 

vehicles is reduced, each by an order of magnitude, DDT&E is not a large 
influence on overall transportation cost. 

(5) In-Flight Operations Sensitivity - Operations cost param- 
eter results are shown in Figures 4-41 and 4-42. As Figure 4-4i shows, 
fixed operations costs above $10^ would be needed to greatly change the 
results. For time-dependent costs, the daily rate would have to exceed 
$10^ to have an impact. Since a man-day might conservatively cost $400, 

$10 would imply 25 people/day or about 8 people on 3 shifts. It seems* 
unlikely that a significantly greater number would be necessary on a 
continuous basis. 

i* Payload Powered Option Sensitivities 

The results presented thus far have assumed the self-power 
(ORS) option discussed in Section 4.A.I. The option considered here 
involves using power from the payload (e.g., a power satellite segment) 
for the up-trip, but using ORV power for the down-trip (ORP option). In 
this option, a fraction of the up-trlp power (calculated as before) is 
assumed for the down- trip. The ORV was charged (mass and cost) only for 
this fraction. Then, the down-trip time is dependent on the available 
power. Three cases for this option follow. 
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Figure 4-39. ORS/LS effect of specific cost of integration 
and test. 
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Figure 4-40. 


ORS/LS effect of specific cost of design, 
test, and evaluation (DDT&E) . 


development. 
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(1) Baseline ORP Power 

° '■> 'M Percent (ORS ba.eUne) of J P^^p^ fPo„ 
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redncea the ORV power source uJZTlL IVi 

trip time Increases, but onlv about 50 “““ “““p- 

trlp time becomes eatremcly long a d ''' <i°wn- 

5»PPtta„t. in the 1 ran e o L " '-tors become 
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- -'3 case, the powr^so::::!:::*.':''--^ -- -urce _ 
“Tth 3 kg/R„ 1„ rhe baseline The T 
V-^Ht Shift from the resuits"' pZtVr T^ 
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^ OpMon Sena-f f- ivit^e^ 

The third mission option fORn't « 

higher orbit, alWs the ORV to remain withTi"'’' “ 

the ORV cost Is then charged to a slnel ' PPylPP** »" orbit, since 

13 higher than for multiple trips. One-lay t'ilp"' 

cchted for 4V, trip time. Initial mass Z i ='"=“1^11100 are pre- 

the benefit of includine th^ ' case, 

IPPlng the power source In the payload Is assessed. 
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Figure 4-44. ORP/LS effect of fraction of ORS 

baseline power for down trip power 
= 10 kg/kW) . 
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(1) ORO AV Sensitivity - Cost sensitivities ere presented in 
Figures 4-46 and 4-47 for Av for nnv id ^ ^ sented in 

nower r ' tor iV for ORV-supplled power and payload-supplied 

from b ""or ™ ^ between the two ORO modes varies 

rom about 20 peroent at d R„/see to about 50 percent at 10 R„/sec. with 

self-powered costs being greater. The difference between the ORS 

ofa™ l”r” “““ to about a factor 

Of 3,5 in transportation cosf tVihd 

e substantial cost benefit 

time for th!^ ^ sensitivities to trip 

. for two power source cases are shown in Figures 4-48 and 4-49 
For trips of 100 days or more, the differences are of the » , 

s- 7 5 tl Th ’ looreases significantly <Co 

7.5 times the ORS cost of Figure 4-18). This different nature for 

»• 

«..u. j.”. .r 

Figures 4-50 and 4-51, The f fov i- ^ ^ ° ^ 

. u difference between the ORS option (Figure 4-1 

and the ORO ORV-powered option is about 3.5 to 4 0 How. 7 

at lo' b : “IF about 10 percent cheaper 

g, but is about 40 percent less at 10^ kg. 

cost 1 5“=ltlvlty - Cost sensitivity to launch 

cost is shown in Figures 4-52 anrf A •i'j c 

ORV-powered and payload- 

powered modes, respectivielv a j Payj-oad 

.U etlvlely. A comparison of Figures 4-20 and 4-52 

hows that the one-way trip becomes about 3 times more expensive for 

e baseline value. If power is charged to the payload (Figure 4-53) 

he transportation cost will be about 10 percent less at SlOO/kg and ' 

30 percent less at 81/hg than the self-powered mode (FiguL 4-52) 
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Figure 4-47. 
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Figure 4-48. ORO/LS effect of trip time, self- 
powered . 
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-51. ORO/LS effect of initial mass, payload 
powered = 0) . 
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Figure 4-53. ORO/LS effect of specific launch cost, 
payload powered (C = M = 0) . 
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Ralalnt/Shuttle Er» Baseline Mission 

The m era baseline is intended to provide a link between present 
30-om projeoted applications and bS applications, which are probably 
several decades away. The differences between the IS and «S bahelines 
are c aracterlzed by the seven paraseters listed in Table 4-5 with all 
t er values as defined in Table 4-3. The new values were selected on 

oro1 Td r“ 1-“' Advanced than that 

t , ='■““1- 

(this limits the initial n^ss to about 25.000 kg and Increases the 
specific launch cost). ( 3 ) fur snaller systems, lower power modules 
should be considered, and (4) the ORV fleet size would not allow 
production-line assembly, which would raise integration costs, 
a* OR/MS General Restl^^g 

The general characteristics of the Ok/MS baseline results are 
Shown in Figures 4-54 through 4-58. The effect of module power t/ 
s own in Figure 4-54 for the range from 5 to 80 kW. The most obvl'ous 
ange from the 0R/b8 results shown in Figure 4-12 is the order of mag- 
nitude shift in transportation cost to about 81.000/kg. other diffeL 
nces rnclude reduced sensitivity to module power level and increased 

TusTJT""' '-on- 

to mini 

mum cost than with argon. At about 5,000 sec, a 80-kW, 60-cm 
^hruster would be close to optimum. Although the sensitivity to module 
ze appears to be small, the compressed logarithmic scale is deceiving 
n he region of minimum cost, the sensitivity is about 80 . 44 /k 8 per cm' 
e t roster diameter. For a payload of 18.000 kg (Figure 4-56)! the 
a solute cost is about 88 z lo^/cm. Thus, the difference between a 

82 T" ^ t'-PO'tation cost difference o, 

. he coot difference is about 85.4 x lo« for a 30-cm to 100-cm 
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Table 4“*5. OR/MS Baseline Cost Model Parameters^ 


^mod 

S8 

20 kW 

Module power 

M 

o 

- 

25 X 10^ kg 

Initial mass in LEO 

C . 
ori 

= 

$40/kg 

ORV integration and testing 

c’ 

rw 

- 

$100/W 

Cost of ORV associated power 

a 

ps 

a 

6 kg/kW 

Power source specific mass 


a 

$300/kg 

Launch cost to LEO 

Propellant 

a 

Xenon 



All other model parameters are the same as those in Table 4-3. 
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Figure 4-56. OR/MS mass breakdown. 
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transportation cost per kg of payload, s/ka x iir 
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Figure A-57. OR/MS cost/kg versus payload. 
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TOTAL transportation C0ST,$ x 10" 




difference. Since thla coet le per roundtrlp, the Investment In 
developing a larger thruster would be well rewarded. 

Cost and mass breakdowns are presented In Figures 4-55 and 4-56 
for the OE/MS baseline. Iransporutlon cost Is now dominated by launch 
coat and power coat with PPU cost third. Total cost is of the same 
magnitude (62 x 10 ) as for larger payloads In the OR/LS bssellne. The 
mass breakdown is similar to the 0R/L8 case with propellant and thruster 
MSS dominating at low 1^^, and power and PPU masses dominating at high 
I . Cost as a function of payload mass, shown In Figure. 4-57, Indicates 
the most cost-effective region to be around 5000 sec. As a reference for 
thruster requirements. Figure 4-58 shows the number of thrusters required 


as a function of and module power. Also shown are Hues of cons'«nt 
minlmuja thruster diameter. 

Power Sou rce Senaitlvltles 

Since power is a major contributor to cost and mass, varia- 
tlOM in these areas are also significant, as shown In Figures 4-59 and 

at 5M0 power costs, 

at 5000 sec. Is about 63.5/kg per 6/W. With an 18.000-kg payload, the 

nsltlvlty Is 66.3 , 10 per 6/W. Thus. If the power cost uncertain- 

ties were on the order of SSO/u hVio f-v-ss... ^ j 

9 u/W, the transportation coat would be 

uncertain by about S3 x 10 ° a 4.11: 

y Dout X 10 (l.e., ±15 percent of the total cost). 

Specific mass sensitivity, at 5000 sec, is about $3/kg per 

g W of specific mass. Thus, an uncertainty of 2 kg/kW with an 

. 00 kg payload results in a transportation cost uncertainty of 

a out $1 X 10 . In this case, specific mass is probably not as impor- 
tant as power cost. 


Initial Mass Sepsitivities 

In the range of lo'* to 10= kg. as might be provided by one or 
a few Shuttles, the transportation cost Is only moderately sensitive to 
initial mass as shown In Figure 4-61. At 5000 sec. the cost sensitivity 
is about 60.002/kg per kg of Initial mass. With an 18,000 kg payload, 
this translates Into s sensitivity of about 630/kg of Initial mass. Thus 
a Shuttle could deliver only 2.5 x 10 « kg per trip, a cost-effective 
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Figure 4-61. 031/MS effect of Initial mass. 


approach would be to deliver several Shuttle loads to LEO before transport 
Using the above values, transporting four 2.5 x 10^ kg loads in one ORV 
trip would effectively reduce the total trip coat by $30/kg or by $3 x 10 

d. Launch Cost Sensitivities 

The effect of launch cost on total transportation cost is shown 
in Figure 4-62 to be quite significant. The sensitivity around the base- 
line at 5000 sec is about $1.6/kg per $/kg of launch cost. With an 
uncertainty of 100 $/kg in launch coats, transportation costs for an 
18,000 kg payload would be uncertain by about $3 x 10 . 

e. Integration and Testing Sensitivity 

With the higher baseline values for power and launch cost and 
the low initial mass, variation of integration costa become relatively 
insignificant, as Figure 4-63 shows. Integration costs significantly 
above $100/kg would be required to affect overall transportation cost. 

B. ON-ORBIT MISSION 

The basic philosophy used in the orbit-raising analysis is applied 
to the on-orbit mission cost modeling. With only minor exceptions, the 
mathematical formulation follows the orbit-raising equations. However, 
the results of the on-orbit modeling are expressed in terms of the total 
cost of providing propulsion per unit of satellite mass per year. This 
figure-of-merit, which includes the cost of transportation from Earth to 
GEO, produces a convenient normalization and a final number of tractable 
size (typically in the range of 5 to 500 $/kg/yr) . 

Applicable portions of the orbit-raising analysis are not repeated 
in this section. Specifically, the propulsion system performance, 
redundancy, and degradation sections are not duplicated. In addition, 
in describing the on-orbit computer model logic, the mathematical por- 
tions are referred to the corresponding relations in Section 4. A. 

Terminology and general definitions used in the on-orbit analysis 
are presented in Figure 4—64. This diagram is identical with that 
shown in Figure 4-1 except for terminology. For the on-orbit mission, 
the satellite stationkeeping and attitude control functions are pro- 
vided by the EP system (i.e. , the OOP). A power source mass allocation 
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Figure 4-64. 00? block diagram. 
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is applied to the OOP, although in most cases the power source physically 
would be part of the basic satellite. Depending on the satellite design, 
the power source might be required only part of the time and could be 
used by the satellite during non-propulsion periods. 

!• Mission Description 

The on-orbit mission is illustrated conceptually in Figure 4—65. 

The OOP is shown as a unit, but in most applications thrusters would be 
placed at widely separated locations to provide moment arms for attitude 
control. In general, the number and location of thruster^ ,1ns tailed 
would depend on satellite si^e and control requirements,. Since such 
satellite design variables cannot be easily modeled, an alternate 
approach using estimated iV requirements is included in the on-orbit 
cost model, as discussed in Section 4.B.2. 

As in the orbit raising analysis, specific details of the orbit 
(such as altitude, inclination, or eccentricity) are not directly 
included in the model. It. is assumed that the AV selected will ade- 
quately characterize mission requirements. The sensitivity of results 
to AV is discussed in Section 4.B.5. Choosing a conservative AV should 
produce a conservative final result , since the total propulsion system 
and propellant scale with AV. For large satellites, the forces and 
torques illustrated in Figure 4-65 are significant contributors to AV. 

An on-orbit AV model is discussed in Appendix E. 

2* On-Orbit Mission General Model 

The logic diagram for the on-orbit mission is shown in Figure 4-66. 
Comparing Figure 4-66 with Figure 4-3 indicates that the only signifi- 
cant difference in logic is the down trip option included in the orbit- 
raising analysis. The mathematical formulation, except for details of 
the cost and mass models, is the same as that presented in Section 4. A. 2. 

The starting point is the definition of the mission in terms of 
mission life, AV, and satellite mass. Satellite mass and size (area) 
must be considered in estimating AV (see Appendix E) . The basic vari- 
able, I^p, is used in a looping process to develop cost results as a 
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On-orblt model logic diagram 
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mass of thrusters (Eq. 4-19) 
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mass of ppus (Eq. 4-22) 
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op 
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propellant mass 
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at 

satellite Initial mass 
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« 

propellant tankage mass (Eq. 4-23) 

M 

ow 

«c 

power source mass (Eq. 4-24) 

M 
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as 

subsystems mass (Eq. 4-25) 

M 

osm 

m 

structure and mechanisms mass (Eq. 


4o Cost Models 
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table 4*-6* Mass Models Used in On-Orbit 
Cost Model i 


+ M^ + M^ + M + M + M 
Qt opp op opt OSS osm ow 


- 2 (o.078 N 


**opp 

Where, 


O.OOl (1 •»: R ) 

pp ps p 


a 

PP 

- 

jQ [-0.503 log + 1.153] 

M 

op 

■ 

[1 - exp (- AV/g^l 3 p)] 

M . 
opt 

^ m 

*^rpt ^op 

M 
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- 

100 [d + 10-*Pj„j) + (N„c/N„p) - x] 

M 

osm 

as 

“am [“ot + “opp + % + “opt “oas] 

M 

ow 

- 

“rw W 






Table 4-7. Cost Formulae Used In On-Orbit Cost Model 
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where 


ot 


opp 


opt 
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OOP ^ot ’*■ ^opp ■** %t: ®ow ^oss 

“osm “to + “ol + “ops * “odd* 


- cost of thrusters (Eq. 4-32 without last term) 
= cost of PPUs (Eq. 4 - 33 ) 

propellant and tankage (Eqs. 4-34 and 4-35) 
power source cost (Eq. 4-38) 


c 

OSS 

s; 

subsystems cost (Eq. 4-37) 

C 

osm 

« 

structure and mechanisms cost 

^to 

sa 

transportation cost 


- 

C' M 
to 08 

“io 

= 

transportation cost from Earth 

“ol 

= 

integration and testing cost 


= 

C'a (M - M ) 

ol oop op** 

C '4 

ol 

= 

Integration cost per kg of dry 

C 

ops 

B 

operations cost 


5: 

C t 

ops ^f 

^f 

9 

mission time 

^odd 

B 

DDT&E cost (Eq. 4-40) 


N 


ov 


number of on-orblt vehicles over which DDt&p 
are amortized wnicn ddx&E costs 


(4-43) 
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^-Orblt/Large Syste m Baseline Mission Reall^^a 

In a fashion similar to the orbit-ralslng mission, two baseline 
sets of parameters are used in developing on-orbit results; large space 
system era (00/LS) and matured Shuttle era (00/MS), Again, for each 
baseline, general results are presented to illustrate the characteristics 
Overall cost sensitivities for selected variables are then analysed to 

assess the Importance of various parameters, particularly thruster 
parameters. 

Orbit /Lar ge System Baseline Parameters 

Tha Urge system baseUne (00/LS) parameter see Is shovm In 

Table 4-8. The parameters that are different from the OE/LS baseline 

reqnlre some discussion and Inelnde AV. mission time, and transportation 
cost* 

An analysis for estimating AV Is presented In Appendix g 
Clearly, without a apeclflc satellite design, such an analysis Is ' 
extremely approximate. However, since the primary purpose of the ansly- 
s s was to determine the order of magnitude of AV to be considered, the 

absolute accuracy Is relatively unimportant. A baseline AV of 

1.65 km/s/yr was selected. 

A baseline mission time of 10 years was chosen, based on 

several factors. Although longer times may be needed to make certain 
satellites (e.g.. solar power satellite) cost effective. It la probably 
unreasonable to store the necessary propellant for such periods. Thus. 

If propellant Is replenished, propulsion system modules could be replaud. 
These additional concept variables would substantially complicate the 
present problem. Mission times of 3 to 10 years are currently being 
achieved with exUtlng satellites. However. sUce current ssteUltes 
are relatively simple compared to those Imagined U the Urge system base- 
line. a nominal 10 year mission Is a reasonable compromise. 

^ Transportation cost, as It Is used In the on-orblt cost analy- 

sis. Includes all costa from Earth to the satellite orbit. This cost Is 
Just the value computed U the orbit raising mlsston. A value of MO/kg 
was selected for the 00/lS baseline based on typical OS/LS baseline results 


4~8» 00/LS Bdscllns Cost Modol Parainetors 



Argon Propellant 

^od ■ 100 kW 

Thruster Life = 10^ hr 

AV =16.5 km/sec 

= 10 years 

Initial mass = 10^ kg 



15 2 

flux ® 10 dGctrons/cm ^ using curve 

c* 

ol 

= $10/kg 

OOV Integration and Testing 

C* 

ow 

= $500/kW 

Cost of OOP Associated Power 

CV 

ot 

= $4.2 X 10^ ■ log (D) 

First Unit Thruster Cost 

C 

op 

“ $4.3 X 10^ (kW)°'^^ 

First Unit PPU Cost 

\ 

= 0.2 

PPU Redundancy 

R 

P 

= 0.2 

Thruster Redundancy 

C ^ 
opt 

= 0.4 (1 + 0.1 X 100) 

Propellant at $0.4/kg 



Tankage at $100/kg with 

Tankage factor = 0.1 

C' 

ops 

= $10^/day 

Time Dependent Cost 

Y 

rs 

= 10 

Life of OOV 

C 

osm 

= $10^ (M 

osm' 

- ^ * “■ (M -M )“-0 

10 oop op 

Cost of OOV Structure and 
Mechanisms 

^odd 

DDT&E Cost Amortized over 

10 OOPs 


= $1 X 10^ (M )°*^ 

OSS 

Cost of COP Subsystems 

a 

ps 

= 3 kg/kW 

Specific Mass for Power 

C’ 

to 

= $50/kg ' 

Orbit Raising Coat per kg 




Ba seline. General Reanl^a 

The general results of the 00/LS baseline are shown in 
Figures 4-67 through 4-71. sp.ol,u cose «/hg/yr) ... g„„etlo„ of I 
or . r«.ge of module power 1. show, lu figure 4-67. At low I , these 
here . shepe simller to the orbit rsieiug results. Howev^?. .t 
gh l,p, the curves eoutinue to decre..e or Just flstten out. The 
eupwion for this hehevior is Ulu.tr.t«i in figure. 4-68 .„d 4-69, 
which .how cost .ud .... hreskdewns. At low I cost is high (due 
to the Urge tr.usport.tton cost) .ud psyioud ..L is low (due to 
Urge propellent, thruster, .mi structure messes). At high I . propel- 
l«.t .ml most other messes ere decreesing. Power source m.s.“i. 

ucreMiug hut is . ...11 fr.ctiou of the totel. Although tot.1 cost 
is increesiug et high its t.ti„e„oe „„ . 

the increasing payload mass* 

specific cost end totel cost es . function of the peyloed 
was. ere shown in Figure 4-70. for the bssellne conditions, minimum 
specific cost ,s weU .. minimum total cost, occurs in the region of 
12,000 sec The sensitivity of specific coat to module pouer at this 

tra’ “• “l‘» 

e 8 . 10 kg payloed, the cost/yr is .bout $10^ per kH of module 

TM^hw «a ”r “ 50-kW (30-cm) module system end 

00-kW (60-cm) module system would be about $1.5 r loV or $15 . lo‘ 
over the mission life. Tgu 

number of thrusters regulred as a function of I 1 . shown in figure 4-71. 
The dashed lines show minimum thruster diameter!^ 

using the baseline results as . stating point, several sensitivity 
studies were performed to assess the importance of various parameterT 

Table IT' *” «■= '-lls-'l-g sections, are shown in 


sp 


207 






X of Total Cost 
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CURVK 'RlEMtENT 

1 . cor 

Tlhruaccr* 

2. COR 


3, COPT 

Propolluit 4 tankage 

4, COSH 

Struct. & aMclia 

S. COSS 

Subayatiima 

6 . cro 

Tranaportatlon *■ 7 

7. COW 

Power eourca 

s, COPS 

Operationa 

9. COPO 

DDT&K 

10, coi 

Integration 

Total Cost 

- A 


Figure 4-68, 00/LS cost breakdown 
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CURVE ELEMENT 

1 . MOT Thrusters 

2. MOPP PPU 

3. HOP Propellant 

4. MOPT Tankage 

5. MOSM Structure S mechanisms 

6. MOSS Subsystems 

7. MOW Power 

8. MOSP Payload 










c* MsBlon Design Sensltivltipa 

Four genoral uIobIou factors affecting coat are discussed In 

Alth u'"'’ transportation cost. 

^ oug eac as a significant Impact on coat and could Influence 

mission design, none appears to substantially change the general shape 
Of the curves or thruster considerations. 


Table 4^9. OQ/LS Data Run and Sensitivity Study Categories 


- P 
sp* mod 


1 1. General 

Cost/kg/yr versus I 
Cost breakdown 
Mass breakdown 
Cost/kg/yr versus ps.yload 
Number of modules 

2. Mission design sensitivities 
AV 

Mission time 
Initial mass 
Transportation cost 

3. Thruster sensitivities 
Life 

Efficiency 
Mass 

First unit cost 
Redundancy 


4 . Propellant sensitivity 
Cost 

5 . PPU Sensitivities 
Mass 

First unit cost 
Redundancy 

6. Power source sensitivities 

Mass 

Cost 

Degradation 

7 OOP system design 
sensitivities 

Structure and mechanisms 
Subsystems 

Integration and testlns 
DDT&E 
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(1) AV Sensitivity Results for a range of AV requirements 
are shown in Figure 4-72. For a given mission time and initial mass* 
such a range implies a' range of satellite area-to-mass ratios or uncer- 
tainty in the AV estimate. The minimum cost region around 12,000 sec 
should not bo significantly shifted over this AV range. 

(2) Mission Time — The effect of mission time on on-orblt 

propulsion costs is shown in Figure 4'-73. Mission times significantly 
shorter than 10 yr Increase special costs dramatically and tend to reduce 
the optimum Times longer than 10 yr will decrease the specific 

costs almost linearly. 

(3) Initial Mass The initial satellite mass, Including the 
OOP and propellant, has a large impact on specific cost, as shown in 
Figure 4-74. Although this effect la accentuated by the constant AV 
used in producing these curves, the general sensitivity Is still sig- 
nificant. The reason for this trend is that, since a reasonable fraction 
of the OOP costs are Independent of mass, specific costs Increase ae 
initial mass decreases. 

(4) Transportation Cost — The sensitivity of specific cost to 
transportation specific coat is prcaented in Figure 4-75. Transportation 
cost begins to have a significant effect on specific cost above about 
$10/kg. As discussed in Section 4. A. 8, transportation will probably be 
in the $30 to $50/kg region for the large system baseline, and on-orbit 
costs will certainly be Influenced over this range, 

d. Thruster Sensitivities 

The thruster sensitivity factors considered for the OO/LS base- 
line are shown in Table 4-9. Since it was found In the OR/LS cases that 
the "R” ratio has only a minor impact on cost, this parameter was omitted 
for the on-orblt work. 

(1) Thruster Life — Two models of thruster life, constant and 
beam dependent, were evaluated as shown In Figure 4-76. Constant life 

4 

greater than about 10 hr has a relatively small Impact on cost. However, 
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Cost 



Flg„„ 4.72. oo/LS effect of 4V 
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shorter life or a beam-current-dependent wearout characteristic would 
substantially increase specific cost. In addition, the sensitivity to 
becomes more pronounced for the variable life model. The optimum 
occurs between 3,000 and 6,000 sec and requires a large>diameter thruster 
(larger than 100 cm). However, if thruster diameter is limited to less 
than 60 cm, the coat sensitivity to I^p is extremely small (for either 
life model) . 

(2) Thruster Efficiency — As shown in Figure 4-77, thruster 
efficiency variations over a reasonable range have little effect on 
specific cost. Efficiencies below O.S would have a larger impact since 
more modules would be needed. 

(3) Thruster Mass — Thruster mass has only a moderate Impact 
on cost, as Figure 4-78 shows. Only increases of several times would 
significantly change the cost level. 

(4) Thruster Cost — The specific cost results are only 
moderately sensitive to thruster first unit cost as indicated in 
Figure 4-79. Again, increases of a factor of two or more would be 
required to change the specific on-orbit costs significantly. 

(5) Thruster Redundancy — The baseline redundancy of 20 per- 
cent should provide a reasonable reliability. Even if a redundancy of 
40 percent were used, the specific cost would not be strongly affected 
as Figure 4-80 shows. 

e. Propellant Sensitivity 

Based on the orbit-raising results, only propellant specific 
cost was considered here. The results are shown in Figure 4-81 for a 
range of propellant costs. Except around 6,000 to 9,000 sec, propellant 
cost has almost no Influence on overall specific cost. Other propellants 
should produce similar results with only a shift in thruster diameters. 

f . PPU Sensitivities 

The sensitivities considered are indicated in Table 4-9 
and Include mass, cost, and redundancy. 
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Figure 4-77. 00/LS effect of thruster efficiency. 
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(1) PPU Mass — As Figure 4-82 shows, PPU unit mass 
varla£?,ons around the baseline have a small Influence on cost, 

(2) PPU Cost — PPU cost is somewhat more important as shown 
in Figure 4-83. However, modest variations around tne baseline would 
not dramatically change the overall cost. 

(3) PPU Redundancy — Results shown in Figure 4-84 Indicate 
that PPU redundancy is also not too important in overall cost. A factor 
of two over the baseline should be an upper bound on the required 
redundancy. 

8* Power Source Sensitivities 

Results of the three major factors studies (Table 4-9) are 
shown in Figures 4-85 through 4-87. Power source specific mass 
(Figure 4-85) and degradation (Figure 4-87) show only small overall cost 
sensitivity. Power source specific cost ($/kW) is a more Important 
factor, as Figure 4-86 shows. In addition to changing the overall cost 
level source, specific cost changes the curve shape and hence changes 

the Igp optimum. At even higher power costs, this effect would be more 
pronounced. 

h. OOP System Design Sensitivities 

Representative system parameters (shown in Table 4-9) were 
investigated. Although many others could have been considered, the four 
selected should provide an indication of the sensitivity of overall cost 
to system factors. 

(1) Structure and Mechanisms Sensitivity — As Figure 4-88 
shows, structure costs below $10^/kg have little effect; costs of 
$10Vkg will increase the overall specific cost level but should not 
significantly, influence thruster selection. 

(2) Subsystems Sensitivity - The subsystems cost parameter 
essentially reflects complexity: the greater the complexity, the 
higher the cost. Figure 4-89 shows that any significant Increase from 
the baseline ($10 /kg) will have a pronounced effect on overall cost. 
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Figure 4-87. OO/LS effect of power source degradation. 
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(3) Integration and Testing Sensitivity ^ Th« 
test parameter illustrated in Figure 4~90 is al i tegration and 

costs above ^O/kg shift the overall cost level IndTT 
thruster and I aalactlon. ^ 

sp 

«) DBTSE Senaitivlty - As ahown In Figure 4-91 .. 
from the bas^Hno ^ s a an Increase 

be n magnitude in DDT&E cost would 

be necaasary to have nmch etfact on overall oo«t n 

DDTSe cost vould not slanlflcantl i «i ’ Wsher 

dot significantly Influence thruster or I selection. 

- ‘^-"■^■‘“^Mdtnred Shuttle Er. Mlhslon °'’ 

The on-orbit/matured Shuttle (00/MS) basellnin 

rr“:“ -■ *•“ = 

o-.. «. .,r,r.:!T :;™ “ r“ 

in Table 4-5. ^ vaj^ent to the values shown 

00/MS Ge neral Reall^^p 

thro H a ='““■> ^isores 4-92 

tba ratge oTlntl:r;Tt:to'’T"“ ^ 

«.000 sec. Further«,r . as n“ rel-BTsW Y 

snail thruster dlaneters <Je : La ^30 LToy- “ 
b:;.. significantly higher Cahoot 20 tlnesl than'that ofTLlT 

A cost breakdown and total cof^^ aa ^ 

in Figure 4-93- o h function of I are shown 

8 93, a mass breakdown is presented in Figure 4-94 Sner^f4 

cost and total cost are she™ m Figures 4-94 and 4-95 as a Lnctl L 

net satellite mass. Although total cost Is alnlnnut at about 7 000 ” 

e spec Ic cost nlnlnun occurs at about 10.000 sec for a payload na“s 
of about 18,000 kg. ^ yj-oaa mass 

Power Sou rce SenBi^!■^v^^a.o 

ores 4-96 I^d ““ ““ ^^0' 

4 96 and 4-97. respectively. Po„er source cost has a relatively 
large influence on overall ^ 

Otars less than 30 cu are of 11 tL L TZ"’ 

interest, the sensitivity is 
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Figure 4-90. 00 /LS effect of integration and test 

specific cost. 
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Table 4-10. 00/MS Baseline Cost Model Parameters* 


^mod 

1 

M 

20 kW 

Module power 

M ■ 

i' ■ 03 

- 

25 X 10^ kg 

Satellite Initial mass 

C*. 

oi 

- 

$40/kg 

Integration and test cost 

C* 

ow 

<8 

$100/W 

Power source cost 

a 

ps 

S 

6 kg/kW 

Power source speoHlc mass 

C 

to 

St 

$ 1,000/kg 

Transportation coat 


m 

Xenon 


*A11 other 
Table 4-8. 

cost model parameters 

are the same as those in 

■ , 
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Figure 4-9/|. 00/MS mass 










somewhat reduced. An Inereaae la power coat f 
slsnmcant than a decreaae. The e«act“rol 

le shown In Figure 4-97 to be of M ^ ^ specific mass 

menge of a's. over a rather large 

-Initial M ass 

The effect of Initial ebso on overall coar t . 
ure 4-98. On-orbIt specific costa varp bp IblVa f 

s given 1 ^^) over the mass range considered Higher " 

lower cost (as before). Thus i-h ’ Sher mass results in 

ing one large satellite rather thaT " advantage in build- 

situation occurred. ones if such a 

^^^gagEorta^ Cost Sen.f ^^,r.s■^. 

oost is relatively sensTwvrtl^trl^I interest, specific 
nre 4-99. the $1.000/kg baseline is rTpTerenrrtr'‘f """ 

«sults. but Iwprovements 1 „ orbit raising . 

benefit the on-orblt mission. “*«"i«osntlp 

^3agratla^nd_l estlng Senslt..H-,. 

Integration costs are a small fraction nf 
relatively large changes have essential! T 
This result is shown in Figure 4-100. ^ ^ specific cost. 

C. COST MODELING CONCLUSIONS 
^^^Eustersi ^ecommendar^ov,. 

»£ costs snd thruTtsr Lrt^lol^." a”^ ““ ^ 

Shown in Table 4-ii _ n™iary of general conclusions is 

key parameters that ‘provide aZToTT 

and thruster development difficulties system costs 

— casts srs wlolnicsd sronnd MOO arcr^^r 
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Figure 4-98. 00/MS effect of initial mass 
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Figure 4-99 . 00/MS effect of specific cost of transportation 

from Earth to GEO. 




Table 4-11, Siunmary of Cost 


Sensitivities 


Modeling Recommendations and 


Parameter 


Recommendations 

^ap* 

Thruster diameter, cm 
Module power, kW 
Propellant 

Thruster sensitivities 
Module power 
Life 

Efficiency 

Mass 

Unit cost 
”R" ratio 
Redundancy 

Propellant sensitivity 
Type 
Cost 


Mission Option 


OR/LS 


8,000-12,000 

40-70 

50-200 

Argon 


High 
Moderate 
Moderate 
low 
. Low 
Low 
Low 


OR/MS 


4,000-5,000 

30-70 

20-80 

Xenon 


Moderate 

N/S 

N/S 

N/S 

N/S 

N/S 

N/S 


00/LS 


10,000-15,000 

40-70 

100-200 

Argon 


High 

High 

Low 

Low 

Low 

N/S 

Low 


00/MS 


4,000-7,000 

30-60 

20-80 

Xenon 


PPU sensitivity 
Mass 
Cost 

Redundancy 

Power source 
sensitivity 
Mass 
Cost 

Degradation 


Moderate 

High 

Low 


High 

High 

Moderate 


Moderate 

High 

N/S 


Low 

High 

Low 


Low 

Moderate 

Low 


Low 

Moderate 

N/S 


6 ] 
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Table 4-11. Summary of Cost Modeling Recommendations and Sensitivities (Continued) 


Parameter 

Mission Option 


OR/LS 

OR/MS 

00/LS 



00/MS 

Mission design 
sensitivity 

AV 

Mission time 

Initial mass 

Launch cost 

High 

Moderate 

High 

High 

N/S 

N/S 

Moderate 

High 

High 

High 

High 

High 

N/S 

N/S 

High 

High 

System design 
sensitivity 

Struc. & mechs. 
Subsystems 

Integration 

DDT&E 

Flight ops 

Low 

Moderate 

Moderate 

Moderate 

Low 

N/S 

N/S 

Low 

N/S 

N/S 

Low 

Miderate 

Moderate 

Low 

N/S 

N/S 

N/S 

Low 

N/S 

N/S 

Self powered 
sensitivity 

Power fraction 

Moderate 

N/S 

N/A 

N/A 

One-way trip 
sensitivities 





AV 

Free power 

Trip time 

Initial mass 

Launch cost 

>1 

High 

High 

Moderate 

High 

Moderate 

N/S 

N/S 

N/S 

N/S 

N/S 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A = Not applicable 

N/S = Not studied 
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power. However, 6,000 sec results In large thruster diameters (or high 
cost). To stay within reasonable bounds (less than 100 cm) and retain 
low system coats, a minimum of about 8.000 sec la needed. With high 
module power, costa are not too strongly affected up to about 12,000 sec. 
With these bounds on 1^^, the range of thruster diameters and power are 
established. A similar rationale was applied to the other three mission 
options. 

Having established the acceptable ranges of key parameters in 

Table 4-11, the results were reviewed with an eye to defining a single 

thruster design" that would satisfy all four mission options. Thruster 

diameter is then bounded between 40 and 60 cm. To satisfy both large 

system baseline mission options, an I in the range of 10,000 to 

sp 

12,000 sec is needed. This range with argon corresponds to a 
thruster diameter range of about 50 to 40 cm, respectively, and a module 
power of about 100 kW (see Figures 4-4 and 4-6). A 40-cm thruster 
operating on xenon at 5,000 sec (the overlapping point for the MS base- 
lines) would have a module power of about 30 kW; a 50-cm thruster would 
have a module power of about 45 kW. 

The beam voltages for argon at 10,000 and 12,000 sec would be 
approximately 2,400 and 3,500 V, respectively. For xenon at 5,000 sec, 
the beam voltage would be about 2100 V. If a 50-cm thruster were designed 
for a beam voltage of 2400 V, it could be operated on argon to produce 
up to 100 kW of beam power or on xenon up to about 50 kW at 5400 sec. 

Thus, a 50-cm thruster is a good compromise for the four mission options 
and also allows for a common beam voltage. 

A further consideration in selecting between a 40-cm and a 50-cm 
diameter thruster is the uncertainty of wearout life. Since module 
power is important to system cost effectiveness, if sufficient life 
could not be obtained with a 50-cm thruster, operation at lower beam 
current would be less significant than for a smaller thruster. At 50 cm, 
the options for maintaining high module power and some flexibility in 
wearout life are substantially greater than for a 40-cm design. 
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Sensitivity ConcliiB<ftnn 


R««lta «£ tha .wsltlvlty studies sr. su-srlMd qusUt.tlvely In 
Table 4-11. From a thruster standpoint, nodule power, life, and e*fl- 
clency are noat Important. The l.port<u.ce of high power was reflected 
to the previous recoanendatlona. U£. la Important, particularly for 

throuL°''''d'^ To a certain extent, life can be compensated for 

through retodancy. However, a sdnlnun life of leas than about 5.000 hr 
would penalise the syaten. Hlghar efficiency than that used In the 
"Odellng would not significantly improve the coat picture, lower effi- 
ciency would not be desirable. The fact that thruster mass has a anall 
influence Indicates mass ccdd be traded for life, efficiency, or cost, 
t cugh reliability la covered by redundancy assumptions, a trade of 
ruater mass for a simpler, more reliable thruster would also be prudent, 
Somewhat unexpectedly, the type and coat of propellant la relatively 
unimportant. Propellant selection certainly influences thruster oper- 
ating parameter selection, but the flexibility added by the availability 
of several propellants la quite significant. 

PPH and power source coat and mass were found to be extremely 
wportant In overall system cost effectiveness. Coats cj the magnitude 
use n the model should he sought at the possible expense of efficiency 
reliability, or mass. (PPH efficiency was effectively Included under ’ 
ruster efficiency sensitivity.) without light-weight, low-cost power 
sources, the overall system cost will be Increased significantly. 

From a system design point of view, this work verifies the fact ' 
that system complexity, as represented by subsystem, integration, and 
development costs. Is relatively important. Altbough the paraa«ter 

values are somewhat arbitrary, the point at which these parameters 
become Important Is defined. 

Of the mission design parameters considered. It la clear that 
transportation cost la the most significant since the other parameters 
are determined hy the mission selected. Since the orbit-ralalng mission 
coat la strongly dependent on launch cost, and on-orbit cost Is dependent 
on orbit-ralslug cost, low-cost transportation to IHO is a fundamental 
requirement for a cost-effective near-garth space tranaportatlon ayatem. 


252 


SECTION 5 

TECHNOLOGY ASSESSMENT 


The source ct the technology-relst.s 
described 1„ Sections 3 end 4 was a thr “ 

me areas covered Included 1 “ assessment 

»al agitations, wearout Ills, thr„steVIcaU„V”“’ 

operatlonal simplifications to reduce ' '■'“Pliant type, and 

~s. ibese areas are dlscusrinlls I'l"-- " ^ 

A. ION OPTICAL SYSTEM 
Perveanrp 

All other thruster'deslgn pTZtZ IrH a 

“ °P““ design. Having specified . ‘he selection of 

and beam voltage required to produce the ^ hd® nnrren.t, 

apeclflc impulse, the first question ‘hrust and 

the size of the thruster (“grid area) ^ answered Is that of 

flcatlons. me lo„ thrus.er being co"7°°°'^ ‘hese specl- 

‘Tseem that has numerous allgneTa “““ “ ‘»-leo«ode 

and focus the Ion beam. Beam curl T 
Tnantltles 1 „ that the mamlmum clU t '' 

by apace charge and Is related to the el amtracted Is limited 

<^«d.s law, mis limitation dictates I'dT^ 

tlon electrode assembly. Analysis of th ‘he accelera- 

P«cesses has been performed by several ““ -“deration 

and presented 1 „ well-referenced form h 7 « “"-"‘^‘ad 

analysis here will be llmltedl IlmLg ‘ha 

‘hruster (lo„ beam, diameter to beam voUage'lT" 

Itage and current requirements. 


First, the beam current Ig and voltage Vg are related by 

3/2 


^b - " \ 


(5-1) 


where 


V. 


V. + V 
b a 


(5-2) 


and Is the negative voltage applied to the accelerator electrode to 


prevent electron backs treaming from the neutralized ion beam into the 
interelectrode acceleration gap. The perveance P, constant for a 
gi\<in electrode design, la a function of aperture and interelectrode 
dimensions, and can be written 


P » C, N P 


1 a a 


(5-3) 


where is the number of apertures in the acceleration electrodes, and 


P^ is individual aperture perveance. Figure 5-1 illustrates the aper- 


ture configuration and gives dimensions of the 700-900 series EMT 
design. In terms of the symbols used in Figure 5-1, 


,-9 


4.4 X 10 K 

TT2 ^ ■ 

).336 +4 ' 


( 0 . 336 . 


(5-4) 


for singly charged mercury lons,®^ and 


H + t 
8 s 


The dimensions shown in Figure 5'1 yield 


P^ = 6.4 X 10“® 
a 


The factor C- accounts for the variation in current as a function of 


beam radius that la a consequence of a radially varying discharge plasma 
density distribution. Essentially, 






4/^ 27rrJ (r)dr 


J ttD 
max 


(5-5) 
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CM IN. 
d,- 0.191 0.076 
d^> 0.162 0.060 
S - 0.211 0.087 
- 0.066 0.022 
t, « 0.0^ 0.016 
t, « 0.061 0.020 


Figure 5-1. Aperture configuration for EMT ion optics. 


and, for typical 30-cjn operation, Cj^ >• 0.52. With the number of 

apertures set at 15,173 for the EM design, the perveance computed using 

Eq. 5-3 Is P = 5.1 x lO”^. However, the maximum perveance measured 

-5 

experimentally with the EM Ion optics design Is only 4 x 10 , Implying 

that the assuiiq>tlons used in obtaining the analytic value of F may be 

slightly In error. To provide a realistic estimate of thruster size 

requirements, the empirical value can be used to determine an effective 

individual aperture perveance <P > by setting 

8 


<p > 
a 



4 X 10~^ 
15,173 


« 2.6 X 10 



(5-6) 


Fur a hexagonal close-packed array of apertures, 


N - 0,9^, (5-7) 

where D Is the diameter of the array (beam diameter) , and S Is the 
center- to-center spacing of the apertures In the array. Eq, 5-3 
can now be written 


P “ 


2.34 X 10 



(5-8) 


and 


I 


b 


-9 D 

2.34 X 10 ^ 

S 



(5-9) 


Considering tue extensive empirical parameter variation that has produced 
the EM Ion optics design, the aperture parameters listed in Figure 5-1 
are probably near optimum and are also a realistic basis for predicting 
beam diameter requirements. Therefore, Eq. 5-9 becomes 

- 4.79 X lO"® . (5-10) 
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Eq. 5-10 can now be used to ku v 

for a tl,rust«r operated on mercury ee e flcur, Tr“ 

totel entrection volc<«e. stnce It le tb b "" 

apecmed. e turtber consddeLL de L o^: '"‘"T ““ 

•mplrleel resulte from EM testing elect k T' 

occur If *’ backstremtlng does not 


\ - Vj. • (5-11) 

Equation 5-ll la somewhat conservative In t-ha#- t-u 
R - 0.7 to 0.78 has been demonstrated without evld”^^^*^ operation with 

however, most testing hes been done under the condir 

Hence. E,. 5 -10 csn be written " "-“• 


^SO 7 

17 m U, f , 


II “ 8.19 X lO^^D^V 


(5 -.12) 


- 1.22 X 10^ T V ”3/2 

b b • (5-13) 

Using Eq. 5-13, the beam requiiements are illustrated 

-srcury In Figure 5-2 end for ergon In Figure 5-3 These 

bssls for determining thruster technology regulrewnts 1'“"“ 

tbrust end specific Impulse to hesm curL ^"^“‘'tncrr"? 

a good technology base for Ion optics designs up to 3ol dt ” 

Figures 5-2 end 5-3 Identlfv tb„ -Haaater. 

that regulre edvsnced telolot ZT 

P«ves„ce date obtained with sZe-otZaTlV" 

In the 0.5 to 4 A current r». -w 50-cm thruster Ion optics 

Oonscuently. the hZ Z 
as a first-order estate 
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BEAM CURRENT, A 


Figure S- 


. Beam diameter required to produce mercury Ion beam 
current . 




beam current, a 


Figure 5-3. 


S curSnT “8"” ion 



Curved Electrode Thrust Losses 


Since thruster diameters other than 30 cm are being considered in 
this study, the impact of grid curvature on thrust losses must be con- 
sidered. Thrust loss was calculated as a function of thruster diameter, 
grid spacing, and radius of curvature (dish depth) for various dished 
grid geometries. For all these calculations, the compensation was limited 
to provide a maximum vectoring angle of 16». Experimental results indi- 
cate that accelerator grid direct ion Interception increases rapidly for 
angles greiiter than 16 Although a relatively large dish depth (small 
radius of curvature) is desirable from a structural standpoint, the 
opposite is true from a thrust loss viewpoint. The results presented 
here show quantitatively the trade-off between thrust loss and dish 
depth (radius of curvature) as a function of thruster diameter. To all 
the thrust loss calculations presented here, the intrinsic beamlet 
divergence, which typically amounts to 1 to 2 percent, must be added. 

e* Grid Compensation 

Since thrust loss is a function of grid compensation, this 
parameter is discussed before the thrust loss calculations are made. 

Figures 5-4 and 5-5 show the geometrical parameters associated with a 
dished grid set. The basic condensation problem arises from the need 
to cancel (1) the effects of the dish angle and (2) the vectoring 
caused by the misalignment of the holes for noncompensated optics that 
occurs when the grids are separated. In the ideal case, the screen 
electrode hole pattern is "shrunk" an amount (compensated) such that all 
bea-nlets are vectored by an angle sufficient to exit parallel to the 
thruster axis, and, in that ideal case, the total compensation C 
required to accomplish this can be broken into two parts : 

^tot “ • (5-14) 

Neglecting the electrode thicknesses, the compensation C(£g) 
required to align the holes is given by 

- (2«gV tan , (5_15) 
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Figure 5-5. 


Maximum dish angle versus 


thruster diameter. 


% 
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where Is the grid spacing, D/2 is the beam radius, and is the 

dish semi-angle at the grid edge. For the present 30-em di^d grid, 
Ag - 0.056 cm, D - 29 cm, • 16% and therefore 

^ (16-) = O.m *. (5,16 

Which amounts to ^<25 percent the total 0,4 percent compensation 
presently used. If thermal loading causes the grid spacing to decrease 
by causing the screen grid to move toward the accel grid, the compensa- 
tion available for vectoring will be increased. 

The compensation C(a) required to vector the beamlets an 
amount a depends on the vectoring sensitivity measured in degrees per 
normalized screen aperture offset. There have been efforts to deter- 
mine this vectoring sensitivity using digital computer trajectory 
calculations, but the computed and experimental results did not agree 
Very well. The basic problem is that the computer programs must be run 
with planar symmetry although the real problem has no symmetry. The 
approach taken here is to use the experimentally measured sensitivity 
for the present 30-cm grid set. Figure 5-6 shows some experimental 
data taken with the coJllmated ExB probe. This data shows that there 
is a 2.4* difference between the r - 0 and r - 10.7 cm profiles. If 
the grid comrenaation were perfect, these profiles would be on top of 
each other. Since the dish semi-angle o at this radius is M2. 4", the 
net vectoring achieved was 

For a grid conpensafclon of 0.5 >ercent» the net grid offoet at thla 
radius is 

4^j. - C.00S) (10.7) - 0,06 tan 12.4* - 0.04 cm ^ 

The measured vectoring sensitivity Is thus 

10 * 

" TbV cm " 250“ /cm , 


263 



Figure 5-6. Beam profiles measured with ExB analyzer; 30-cm 
thruster using 0.5% compensation (Ref* 59). 
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or. in not«aU«d fom. since the screen enei-tar. « . 

®nora -IVaapercent ot screen spertnre dlsneter oflI!!t” '''^*'** 

iz7nT:::r.z:T ^ - “• 

««t.d . sensitivitp - 

described in the folloslnn aecti .c * “Iculetlons 

30^ EHI „ss nsed! ’ -aured »lth the 

**“ l oss Caleiilatlanw 

These calculations were made to aid in ^ 
of non, thruster design concepts. Set rei 2.ZT 

ta aaking the thrust loss celculstlons: » “"“aP^aao -are used 

• Ih. current density scross the grid varies ss 


■JCr) 


J cos 
o 2r 


A naxlmum beam vectoring of 16* at tbe be«. edge 

sensitivity -1V2.1 percent/ 

Using the first assumption rho a.i. 

^ICSS <*" paroant) is glv^ ^ thrust vector cos 6 snd thrust 


cos 0 


where 


i° cos cos e(r) dr/^ cos dr . ( 5 - 17 ) 


(5-18) 


^loas “ “ cos 0 J (X) \ 

and e(r) is the difference between the dlab angle a and the t- 
angle. The net angle 0fr^ n«os< 4 , "8 a and the net vectoring 

ia Che dlfferencetL 

Y(r) = ^ '"t ”« vectoring angle 


S(r) « a(r) - Y(r) , 


(5-19) 
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where 


Y(r) ■ “ d tan o] , 

A BKlBum vectoring angle o( 16“ lo eeeumed to prevent direct Inter- 
ception BO that the C<o) conpenaatlon amounta to a maxlmuin of 16 x 
2.1 percent - 33.6 percent of the acreen apertore diameter at the beam 

edge. Onder this assumption, the compensation used for the thrust lose 
calculation was 


■’tot 


■H 


_ + tan i}i 
® max 




(5.20) 

where 0 » 16® if a m „ 

*max 6*6 If a <lfi® Tu-t 

is plotted in t expression 

in Figure 5-7 as a f™.ctlon of thruster diameter with r.dl. . 
o curvature as a parameter. The present 30 

Increased to n » P«eent 30-cm compensation should he 

ncreased to 0.6 percent to be perfectly compensated. 

intrJsVr r* ’ ° 

tatrlnslc beamlet divergence, which typically amounts to 1 or 2 percent) 
e amount of thrust loss calculated using Eqs. 5.16, 5.17.and 5.19 
e P otted In Figure 5-8 as a function of thrust diameter with grid 
radius of curvature as a parameter. Figure 5-8 shows that the present 
20-ln grid curvature could be used up to a thruster diameter of .50 cm 
th less than 1 percent thrust loss. Figures 5-4 and 5-5 show that 
this diameter corresponds to a maximum dish angle of .30* and dish depth 


practic al Considerationn 

the calculations discussed above show the relationships 
between beam divergence, thrust loss, and grid compensation parameters 
Ideally spaced and aligned grids. In practice. Ion optics elec- 
trodes cannot be spaced completely uniformly, and alignment to obtain 
the pet required Offset In all radial directions presents a very 
moult «sembly task. If -.brusters are required In very large nus. 

. requirements for such extremely high accuracies represent a cost! 
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COSINE CURRENT DENSITY DISTRIBUTION 
16 DEG MAXIMUM VECTORING ANGLE 
GRID COMPENSATION SHOWN IN FIG. B-7 



25 BO 

THRUSTER DIAMETER. CM 

Figure 5-g. Ihruet loss versus thtusi 
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and formidable manufacturing obligation. Since 0.4 percent is 
considered the empirically determined practical limit tor 30-cm 
thrusters, the calculations given above are only useful in a relative 
sense and show what could be achieved Ideally. Achieving relatively 
small thrust losses would require developing fabrication processes to 
achieve adequate aperture alignment and electrode spacing. 


B. CATMODES 

The hollow cathode configuration shown in Figure 5-9 is the 
accepted cathode design for 30-cm thrusters. This cathode geometry 
as evolved from extensive development work at HASA toRC and elsewhere 
and a complete referenced discussion of this development has been pre- 
sented by Kaufman. Essential features of this cathode design are 

• A hollow cube with a gas-flow-llmitlng orifice 

M lM0*r^ temperature 

• snrf«r°^ “““ a low-work-function 

tonJjrC »Parating 

• Operational inodes with characteristics that donond n 

and emission conditions. that depend on gas flow 


Altuough several possible explanations of hollow-cathode physical 

processes are considered in Kaufman-s review article, the discussion is 

by no means conclusive. Some conclusions tor which there is general 
agreement are that 


deLltf °n characterized by an electron 

dis trihut 1 ^* t-h of 10l8/m3 with a Maxwellian energy 
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4 » 77-13 


SWAGED HEATER> 


BARIUM AI.UMINATE IMPREGNATED 
POROUS TUNGSTEN "INSERT" 



KEEPER ELECTRODE' 


Figure 5-9, Hollow cathode configuration used in 30-cin EMT. 
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Cathode orifice wear is proportional to emission current 
reqi^irements and orifice dimensions. 

The nominal emission current rating of the 30-cm thruster hollow 
cathode is 12 A for 2 A of beam current. Life tests have been conducted at 
this current that show negligible cathode wear for 18,000 hr of testing. 
9iuce the 30-cm thruster cathode was scaled from the 15-cm thruster cath- 
ode (SERT II) on the basis of empirical results, we can suggest with some 
confidence that scaling a cathode for higher emission currents should not 
present a basic technology problem. In previous sections, beam currents 
as high as 15 A were considered a possibility, and such beam currents 
imply emissloh currents of 90 A. Emission currents of 20 A have been 
achieved without difficulty from the present 30-cm cathode design, 
but extrapolation to higher currents would require experimental verifi- 
cation and probably some additional engineering to achieve the thermal 
environment considered to be a requirement for long lifetime. 

The discussion presented above is relevant only for operation with 
mercury vapor as the propellant or working gas. Achieving the requisite 
low-voltage discharge conditions has not been adequately demonstrated 
With other gases, is not well documented empirically, and was not a sub- 
ject of this Investigation. Empirically, xenon appears most similar to 
mercury vapor with respect to discharge characteristics and cathode 
operation. Consequently, for operation with propellant other than mer- 
cury, a basic technology investigation and development program is 
considered to be a requirement. 
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C. THERMAL LIMITS 

As discharge power is increased to obtain higher beam currents, 
the temperature distribution throughout the thruster changes. Power is 
not distributed uniformly, and some parts of the thruster will be heated 
more than others. To assess the relationship of thermal limits to beam 
current, a thermal modeling analysis was performed. An 11-node thermal 
model for the 30-cm thruster has been developed and has been used to 
predict the temperatures of critical thruster components for high beam 
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current operation. The starting point for this effort was a thermal 
model developed by Oglebay at NASA LeRC^° for the 30-cm thruster. The 
nodal conductances, view factors, and power distribution for this model 
were received from NASA LeRC and were Incorporated where applicable in 
the development of this simplified thermal model. A description and the 
results of this model are presented below. 


1. Thermal Model 


The thermal analysis was carried out by representing the 30-cm 
thruster as a lumped parameter network of 11 nodes with 22 radiation and 
conductor resistances. A schematic showing the thermal model and nodal 
numbering is shown in Figure 5-10. Each of the principal thruster coim- 
ponents was represented by a node that was connected to adjacent nodes 
by radiation and conductance resistors. These parameters., along with the 
nodal area, view factor, and the emissivlty assumed,, are summarized in 
Table 5-1. The radiation exchange factors were calculated using 


^/ij 




Vij 


-1 


(5-23) 


and are also shown in Table 5-1. 


2. Power Inputs 


The nodal power inputs were calculated as a function of beam cur- 
rent on the basis of a discharge chamber efficiency of 185 eV/ion and a 
fixed loss of 21 W. The actual effective heating power produced by the 
discharge chamber efficiency was assumed to be reduced by (37 + 10) 
eV/ion to account for the kinetic energy and heat of ionization 
taken away by the beam ions. This results in an effective discharge 
heating power as a function of beam current of 
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Table 5-1. Suntmary of Nodal Areas, View Factors, Emissivltles 
and Form Factors Used in Thermal Model 


Nodes 

Connected 


2 and 1 

3 and 2 

3 and 6 

5 and 7 

6 and 7 

8 and 11 
10 and 11 

9 and 11 

4 and 11 

3 and 11 

4 and 3 
4 and 9 
4 and 10 

4 and 5 

5 and 6 

5 and 8 
9 and 10 

2 and 6 

6 and 8 

3 and 5 
5 and 9 
3 and 9 
3 and 10 


Connector 

Number 


231 

120 

120 

352 

396 

81.1 


195 

231 

74.4 

396 

396 

81.1 


111 

195 

91 

195 

171 

120 

155 

120 

171 

120 

171 

114 

171 

113 

171 

352 

215 

396 

215 

81. 

114 

114 


I 0.15 
0.20 
0.15 
0.15 
0.15 
0.15 
0.15 
0.15 
0.25 
0.25 
0.25 
0. 25 
0.25 
0.25 
0.15 
0.15 
0.15 


1 

0.20 

0.15 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 


0.5 
0.47 
0.03 
0. 175 
0.50 
0.50 
0.50 
0.235 
0.05 
0.16 
0.17 
0.06 
0.06 
0.50 
0.325 
0.02 
0.265 


0.25 0.15 0.053 
0.25 0.15 0.056 


A^F^, 
in. 2 


30.1 

14.6 
2.5 

30. 9 

51.6 

10.6 

14.5 

9.2 

7.4 

16.8 

10.1 

6.07 

6.05 

22.05 
18.2 

3.0 

7.5 

0,05a 
0.05a 
0.095® 
0. 119® 

5.06 
5.27 


Conductances in W/‘’C 



where Is in amperes , and P i® j 

to the screen, anode, and hecSlste nodls ,Tf iT” ■‘^=»«“ted 

oHectlve screen grid open area ol 70 percent the’”'l.T“""' “ 

which strike this electrode orod ’ 

ottrode produce an effective heating power of 

/ . _ V 


sc " 


(“0. 70^°) [^7 + loj 1 . 20 , 


- f ^ j b ’ (5-25) 

Where we have neglected the cooling effect of m, , 
with ion wall flux. electrons that combine 

By assuming that all the fixed losses are disc .k 
P late electrode, the power to this electrode was taLn aT ^ 

^back * IB I + 21 

. ^ (5-726) 

where the coefficient nf t 

2-A beam current. The heat to thTlnodl OSlebay’s model at a 

EOS. 5-25 and 5-26 from Eg, 5-24 ras in ““ 

q ^ -J4, resulting in an anode power 


anode * • 


sZZZ 6““"°" -e haam Ions 


(5-di7) 


r 

anode 

P 

screen 

p 

back 

P, 

lost 

Total: 


105 r 

D 

20 

13 + 21 

47 L 


185 Ijj + 21. 


The accelerator grid power (I y '> 
of beam current by linearly extranortl^*^ calculated as a function 
accelerator current at 1 and 2 A h ^ typically measured 

was found by solving: current. The total voltage 

I, » P V 3/2 
b X t 


2.75 


as a function of beam current and uslns a tot«i 

g a total thruster perveance of 

P - 4.3 X ^ 

which is typical of present close spaced optics A ^ i 

the resulting nodal powers as a function of b ’ 

Table 5-2. tunctlon of bean, current la shown In 

Results and PlscusiafAT' 

the calculated temperatures as a functl . 

2 ?: 

r- - =.teeme:: -h^Vca^jra”^^ , 

temperatur. are plotted In figure 5-11 as a function of'be“ cur 

7o estimate th** 1 u ‘-wt-rent. 

30-cn. thruster, we will ass IT the present 

than 500-C and that the vaoori temperature must be less 

we will also assl bat I 

temperature and that the lower lln.lt to vaporiser temnert ' 

the rear shield temperature Ihd.v rn temperature would be 

^able 5-3 or figure 5-1 rthlt T bi b from 

tiiG highest besm current io c a 

^intaln the radial magnet temperature at less thsn“oo“0 . “ 

3:i:::rci“ “ ~ ^ 

- o. .aporlsers :ral: r 

(-taen grid,.- if an uppeT lllt o T 

the 30-cm thruster could probably be operated up to 15 ! ’ 

Clearly, the maximum beam cur^en^ 
hy a maslmum operating te.^erature (probably Ton " 
thermal analysis discussed previously were used tlTstTbilsh'r^ " 
emperature-lln,ij.ed beam current, as shown in Figure 5-12 si 
-action Of the radiated heat is rejected through the ion’opti::; Ls"" 
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rt Table 5-2. Summary of Nodal Input Powers Used In 

[ Thermal Model 













Table 5-3. Summary of Calculated Temperatures for 30-cm 
Thruster. Values in ( ) Were Determined 

Experimentally . 




1 

Calculated Temperature, 

“C 

] 

Node 


Beam Current, A 


1 

2 

4 

6 

8 

10 

Rear shield 

2 

(120) 

97 

(153) 

143 

206 

251 

288 

319 

Base plate 

3 


(262) 

258 

(328) 

348 

411 

462 

504 

Anode 

4 

271 

372 

492 

574 

638 

691 

Body 

5 

(174) 

152 

(230) 

224 

312 

372 

420 

459 

Ground Sc. 

6 

47 

(96) 

96 

160 

204 

239 

268 

Mask 

8 

33 

(62) 

79 

137 

178 

210 

236 

Screen 

9 

202 

293 

404 

480 

540 

589 

Accel 

10 

163 

246 

350 

426 

487 

538 


6119 
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Figure 5-11. Calculated temperatures versus I, . 
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BEAM 



MAXIMU^r BEAM CURRENT, A 


Figure 5-12. Estimated maximum beam current (limited by 
temperature) versus thruster diameter. 
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curve was derived by assuming that radiated power is proportional to 
beam diameter and beam current in the form 

" ^2 (^b) max ’ (5-23) 

Therefore » 

^ " '‘3(M"iax (5-29) 

max 

{^b)max " ^4® * (5-30) 

where the K's are constants. 

^max * maximum allowable teraperature 

D = beam diameter 

^b max maximum allowable beam current 

^4 " 0.0167 A/cm^ . 


The constant In E,. S-30 was avaluatad for tha 30-c., thruatar aasnnlng 
a maximum current of 15 A. 


0. WEAROUT LIFE 

In tha pracadlng dlacuaslons, thruatar dlamatar, or aora ptoparly 
ion baam diamatar, was datarwlnad as a function of spaaitiad baaa cur- 
rant (thrust) and ™itaga (I,p) f„u both marcury and arion propallants 
on tha basis of stata-of-tha-art ion axtraatlon systam capacity only. 
It is now appropriate to axawina thasa rasults with raspact to those 

practical constraints that are presently known and others that can be 
anticipated. 


Consider first the operation of a 
voltage to obtain higher beam current. 


30-cm thruster at higher beam 
With regard to tbe ion extraction 


ystem, results presented ... cnar operation at 

incraasad baaw voltage has ralatlvaly little Iwpact on tha design of 
this thruster coeponant. If beat, current is kept within the parvaance 
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a“clT : T 

ctlng the Ion oytleal propertlee o£ the extraction eyeten. Slmllar- 
y. operation et increased beam voltage Is o£ little or no consequence 
to either discharge chamber or cathode operation or design. Propellant 
lectrlcal Isolators must be adapted to the appropriate voltage Insula- 
tion. as must all the Insulation used, but this Is not a significant 
constraint unless extremely high voltages are specified (>10 kv). 

Increeslng the beam current at constant, thruster diameter, homever 
requires an Increase In beam current density and thus a proportional 
-ncrease In plasma density and discharge current. Increasing plasma 
density at an' approximately constant Ionisation fraction has been found 
ncrease the double Ionisation percentage, and thus discharge chamber 
erosion rates Increase at a rate somewhat greater than linear with 
increasing discharge current (beam current), dn approximate dependence 

partlcl r Of 

Lt Tu u ^ ™ 

that surface can be expressed as 


1 - + 


(5-31) 
(5-32) 

Where S la the sputtering coefficient for the appropriate Ion species. 

d e is the electronic charge. The fraction is a function of 

the total current dentisy i On ^ ^ ^ 

with an EMT 7 nn empirical data obtained 

With an EMT-like 700-series 30-cm Hg thruster. 


'-H- 


» 0.1 + 


40 




(5-33) 


where I Is the total beam current, and the current densities are average 

va ues for the screen grid boundary. B> setting proportional to total 
beam current, Eq. 5-32 becomes 






(5-34) 
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whsrs C is the proportionality constant between current density and total 
current. The sputtering coefficient for doubly charged ions is assumed 
to have the value measured for singly charged lons^^ at twice the ion 
energy (voltage). For the energy range of interest (30 to 40 V), the 
ratio S. 


1 ./®+ approximately 15. Hence, Eq. 5-34 can be written 


q « K I, 




K 1^(0. 2 + 1.8) , 


(5-35) 


where K includes the appropriate considerations for relating current 
density to beam current, sputtering coefficient for the material being 
eroded, and to the energy of the eroding ions. 

If a particular discharge chamber surface is now selected as the 
critical surface for material loss, it is possible to use Eq. 5-35 to 
relate useful lifetime to beam current. By specifying the allowable 
material loss in terms of atoms per unit area N, the lifetime, x, is 
given by 


q T 


N 


constant , 


(5-36) 


and hence 


q^Tj = ^ 2 ^ 2 * » 

where the subscripts denote lifetimes and erosion fluxes for different 
current values; therefore, 

q^KIj (0.2 Ij + 1.8) = T 2 KI 2 (0.2 + 1.8) 

or 

1.8 + 0.2Ij 
1.8 + 0.21^ 

Eq. 5-37 can now be used to extrapolate measured thruster lifetime for 
operation at lower currents to operation at higher current under the 
assumption that the conditions (particularly constant discharge voltage) 
stated in Eq. 5-33 are valid. Thus, conditions are thought to be satis- 
fied for a 30-cm thruster of the 700-seriea EM design if the discharge 


) 


(5-37) 
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power Input in watts per beam ampere (eV/ion) and the true propellant 
utilization efficiency are held constant. As an example, the screen 
grid lifetime was 60 percent consumed in the 10,000 hr endurance test 
of thruster S/N 701 (NASA CR-135011) and the average current was 
1,4 The lifetime T at current can be estimated by 


t(I) 


4.85 X 10^ 
r (1.8 + 0.2 I) 

D 


(5-38) 


for operation at a discharge power input of 185 eV/ion and a corrected 
discharge propellant utilization efficiency of 88 percent (as shown 
graphically in Figure 5-13). Although the applicability of this curve 
to thruster operation in general may be quantitatively in question, this 
type of relationship between lifetime and beam current is considered to 
be qualitatively valid. 

The ion extraction screen grid has been selected as the critical 
surface here because, as a discharge chamber component, the ion extrac- 
tion screen grid thickness affects the efficiency of producing beam ions 
(eV/ion). Consequently, if one wishes to preserve thruster efficiency, 
one cannot simply Increase the grid thickness to increase its lifetime, 
as is possible with some other discharge chamber components. Although 
the absolute value of the relationship shown in Figure 5-13 may be 
somewhat uncertain, it is possible to make the observation that, for 
short-term missions, existing 30-cm technology could produce high beam 


currents with useful lifetimes. 

It is anticipated that operation at lower propellant utilization 
efficiency will change Eqs. 5-33 and 5-38 such that a family of 
curves will be generated that arc shifted to higher lifetime. Again, 
drawing on empirical data, the relationship between 
rected discharge chamber propellant utilization efficiency, shown in 
Figure 5-14, can be used to modify Eqs. 5-33 and 5-38 to obtain 
Figure 5-15. The applicability of the empirical data in Figure 5-14 in 
a general consideration la quantitatively questionable; however, the 
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ngiiro 5-15. Estimated thruster wearout versus beam 
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trends prodiitsd (ss she™ In Figure 5-15)should be valid. Consequently, 
the observation that a relatively smll change in propellant utilization 
C5 percent) can produce a large change in screen grid wearout llietine 

(100 percent) is considered valid (except for the absolute values of 
the changes). 

It is useful at this point to consider the constraints placed on 
thruster design by the implications of figure 5-15. If discharge chamber 
wearout lifetime is determined by the ion extraction screen thickness, 
unless the thruster configuration and operating characteristics are 
significantly changed, then Figure 5-15 can be used to restrict the 
relationship between beam current and beam diameter. Supposing that the 
wearout lifetimes shown in Figure 5-15 may be as much as a factor of two 
too low, the maximum beam current that can be considered practical for a 
minimal lifetime (3000 hr) is still only about 10 A. This limit is 
essentially independent of whether the propellant is argon or mercury, 
because the sputtering coefficients are essentially the same for both 
elements and doubly charged ions are formed in argon about as readily 
as in mercury. Extension of this limit would require new. advanced 
technology that is not directly identifiable from vjnsideratlon of 
present 30-cm thruster development. Figure 5-16 shows the modifications 
of the beam-current/beam-diameter relationship with the current density 
limitation applied. The shaded areas of the charts are unavailable 
because the current densities required would not be consistent with 
lifetimes considered to be useful. The information in these relation- 
ships is applied as Illustrated in the following example. If a given 
mission analysis Indicated that a 40 A beam of argon ions at 2 kV beam 
voltage would be optimal for each thruster module, then Figure 5-16 
indicates that the minimum beam diameter for useful lifetime would be 
60 cm even though the ion optics could be designed to produce 40 A at 
2 kV with a 50- cm- diameter beam. 

E. THRUSTER SCALING 

In the preceding sections, the factors that limit the power and 
thrust capacity of existing 30-cra mercury ion thrusters were discussed. 
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It was shown that the upper ll:.it on heae current Is 10 or 15 A depending 
on whether wearont lifetime or therwal constraints are considered to be 
^st Iwpurtant. For a Blsalon of 125 days or lore, screen grid wearout 
fetloe dominates and the current limit Is 10 A. One option Is to con- 
sider using several thruster modules with a single power processor to 
Obtain the tegulred lifetime. Ihe Ion optics thermal limit then becomes 
most important because any attempt to alleviate thermal problems In this 
component would constitute a major development task, whereas magnet 
temperature limitations might be alleviated by changes in magnet material 
or reduction of heat transfer to the magnets. Hence, 15 A would be the 
upper operational limit for a 30-cm-dlameter thruster, based on the 
assumed Ion optics thermal considerations. This limit is Independent of 
propellant type or specific Impulse since the discharge power Input 
requ red per unit of beam current Is essentially a constant and would be 
appromlmately 3 kH. Consequently, any Improvement In 30 -cb thruster 
capacity would have to address the problem of rejecting this order of 
thermal Input without the temperature of the accelerating electrodes 
exceeding lOO'C. An advanced concept 30-cm thruster would therefore 
upend on development of new technology capable of higher power density 
operation and designed specifically for heat rejection. Some of the 

fundamental requirements that can be predicted for specifying such u 
thruster concept are 


• High specific Impulse operation to permit adequate electrode 
thickness for heat rejection by conduction 

• Improved beam and discharge plasma uniformity to distribute 
thermal Inputs more evenly 

• Incorporation of materials capable of higher temperature 
operation 

• Efficient operation at low discharge voltage to minimize ion 
sputtering erosion. 

Increasing the thruster beam area while keeping power (or thrust) 
densities at or below the limitations Indicated above can be approached 
as an engineering problem rather than one requiring technology. As in 
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determining thruster limitations, the appropriate starting point in 
considering larger thruster modules is to examine the ion extraction 
system requirement • 

One way to increase the beam area would be to increase the beam 
diameter. To maintain high current density at low beam voltage, the 
dimensions of the electrode thickness and interelecfcrode spacing must be 
kept as small as practical. This dimensional limitation Implies the use 
of curved or dished electrodes to achieve mechanical stability. As 
discussed in Section 5. A. 2, electrode curvature introduces a thrust loss 
because of non-axlal beam trajectories. The technique used to vector the 
off-axis bearalets to ’'compensate" for electrode curvature approaches the 
direct interception limit at 30-cm beam diameter and therefore larger 
diameter thrusters might have to operate with uncompensated apertures 
at the periphery. Thrust losses resulting from this limitation were 
discussed in Section 5.A.2.b. An additional constraint to increasing 
beam diameter is the availability of electrode material. Molybdenum 
sheet is currently only avllable in widths up to 24 in. (61 cm), and 
at present no other material offers a viable substitute. Therefore, 

50 cm is at present the practical fabrication limit on beam diameter. 

F. PROPELLANT-TYPE CONSIDERATIONS 

The 30-cm EMT has been developed primarily for operation using 
mere../ as the propellant. It is a relatively straightforward matter to 
operate a thruster on other gases since the discharge mechanisms are 
not dependent on the mass or ionization potential of the gas used. 
Performance characteristics are another matter, however, and the highly 
attractive characteristics of ion thrusters based on operation with 
mercury may not be maintained with an arbitrary choice of propellant. 
Because of the relatively large quantity of propellant required for 
some large earth orbital satellite applications (-v 10 kg), an inert gas 
(eg., argon) appears to be an attractive propellent possibility. The 
first consideration is the ionization potential: 15.75 V for argon as 
compared to 10.4 V for mercury. This means that more discharge energy 
goes into ionization of each argon ion and also that discharge v iltage 
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must be higher to obtain the same number o£ ionizations per primary 
electron (cathode emitted). These facts present a trade-off between 
operating at increased discharge voltage (thereby increasing discharge 
chamber ion sputtering erosion) and limiting discharge voltage (thereby 
reducing propellant efficiency and Increasing charge exchange ion sput- 
tering erosion of the accel electrode). Thi«! trade in erosion rates 
will be assesed in more detail at a later time. 

The mass of the argon atom is only 20 percent of the atomic mass of 
mercury, and therefore operation on argon implies higher specific 
Impulse or higher power per unit of thrust. To approach the 30-cm EtfT 
thrust and lifetime specifications (for mercury operation), an EMT design 
thruster operated on argon would have to be operated at a specific 
Impulse of about 8000 sec and at a 12-kW power level (5500-V beam volt- 
age at 2-A beam current). This power requirement could be prohibitive 
in some applications. 

Next we consider some of the advantages and operational implications 
of using argon as the propellant. First, argon is a gas under anticipated 
operating conditions and no heaters or vaporizers are required for pro- 
pellant supply. Flow rates can be accurately controlled using flow 
Impedances and constant pressure regulators. This implies that thruster 
control must be established using only electrical parameters (not pro- 
pellant flow) unless a variable flow Impedance is developed. Thruster 
turn-on time can be essentially eliminated since no vaporizer warm-up 
is required. Because of reduced ionization efficiency, double ioniza- 
tion processes are not anticipated to be as important to thruster dis- 
charge chamber life in argon thrusters as in mercury thrusters. Thus, 
it may be possible to Increase propellant utilization efficiency by 
increasing discharge voltage and yet maintain constant wearout life. 

The more significant Implication tc advanced thruster technology, 
in so far as argon operation is conceimed, is thruster lifetime and con- 
trol. In applications where argon is nearly mandatory because of the 
quantity of propellant required, the number of thrusters will be large - 
even if the module size becomes significantly larger. Moreover, it 
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would be extremely desirable for an ion thruster to turn on and off with 
the simplicity of a fluorescent lighting fixture and to be capable of 
operating as long, or longer. Current hollow cathode technology provides 
long operating life but is relatively complex from the standpoint of 
cathode Ignition and control. Advanced cathode technology has been 
investigated by personnel at NASA LeRC, Colorado State University, and 
HRL, but no results have been reported, at this time, that show possi- 
bilities of simplified long-life thruster operation on argon or other 
gaseous propellants. Thus, there is a clear opportunity for thruster/ 
cathode technology innovation. 

G. OPERATIONAL SIMPLIFICATIONS 

The 30-cm mercury ion thruster, as presently developed, requires 12 
power supplies, at least 8 of which are required to be more or less con- 
tinuously variable in output. The 12 power supplies are 

Variable output power supplies 

• Beam or screen supply 

«i Discharge supply 

• Accel supply 

• Magnetic baffle coil supply 

• Cathode vaporizer supply 

• Main vaporizer supply 

• Neutralizer vaporizer supply 

• Neutralizer keeper supply 

Fixed output power supplies 

• Cathode tip heater supply 

• Neutralizer cathode tip heater supply 
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Isolator heatev^ supply 


• Cathode keeper power supply 

Of the variable output power supplies, the only obvious possibility for 
simplification Is the combination of beam and accelerator supplies. 
Feasibility of this combination was demonstrated under the 2.5-kW 
Advanced Technology Ion Thruster program (contract NAS 3-17831) during 
1975, ' Similarly, all of the fixed power supplies listed, except the 
cathode keeper supply, are used only during start-up (or for very low 
power operation) . Since the discharge power supply is variable and not 
in use during start-up until the heaters are no longer necessary, it can 
be used to supply heater power by appropriate heater design and inter- 
connections of heaters. The power processor being developed for the 
30-cm EMT uses the discharge supply for powering the Isolator heaters 
prior to discharge ignition. This capability was also demonstrated 
during 1975 under the above-mentioned program for the cathode tip and 
the Isolator heaters. Similarly, the cathode keeper discharge can be 
maintained after keeper Ignition by using a dropping resistor from the 
anode, thereby eliminating the steady-state keeper supply requirement. 
The stability of the keeper discharge In nhls mode of operation has also 
been demonstrated fur all anticipated transient conditions. Thus, It 
would be possible to eliminate* one variable output and three fixed output 
power supplies without altering any of the steady -state control concepts 
now in use. 

To Identify further simplifications, it is appropriate to examine 
the function of each of the eight remaining power supplies as related to 
thruster operation requirements. 

1, Beam (or Screen) Power Supply 

This supply provides the voltage necessary to accelerate and focus 
the ion beam. It Is operated as a constant -volt age supply with a cur- 
rent limit. Beam current Is determined by discharge and vaporizer 
parameters. 


2. Discharge Power Supply 


This power supply ionizes the propellant and is operated as a 
constant current supply. Discharge current is an independent param- * 
eter determined by the beam current desired. Cathode propellant flow 
determines the discharge voltage, 

3. Magnetic Baffle Coil Power Supply 

This power supply excites the magnetic baffle coll and is an effec- 
tive adjustment on the discharge volt age/ cathode propellant flow 
characteristic. Variable output for this power supply is a convenience 
or flexibility option. 

4. Cathode Vaporizer Heater Supply 

This power supply provides cathode vaporizer propellant flow and its 
output is controlled in direct proportion to discharge voltage (closed 
loop), as compared with a reference voltage (independent parameter). 

5. Main Vaporizer Heater Supply 

This power supply provides main vaporizer propellant flow and its 
output is controlled in proportion to beam current (closed loop) and 
compared with a reference signal (independent parameter). 

6' Neutralizer Keeper Power Supply 

This power supply initiates and maintains the neutralizer hollow 
cathode plasma. The power supply is operated as a constant-current 
source with keeper current as an independent variable. Keeper power 
supply voltage varies with neutralizer vaporizer propellant flow rate. 
The neutralizer keeper currt.-.t set-point is automatically changed when 
the beam current is turned ofj to maintain the total neutralizer emis- 
slon current essentially constant. 

Neutralizer Vapo rizer Power Supply 

This power supply provides neutralizer vaporizer propellant flow 
and is controlled in direct proportion to the keeper voltage Culoeed 
-OOP) in comparison to a reference voltage (Independent variable). 
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Neutraliz er Cathode Tip 

This power supply Is used only during start-up periods to obtain 
Initial neutraliser cathode Ignition. Combining this function with the 

Bsln cathode and Isolator heater presents an Insulation problem for 
Steady state operation. 

Finding further simplifications that would eliminate power supplies 
involues the defining of new control concepts or some loss of fleslblllry 
in parameter selection. For exsmple. operation at a single beam current 
could reduce the discharge supply from a variable to a fixed-output power 

sZl r "“I" vaporiser supply somewhat. 

Similarly, a thruster modification could eliminate the magnetic baffle 

lldTr™"' 

would eliminate vaporizer power supplies. 
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SECTION 6 

ADVANCED THRUSTER CONCEPTS 
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50 cm Is probably near the limit of high performance 
operation with a single cathode. At thla^sS" 
plasma nonuniformity with a large central cathode 
extensive experiments for magnetic 
field shaping (e.g,, cusped fields). 


A 50-cm thruster Is in the region of Interest 
large vehicle missions In Earth orbit. 


for 


The design shown in Figure 6-1 utilizes the 30-cm thruster struc- 
tural approach and the 30-cin thruster length. To a first approxlnatlon. 

e length need not be scaled since the neutral atom naan free path 
before ionization is less than the chamber length. Additional insulators 
and supports are provided for the optics and ground covers. With the 
two-point glmbal mount, the structure would probably require additional 

St <“‘18 ="<1 thickness compared to the 30-cm design since loads are 

carried further. In addition, the fabrication tolerances and assembly 
procedures would be tighter than for the 30 -cm thruster because of the 
greater difficulty of defining the ion beam axis relative to tbe gi-K.i 


The cathode and cathode pole shown are the 30-cm thruster design. 
The propellant plenum was arbitrarily placed at half the anode radius 

Clearly, the final detailed designs of these and other parts would 
require prototype testing. 


ADVANCED 50-CM THRUSTER CONCEPT 


The 50-cm thruster design discussed in the previous section has 
several fabrication and structural weaknesses. The concept Illustrated 
in Figure 6-2 should improve on these aspects of the design. Since a 
major portion of the thruster mass is associated with the optics 
assembly, use of the shell for structure is inefficient. By bridging 
around the shell at three locations, the optics loads can be transferred 
more directly to tbe glmbal mounts. In addition to improving the load 

path, the concept shown in Figure 6-2 simplifies the optlcs-glmbal 
tolerance problem. 
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Figure 6-2. Continued. 
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The 30-cin EMT shell design resulted from a series of changes 
eslgned to make the thruster adequate dynamically without major re- 
design. However, if the shell design task were Initiated anew, a 
concept more like that shown In Figure 6-2 would probably result 

Avoidance of complex jigging, as currently required, should reduce 
assembly time and cost. 

C. OVAL CROSS SECTION CONCEPT 

An alternative to increasing the Ion beam diameter Is to Increase 
the beam area by increasing one dimension only, thereby producing an 
oval cross section. An example of this approach Is Illustrated In 
Figures 6-3 and 6-4. The beam envelope becomes 30 cm x 90 cm and 
could be expected to perform at approximately three times the 30-cm 
thruster capacity with essentially the same lifetime and operational 
characteristics. Since molybdenum sheet Is available In any desired 
length (at 60 cm widths), fabrication of the electrodes becomes a 
question of designing the proper hydroforming fixtures. Several main 

problems can be anticipated In developing such an advanced concept 
thruster: 

• Achieving structural stability In the ion 
extraction electrodes and mounting assembly 

• Achieving a relatively uniform discharge 

plasma In the elongated geometry 

• Achieving adequate structural strength. 

The structural questions present an engineering problem that Is probably 
tractable If the weight constraints Imposed are not too severe. The 
plasma uniformity problem Is not so easily dismissed, however. The 
proposed solution to this problem makes use of two thruster cathodes 
and central propellant Introduction. Integration and control’ of two 
hollow cathodes as shown represents a development task thar has not been 
approached before! however, there should be no fundamental difficulty 
in such an approach. Solution of this problem will have to be primarily 
empirical because discharge modeling Is not sufficiently analytical to 
permit design computations. Parameters to be explored would include: 
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(a) DISCHARGE CHAMBER VIEW FROM FRONT 


Figure 6—3, Oval, 30 x 90— cm thruster design. 










• Baffle size and shape ' 

• Cathodt:: location 

• Pole piece geometry 

• Propellant introduction 

• Methods for controlling cathode 
propellant flow. 

In the configuration shown in Figure 6-3, the magnetic field line 
geometry is still the divergent or SERT II type. Improved plasma uni- 
formity has been obtained and reported by workers making use of boundary 
magnetic field geometries that have a cusped shape. A single cusp ' 
configuration, as shown in Figure 6-5, is the most readily adapted to 
existing technology. In this configuration, the plasma volume is not 
entirely free of magnetic field, but improved plasma uniformity has 
been demonstrated. A multiple cusped geometry, shown in Figure 6-6, 
introduces considerably more complexity. Since the latter geometry 

might be operable with a single cathode, however, it cannot be ruled out 
a priori. 


To e\amine some of the properties of an oval thruster in terms of 
EMT technology, two assumptions were made: 

• An approximate thruster weight can be estimated 
by assuming that the weight of discharge chamber 
elements can be scaled in proportion to the cir- 
cumference of the chaiTber, and that optics 
assembly and backplate weight scales in proportion 
to beam area. 

• A plasma and beam current distribution can be 
established across the shorter beam dimension 
that matches the riidial density distribution 
of a 30-cm thruster. 


Given these assumptions, the scaled thruster mass would be about 
2.4 times that o*f_ the 30-cm HSlT, and ttit beam area ratio would be 3 55 
times that of the»|^0-cfli^,EM!n If^one assumes a uniform beam current 
density, then the larger .j^hrus.ter' could provide 7 A beam current with a 
lifetime of greater than lu hr or a maximum beam current of 53 A for 
shorter lifetime. 
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Figure 6-6. Multiple cusp magnetic field discharge 
chamber configuration. 
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If one assumes a current distribution that has the form 


J(r) - cos f- , 


then the current capacity for the oblong thruster would -be about 4.5 
times that of the EMT, and, consequently, the thruster could provide 
9-A beam current with greater than a lO^-hr lifetime or a maximum beam 
current of 67 A for shorter periods. Thus, the scaled oval thruster 
configuration would have the following specifications; 


• Thruster mass 18.5 kg 


• Beam current 9 A (average) 

• Beam current 67 A (maximum) 

• Lifetime 10^ (A-hr) . 

Preliminary assessment of the thermal characteristics of the 
stretched" thruster shows no basis for estimating significant differ- 
ences since the primary thermal process, radiation through the screen 
grid, has been scaled directly by area. Consequently, providing that 
discharge current and plasma density distributions can be established 
as assumed above, the scaled thruster configurations would have com- 
ponent temperatures essentially the same as those of the EMT. For an 
EMT-type magnetic confinement geometry, this may be a rather pessi- 
mistic assunqitlon, and the cusped configurations may be more realistic 
for achieving a satisfactory plasma distribution in the oval geometry. 

The cusped magnetic geometries may offer more promise for 
establishing a satisfactorily uniform plasma distribution, but they will 
clso carry a weight penalty. We estimate that the additional anode and 
poleplece required for a single cusp geometry would add 1.2 kg (6%) to 
the thruster assembly. 

Assuming a periodicity every 2.5 cm for the multiple cusp 
geometry, 6.5 kg (35%) additional weight would be added to the thruster. 
Since the 6.5 kg increase in weight could require heavier structural 
elements to satisfy launch requirements, this is probably a low estimate, 
Consequently, the single cusp may be the most promising aLtei^ative. 



The most critical single element in this ’’stretched'* configuration 
design is considered to be the ion optics assembly. The concerns are 
more practical than fundamental In that a fabrication process is 
required for manufacturing electrodes in the oblong shape with 
sufficient conformality that uniform Interelectrode spacings can be 
achieved. Conceivably, the processes now used for hydroforming and 
chemically milling the circular, dlshed-grid electrodes can be adapted; 
however, this supposition remains to be demonstrated. Similarly, the 
electrode mounting structure becomes more vulnerable to distortion in 
the Oblong configuration, and that possibility also affects the ability 
to achieve satisfactory interelectrode spacing. These two elements are 
the major concerns for increasing the longer dimension of the "stretched" 
thruster. If low 1^^ (close grid spacing) is not a requirement, then 
the degree of difficulty is lessened. If these unproven factors can be 
adequately handled, the long dimension of the oval thruster can be 
extended or shortened to match thruster module requirements with a mini- 
mal impact on thruster design and performance characteristics. 

An additional factor that should be mentioned is reliability. 
Although the oval concept may require multiple cathodes for best 
performance, operation with a failed cathode would be possible. 
Performance would probably be degraded with a failed cathode (e.g,, 
one failed out of three) , but mission options should be better relative 

to having a completely failed module. 
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SECTION 7 


STUDY CONCLUSIONS 

The most significant products of this work are the cost modeling 
methodology developed for Earth orbit missions, the generalized analysis 
approach for planetary missions, and the technology assessment results. 
Specific conclusions relative to the selection of an advanced thruster 
concept were also developed. These specific conclusions are Important 
and necessary, but their value depends greatly on the input assumptions. 
Although the assumptions and resulting conclusions may be changed by 
additional study, the analysis tools developed here should provide a 
straightforward means of evaluating the assumption/conclusion relationship. 

Earth orbit mission ^ost modeling results produced several reason- 
ably clear conclusions. Under most conditions, module power is the key 
thruster parameter in producing low-coat transportation and on-orbit 
propulsion. The upper bound on module power is not a hard boundary, but 
thruster development considerations will probably restrict the power 
level to about 100 kW. Depending on the mission, thruster wearout life 
is the next most important parameter. When refurbishment is not included, 
wearout life is quite significant in the overall system cost picture. 

Other thruster parameters initially thought to be important (e.g., mass, 
unit cost, redundancy) appear to have or ly second-order effects on total 
cost. Power processor and power source sarameters strongly influence 

total cost and should be carefully considered in selecting a vehicle 
design. 

After evaluating the cost sensitivity of the transportation and 
on-orbit missions to a wide range of parameters, a 50-cm-dlameter 
thruster was suggested. This thruster, operated at a beam voltage of 
about 2400 V, would satisfy a wide range of Earth orbit missions, For 
large systems (SPS size), the 50-cm thruster could be operated on argon 
with a beam power of about 100 kW. For somewhat smaller vehicles, xenon 
could be used effectively with a beam power of about 50 kW. \ 
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The generalized mission analysis approach was used to evaluate a 
broad range of planetary and small-body missions. For most of these 
missions, the existing 30-cm thruster is a good choice. However, payload 
performance for Encke's comet. Earth observatory, asteroid rendezvous, 
asteroid sample return, and out-of-the-ecliptic missions would be 
improved by higher thruster beam current. Beam currents of 4 A or more 
would significantly reduce the number of modules, reduce the propulsion 
system specific mass, and Increase payload. This characteristic is 
similar to the module power depr-ndence discussed previously. The primary 
difference between the Earth oialt and planetary missions ‘is specific 
Impulse. Planetary missions tend to demand lower specific impulse 
(e,g., 3000 sec) for high performance. Lighter and cheaper power sources 
would tend to remove this constraint and add some flexibility to increas- 
ing module power. 

Several thruster design options were considered for thrusters in 
the 50-cm range. In addition to a simple scaling of the conventional 
circular cross section, an oval-shaped thruster was suggested. The oval 
cross section may Increase the scaling range without significantly 
increasing the total technology effort. 
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SECTION 9 

DEFINITIONS AND NOMENCLATUEIE 


A. abbreviations 


B. 


SEP 

EMT 

LEO 

GEO 

ORV 

OOP 

LS 

MS 

OR 

00 

ORS 

ORP 

ORO 

PS 

AU 

PPU 


propulsion system performance 

I 


solar electric propulsion 
engineering model thruster 
low Earth orbit 
geosynchronous Earth orbit 
orbit raising vehicle 
on-orblt propulsion system 
large space system era 
matured shuttle era 
orbit raising 
on-orblt 

orbit raising self-~powered 
orbit raising payload powered 
orbit raising one way trip 
propulsion system 
astronautical unit 
power processing unit 


sp 

T 

P 

ps 
P 

P 


PP 


mod 


P 

I 

E 

m. 


tot 


efficiency) ® ’ utlll,aclon 

thrust 

propulsion system Input power (after degradation) 

power processor Input power 

thruster Input power 

total propulsion system )et power 

thruster beam power (Jet power) 

total Installed power (before degradation) 
beam current 

neutral flowrate equivalent current 
beam voltage 
Ion mass 
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8 o 

- 

gravitational constant 

e 

- 

electronic charge 


- 

Ion beam jet energy 


a* 

Ion velocity 

V 

■ 

average beam velocity 

^ps 

- 

propulsion system total efficiency 

Ht 

- 

thruster total efficiency 

npp 

m 

power processor efficiency 

nc 

M 

cabling efficiency 

^pd 

- 

power source degradation factor 

Hu 

» 

thruster measured propellant utilization efficiency 


- 

thruster electrical efficiency 


a: 

mass fraction of singly charged ions 

^2 

■ 

mass fraction of doubly charged ions 

Y 

= 

thrust loss of efficiency correction parameter 

3 

= 

parameter defined by Eq. 3-* 32 

w 

es 

atomic mass unit 

w 

m 

atomic weight 

\ 

m 

atomic weight ratio 

AH 

s 

ampere-hours 

Aq 

s 

vehicle initial acceleration 

AV 

tt 

vehicle velocity change associated with trajectory 

C 3 

es 

launch energy parameter 


= 

mission propulsion time 


as 

thruster life 

U) 

cs 

thruster life parameter, Eq. 4-20 

D 

- 

thruster beam diameter 

lOTRIP 


orbit raising mission round trip propulsion time 

^eq 

= 

equivalent 30-cm thruster beam current 

^mod 

B 

beam current associated with a single thruster module 

P 

as 

perveance for an electrode design 


s 

perveance for a single aperture 

y.t 


total accelerating voltage 

Va 

s 

accelerator voltage 
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Na 

A 

ds 

da 

S 

J(r) 

Jo 

R 

^tot 

C(Ag) 

C(a) 


a 


e 

A 


net 

net 




<^1 


max 


4 



^dis 

^sc 

^back 

^anode 

^lost 


qr 

K^. K2 
% 


" number of apertures 

= ion optica parameter defined by Eq, 5-5 
“ Ion optics dimension 
■ screen grid to accelerator grid spacing 
» screen grid aperture diameter 
« accelerator grid aperture diameter 
“ aperture center-to<-center spacing 

- screen grid thickness 

» accelerator grid thickness 

beam current density as a function of radius r 
» beam current density on centerllna 

- "R” ratio defined by Eq. 5-11 

- ion optics aperture "total compensation” 

- compensation requxred to correct for grid curvature 
» compensation required to correct for‘ aperture offset 

electrode curvature semi-angle at edge 

■ beam vectoring angle 

- aperture offset distance 

** radius parameter used In defining beam current ' 
profile 

beam maximum vectoring angle allowable without 
direct interception 

= grid dish depth 
“ electron temperature 

■ electron density 

■ thermal model parameter, Eq. 5-23 
“ thermal model parameter, Eq. 5-23 
» thermal model parameter, Eq. 5-23 

■ discharge power 

“ screen grid heating power by plasma 
= back plate heating power by plasm 
a anode heating power by plasma 
" power radiated through apertures 
■= total power radiated from thruster 
- constants, Eq. 5-28 
“ constant, Eq. 5-29 
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^4 

Ab 

^max 

q 

S+. S+4. 

C 

K 

N 

T 


constant, Eq, S-30 
beam area 

maximum allowable temperature 

Ion current density 
sputtering coefficients 
constant used In Eq. 5-34 
constant used In Eq, 5-35 

allowable material loss in atoms per unit area 
lifetime used In Eq. 5-36 


C. MASS DEFINITIONS 


Mo 

Mf 

«P 

Mpjt 

Mps 

MF 

MFR 

Map 

Moop 

Mrt» Mot 

Morp* Mop 
Mrpt» Mopt 
“rpt 

Mrppt Mopp 

«pp 

Mrw» Mow 

Mrss» Mqss 
M rsra» Mqqjjj 

“am 

“ps 

Rt 


- vehicle initial mass 

- vehicle final mass 

" propellant mass (general analysis) 

■ propellant mass flowrate 
“ payload mass 

“ propulsion system mass 
“ payload mass fraction, Eq. 3-5 
" required payload mass fraction 
“ satellite net mass 
“ on-orbit vehicle mass 

■ thruster unit mass 
“ propellant mass 

“ propellant tankage mass 
= propellant tankage fraction of propellant 
“ power processor mass 
“ power processor specific mass 

■ power source mass 
=• subsystems mass 

“ structure and mechanisms mass 

structure and mechanisms mass fraction 
“ propulsion system speclflp mass 
“ thruster redundancy 
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N 


ot 


N 


H 


op 


mod 


N 


ow 


used 


- power processor redundancy 

number of thrusters (with redundancy) 
number of power processors (with redundancy) 

- number of operating modules (thrusters) 
number of thrusters required by wearout 

mass of thrusters worn out during one round trip 


D- COST DEFINITIONS 

^or * total orbit raising cost 

^oef " ORV utilization cost factor 

^1 " ORV hardware cost 


Crt» 

C 

ot 

- 

cost of thrusters 

^rt 


- 

thruster unit cost 

C 

orp 

, c 

op 

It 

propellant cost 

c * 

orp 


m 

propellant cost per unit mass 

%t’ 

’ ^opt 

■ 

propellant tankage cost 

^rpt 


m 

tankage cost per unit mass 

Ot a. 

Pt 


■ 

tankage factor 

C - 
rpp» 

C 

opp 

a 

cost of power processors 

C 

rpp, 

C 

opp 

ar 

power processor unit cost 

^opf 


m 

flight operations fixed cost 

*^opt 


■ 

flight operations time dependent 



m 

launch cost to LEO 



■ 

launch cost per unit mass 

r 


m 

Interest rate 

c , 

rsm’ 

C 

osm 

m 

structure and mechanisms cost 

C a 

rs3* 

C 

OSS 

at 

subsystems cost 

C , c 
rw’ 

ow 

m 

power source cost 


rw “ power cost per unit power 

^ori’ ^ol " integration and testing 

<^01 ‘ Integration cost per unit of ORV or OOP dry mess 

■ (Dm« test, end evaluation 


c , 

ord 


DDT&E cost coefficient 
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Nr 

^res 

^rsm 

^oop 

®to 

^to 

^opa 

Cipa 


learning cnrvn parmeter, Eq. 4.30 

- numbar of OEVa over which DDI4E coat, .r 

■ “'■•■eystena coar coefficient ““«lxed 

“ eed mechanians coat coefficient 

coat of on-orhit propnlaion 

- cranaportation coat, earth to final orbit 

naportatlon coot per unit uaaa 

- ““-Mbit operatlona coat 
operations cost per unit time 
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APPENDIX A 
DOCUMENT LIST 

The following document list was provided by NASA LeRC 
as an initial guide in mission selection and analysis. 
Although most of these reports were consulted in the course 
of the study j not all were used directly. Reports referenced 
directly in hhe text were presented in Section 8. 

1. Future space Transportation Systems Analysis Study 

Transportation Requirements 
Interum Report 11 December 1974 Contract NAS9~14323 

T^hnology 

?.* J^squitements and Considerations in Selecting Space 

Tug Propulsipn Systems •«.«ccing apace 

C, J, Cohan AAS Paper no. 75-lfO 

Satellite Power 

w, B. Lenoir and R. E. Currie, Jr.. February 1975 

6. Mission Roles for the Solar Electric Propulsion 
age (SEPS) with the Space Transportation System 
Northrop Services. Inc. Final Review Presentation. 

January 1975 NAS8-30742 ation. 


neacrlpi"n““‘'''‘' Propulsion Thrust Subeyatem 
JPL TM 701-209, 1 February 1975 


8. Concept Definition and System Analysis 
Electric Propulsion Stage ^ 

Vols. 1-5, Boeing, NAS8-30921. January 1975 


tor So,lar 


Lcke^rwo”' Eleetrio Slow Plyby of 

JPL TM 760-90 Rev. A, 25 January 1974 


Comet 
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10. Interplanetary Spacecraft Design Using Solar 
Electric Propulsion - A Circular 1.0 AU Out-of-the- 
Ecliptic Mission 

J. H. Duxbury (JPL) NAS 7-100 

III EEPC Paper no. 7A - October 1974 

11. U. S, Solar Electric Propulsion Planetary Mission 
Candidates: Out-o£-the-Ecllptic, Small Bodies, and 
Orblters of Mercury and Saturn 

K. L. Atkins (JPL) NAS7-100 

III EEPC Paper no, 74-243 October 1974 

12. Mission Applications of Electric Propulsion 
K. L. Atkins, AIAA Paper no. 74-1085 October 1974 

13. Comet ary Exploration: A Case for Encke 

K. L, Atkins and J. W, Moore, AIAA Paper no. 73-596 
10 July 1973 

14. Solar Electric Spacecraft for the Encke Slow Flyby 
Mission 

J. H. Duxbury, AIAA Paper 73-1126, November 1973 

15. A Study of the Compatibility of Science Instruments 
with the Solar Electric Propulsion Space Vehicle 

JPL TM 33-641 15 October 1973 

16. Feasibility Study for a Solar Electric Propulsion 
Stage and Integrated SEP Spacecraft 

NAS8-27360, 18 January 1973 and 27 March 1972 

17. Thermoelectric Outer Planets Spacecraft (TOPS) 
Advanced Systems Technology Project 

JPL TM 33-589 1 April 1973 

18. Solar Electric Propulsion System Integration 
Technology (SEPSIT) 

JPL TM 33-583, 15 November 1972, Vol. 1, 2, 3 

19. Solar Electric Propulslon/Instrument Subsystems 
Interaction Study 

TRW NAS2-6940 Final 30 March 1973 and Mid-Term 
19 September 1972 

20. Solar Electric Multimission Spacecraft (SEMMS) 

Phase A 

JPL 617-2, 30 July 1971 and 617-4 3 March 1972 
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Voi* IaI*^ 1B°^11 ““^*^i‘Mlsalon s 

» » 11 TRW, jpi. Contract 9523941? 

Solar Electric p. ’ IS March 1970 

®*="4y I'^opulBlon Asteroid Belt M^ ^ 

JPL Vol. 1 5070.21 1 „ »®lt Mission 

5070-21^3 S070-21.2, vol ttt 

JPL Contract 952566 19 ja« 

19 January 1970 

Spacecraft^'^^^*^®*^ ^lyby Mission Using a Sol 

JPI ASO 760-18 1 M w Electric 

'ww AO, 1 March 1968 

31 October 1973 * Allan and P. r. odom, 

^6 Aug., at 1975®”** aas Paper no. 75 - 153 , 

Sp«e mght’‘^“““'’" *” "Comet.ry Explorer." Sodd.rd 


Spaceeraft 
Pln.l Report‘^aL%*!L!>'!‘P" 


.aae 

PUaeta"" ePa«d lour mIso I ona to rh 

Ap Fiandro, Pina, » o the Outer 

april t 970 **»”« EASd Gtoot uo. SGR 45-OO3-0 
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si^ri-'Appucitirn.: 

Report no. M-21, March 1^7C 

33. sul.r Electric Frcpclslon for Jupiter ana S.turn 

Orblter Missions Julvl970 

1.1. TeR.Ie Report no. M-24, July • 

34 B.ll.y'. . Oo»et .nythrough ana Renaervou. Miasiona 

Via Solar Electric Propulsion ^ 

1.1. T.R.l. Report no. T-28, May i.971 

35. Mars surface Sample KeturnMiaaiona via Solar 
Electric Propulsion^ ^ ^ August 1971 

I.l.T.R.l. Report no. Tr/»» *«8'*«‘* 

. s-K*. «i^nl:ua and Prospects of the HAS A 

and^^oputllon Research and Technology 

inryorklng Group on^Space Rower and Propulsion. 

Vol. 1, 2, 3* 30 May 1975 

37. OAST Workshop Report, 1975 

38. HRL Proposal to RFP 3-74063A 

39. Feaslhillty Study Electric .>ropulslon 

Stage, Rockwell, . NAS-8-27360 

40. S.P.S. Study, heRC, Fall 1975 

41. Earth Observatory Satellite (E.O.S.) , TRW, 
iJAS-5-20519 


42. E.OeSe. GE, NAS-5-20518 

43. E.O.Se, Grumman, NAS-5-20520 

44. Status of N.E.P.. Stearns, AAS-75-164 


45. Thermionic S/L Design Study, 
952381 

46 . Mission Concepts for S.E.P. 

Presentation, 1-18-74 


120 kW N.E.P.S. JPl 

j . 

Test Flight, Melsaenger 


ArriSMDix B 

MISSION STUDY DETAILED RESULTS 

In Section 3 ell mission work was condensed into a 
specific format, called the "generallaed analysis approach.” 
However, several analyses were conducted in more detail with 
trajectory runs. The missions analysed included (1) Mercury 
orblter, (2) Comet Encke, (3) Comet Halley, and (4) out- 
of-the-ecliptic. The basic results for the last of these 
four missions were adequately discussed in Section 3. Ad- 
ditional detailed results for the other three missions are 
discussed here. 

A. MERCURY ORBITER (MO) MISSION 

The MO mission was analysed to clarify certain 
Characteristics that bear on thruster requirements. Speclfl- 
cally, efforts were directed toward evaluating trajectories 
that would Sllmlnate potentially severe thruster thermal 
requirements. By constraining the nominal thrust vector to 
be directed 90° to the sun line, thruster thermal loading 
should be controllable well within present thruster llnlte. 
Since a mission "penalty" might be expected from such e 
constraint, this analysis was Initiated to estimate the order- 

magnitude of the effect of the constraint on payload and/ot 
flight time. 

Two series ,of trajectory calculations were made to bound 
the required acceleration levels with thrust attitude con- 
Strained to be 90° from the sun line; 

* Low Acceleration. This case is similar 

^ launch eneigy 

21^7 kwf I* (maximum of 

21.7 kW), injected mass of 3827 kg, but with a 

J®6®J'jJSf| *:hru8t angle and acceleration level of 


331 


• High Power, High Acceleration . ?or this case, low 
thrust propulsion is traded for launch energy* The 
launch mass was raised to 5720 kg for a C 3 of 
4 km2/sec^, power level was 28.4 kW C38jj^3 kW 
maximuitt) , and acceleration of 2.07 x 10“^ g. 

Figure B-l presents the approach velocities for these two 

cases. Although the low acceleration runs are connected by 

a dotted line, they are actually discrete solutions which do 

not exist at intermediate points. Each solution represents 

a discrete number of orbit revolutions* The low acceleration 

level, with constrained thrust attitude, can give similar 

flight times to the unconstrained JPL case Hut at a cost of 

higher approach velocities. Allowing flight times to Increase 

(utilizing coast periods) can greatly reduce the approach 

velocity. Approach velocity is Important because it strongly 

influences retro requirements. 

Figure B-2 depicts both the approach mass to Mercury 
and the net orbited mass. A 500 km by 6 hr orbit was plotted 
for comparison of trajectories. All cases use the same high 
thrust retro propulsion assumptions after ejecting a propulslpn 
system rated at 55.3 kg/kW. (This is the value used by. JFL) . 

The results show that although the higher acceleration and 
power can orbit about 1250 kg, it is possible to orbit 1190 
kg for a much lower power if flight time is increased. The 
relatively small difference probably does not justify the 
added propulsion system expense, unless flight time must be 
constrained to low values. Also, comparing Figures B-J. and 
6-2 shows that reducing the approach velocity to below about 
1 km/sec reduces net orbited mass because the total velocity 
starts to be more effectively delivered with an impulsive 
high thrust retro. 

The total number of thrusters required is shown in Fig- 
ure B-3. The purpose of the figure is only to determine whether 
the 30-cm thruster is adequate for the mission. It is not 
intended to be used to select a mission profile or nominal 
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power level. The range of Injected masses adequate for a 
Mercury orblter mission was found to be 3000 to 6000 kg from 
literature surveys and other HAG studies. The acceleration 
levels appeared to be bounded between 1.6 and 2 . 1 x 10 g. 

HAG EPSTOP computer runs were also made to determine the 
required number of 30-cm thrusters and sparf^.',? using the 
reliability computations with acceptable thruster reliability 
greater than 0.995. The results, presented in Figure B-3, 
show that a maximum of 16 thrusters would be required for the 
seemingly least desirable missions with high power levels 
(though these missions may still be acceptable). The lower 
power profiles would probably require 8 to 10 30-cm thrusters, 
depending on reliability requirements. For these calculations, 
only thruster reliability was considered since the intent was 

to define thruster requirements. 

From this work, we concluded that the existing 30— cm 
thruster can adequately support a range of seemingly acceptable 
MO missions. The 90® thrust attitude constraint significantly 
reduces the thruster thermal level that would be experienced 
in the unconstrained case. The Impact of this constraint is 
mainly reflected in slightly increased flight times. 

The Mercury orbiter mission was also used to Illustrate 
the sensitivity of mission performance to propulsion system 
specific mass. The variation of propulsion system mass with 
thruster lifetime is shown in Figure B-4, which uses detailed 
trajectory computations from EPSTOP. The propulsion system 
mass (see Section 3.B) is plotted versus thruster lifetime. 

The Mercury orbiter mission is a fairly long duration mission 
and typically would require 32,000 A-hr wearout lifetime for 
a 2 A thruster. If the actual lifetime of the 2 A thruster 
is greater than that value, no savings in propulsion system 
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maas over that shown In Figure B-4 would occur. Or the other 
hand, if the lifetime is less than 32,000 A-hr, spare thrusters 
must be carried to replace those that wear out. Every point 
on any of these curves in Figure B-A has similar reliability 
(at least 0.995) and flies basically the same trajectory. 

Assuming 20,000 A-hr as a representative thruster life- 
time. Figure B-A indicates a small mass savings (A thrusters 
or 68 kg) obtained by increasing the thruster lifetime to 
32,000 A-hr. Alternatively, if the beam current is increased 
to A A, with a wearout lifetime of 20,000 A-hr. the propul- 
sion system mass would decrease by almost 200 kg. The rela- 
tive mass savings from increasing beam current from A A to 
12 A is less than 80 kg. which indicates that most of the 

benelft could be obtained from raialng beam current from 
2 A to A A". 

11 „. thrusters Is ludlceted by the deshed 

lines in Figure B-4. these line, ere not streight verticel 
lines .nneuse o£ the fevorehle tradeoff in reliability at 

period 7,:r:T '"'r 

. hus, a large decrease in propulsion system mess 
(end a correspondingly large Increase in net payload) is oh- 

of eZe:: h“‘“r — 

ere. i ”«i"talned relatively constant. In- 

cr s ng the lifetime of the thruster helps somewhat, hut 
results in a far less dremetlc mass improvement even for 
relatively large changes in lifetime. 

B. COMET ENCKE SLOW FLYBY MISSION 

bctweermlLi!n ^^^terest was focused on the relationship 

l^tlcs. Comet EnckVwaTrhV ^^^-ttling character- 

assessing this sensitivity.**^^” ‘=ypical mission for 
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An Knckn alow flyby (4.5 km/aec) trajectory was almulatedj 
the reaults are shown in Table B-l. The reference trajectory 
was that obtained with the EMT power varied at constant bea« 
voltage. For all options except the last two (which Incorpo- 
rate a two-fold Increase In nodule power), the resulting per- 
fornance variation la less than 10 kg. These results suggest 
that a small performance gain and a reduction In power may be 
obtained by operating near the perveance limit. It Is not 
surprising that performance is Improved by increasing efficien- 
cy or decreasing thruster mass. 

On the other hand, the large performance increase caused 
by increasing maximum beam current Is Important, and this la 
the primary result of this portion of the study. The large per- 
formance Increase is not due to Improved thruster efficiency 
but rather to propulsion system mass reductions due to using 
this operating condition. This option could also be compared 
to Increasing specific impulse, which also yields a substantial 
gain In net spacecraft mass. The last cases in Table B-l were 
chosun to yield the jama module power (5.2 kW at 4100 sec and 
2 A) as for the 4 A. 3000 sec case. Both cases give a large 

payload Improvement, but the higher 1 requires a 35% increase 
in power. 


C. COMET HALLEY MISSION 

This effort was performed at a time when the Comet Halley 

mission was being actively considered by NASA. Since it Is 
quit. Ut. to .tort the ,, „„ 

for the record. . 

1 . Comet Science 


The next appearance of comet Halley has an estimated 
perihelion date of February 5-9, 1986{ the uncertainty is 
due to the postulated presence of non-gravltatlonal decelerating 


Table B~l. 


Operating Conditions on 

Bncke Slow Flyby Trajectory 


Condition 

Total Thruster 
Power at lAU»kW 

Relative 

Net Mass, kg 

Constant I ■ 3000 

flp 

19.8 

0.0 

Constant I m 2850 
sp 

19.0 

+.2 

Constant I - 2700 
sp 

18.2 

-.3 

"Farveance Limit Throttling" 

18.2 

+5.5 

Variable 1 

sp 

18.8 

-2.5 

Constant Utilization Efficiency 
and Beam Divergence 

19.4 

+8^8 

Total Efficiency Increased 1% 

19.6 

+3.4 

Thruster Mass Decreased 1 kg 

19.8 

+10.0 

Maximum Beam Current ■ 4 A 

18.4 

+267.2 

Constant - 4100 sec 

®P 

26.8 

+64. 6 

- J 
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forces (mass ejection under the Influence of solar radiation 
heating). Earth-based sightings will recover Halley's comet 
when the nuclear brightness is about 20th magnitude, which 
may occur in 1983 and could refine the estimate of when the 
perihelion will occur. Unfortunately, the along track distance 
will be the least sensitive to optical measurements, so con- 
siderable ephemerls uncertainty will probably exist. 

The Whipple icy congolmerate model of cometary behavior 
is currently most widely accepted; it views comets as "small, 
fragile, low-density structures of frozen gases, well mixed 
with dust and larger debris. Approaching the sun, the gases 
begin to vaporize from this nucleus, forming an extended 
atmosphere of dust and gas (seen as the coma) which streams 
back under the influence of solar pressure and the solar wind, 
to form tails." The burst of activity noticed at about 3 AU 
from the sun may be caused by water vaporization, which is 
thought to be the major volatile component of most comets. 

Imaging the nucleus has a high priority because a nucleus 
has never been seen, only postulated. Mass spectroscopy and 
dust analysis also have high priority in studies of the origin 
of comets and of the composition of the nucleus. Measurements 
of the existence, location, and physical properties of the 
solar wind and bowshock interactions (with potential contact 
surfaces) is of great Interest to plasma physicists. 

The desired mission scenario, a spatial survey, would 
probably satisfy most science objectives. The spacecraft 
would enter the bowshock, traverse the coma, and then image 
the nucleus. Then it would fly slowly through the tall in a 
manner that would maximize the data collection time and 
distance from the nucleus. 

Tradeoffs among science return, instrument development 
requirements and costs, and approach velocity (trip time 
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technology) 
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Figure B“6» Ion drive flight path to Comet Halley. 
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AV. PROPULSION SYSTEM VELOCITY ADDITION, KM/SEC 
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Ag = INITIAL ACCELERATION x IQ-B q 
Figure B-8. Typical propulsion system requirements. 





Figure B-t. Comet Halley mlaa^ftn a 

launch date, performance analysis, 1982 
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Figure B-10. Power requirements, 1982 launch date. 
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Figure B-11. Performance analysis, 1983 launch date. ‘ 
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APPROACH VELOCITY TO HALLEY. KM/SEC 
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Figure B-12. Power requirements, 1983 launch date. 
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are presented. The first is the present engineering model 
thruster (EMT) with present technology power conditioning 
and solar panels. The EMT operates at 2 A and 3000 sec I 

sp 

When spares are added* the propulsion system specific mass 

a becomes 45.5 kg/kW. 
p s 

The next two systems have modified EMTs . The first 

assumes an I of 5000 sec; the second assumes a beam current 
8 p 

of 6 A. Power conditioners are modeled according to present 
technology estimates of masses. In the second case* Increased 
wearout rate of higher beam current thrusters necessifates 
thruster replacement due to wearout. As one thruster wears 
out* a replacement thruster is connected to the power condl?- 
tioner in a manner similar to that occurlng in a thruster 
failure. Specific masses of the two systems are 33 and 29 
kg/kW* respectively. 

The final system modeled contains an advanced concept 
thruster operating at high beam current (6 A) and high 
(5000 sec). This thruster is assumed lighter and more ef- 
ficient than the EMT. Corresponding breakthroughs occur in 
power conditioning and solar panel masses to yield an 
extremely light system having a of only 14 kg/kW. This 

PS n 

is still somewhat higher than the ion drive system with 
a value of only 10 kg/kW (a 4 kg/kW difference). 

Figures B-9 and B-10. present the results of the March 
1982 launch opportunity. The net spacecraft, mass is given 
in Figure B-9 as a function of desired approach velocity and 
technology level and compared to an estimated mass of 450 kg 
required for this type of mission. None of the technology 
levels is capable of rendezvous with the required mass* but 
the advanced technology comes very close (« 1 km/sec) . In 
fact* the extra payload delivered by the ion drive system 
with its 4 kg/kW difference at about a IOC kW power level 
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could achieve the rendezvous* The three present power 
technology levels do show some approach velocity improve-* 
menta with thruster changes at 450 kg (18 km/aec to 12 km/aec), 
but, as seen in Figure B-10, higher power is required. 

Figure B-10 depicts the same trades as Figure B-9 in terms 
of power required for a 450 kg net spacecraft. Whenever a 
technology's performance is greater than 450 kg (in 
Figure B-9), the Shuttle payload is scaled to decrease the 
power requirements (i.e , , Initial mass), but acceleration, AV, 
and propulsion system characteristics are kept constant. 

The penalty for the 6 km/sec decrease in approach velocity 
is now clear — over 30 kW of Increased power is required (i.e*, 
from 50 to 84 kW). At a potential $500,000/kW of solar panel, 
the cost for the 33X extra science time is very high. Simi- 
larly, the power/cost for reducing the approach velocity 
towards zero is extraordinary for even the advanced technol- 
ogy. Reducing the approach velocity from 4 to 1 km/sec more 
than doubles (from 50 to 110 kW) . These values are con- 
sistent with the ion drive results when efficiency differ- 
ences (0.74 versus 0.76), acceleration levels, and a are 
taken into account. Small technology changes require large 
power jumps because of the asymptotic nature of the curves. 

For each curve in Figure B-10, the last few km/sec changes 
in approach velocity (above the 450 kg limit) causes high 
relative changes in power level. A clear message is that 
the trajectory requirements are barely satisfied, and any 
slight change in technology will be extremely costly. 

The 1983 opportunity is depicted in Figures B-11 and 
B-12, in which trends similar to the 1982 launch are exhib- 
ited, The basic difference is the minimum approach velocity, 
which is now about 10 km/sec for the unconstrained arrival 
date. The three present technology curves (EMT and modified 
EMT) are also closer to each other in approach velocity 
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APPENDIX C 

EARTH ORBIT COST MODEL PROGRAM LISTING 

Computer programs were written to model the cost methodology for 
both the orbit-rising and on-orbit cost models described in Section 4 
of this report. These programs were developed and used with the DEC 
System-lO time-sharing service and were written in FORTRAN IV. Copies 
of the listings for the three principal versions of the orbit-raising 
and the on-orbit programs are presented in Tables C-1 through C-6. 

These listings are for Shuttle era bas-^line cases run using Xe pro- 
pellant and an Initial mass of 25,000 kg. 

The three versions of each coat model produce outputs which 
calculate 

• Cost per net payload ($/kg) 

• Normalized cost breakdown for 
each major subsystem 

• Mass breakdown (kg) for each 
major subsystem 

Each of these programs plots the chosen variable (e.g., cost per net 
payload) as a function of 1^^ with any of the other input variables as 
a parameter. Since most of the equations are common to all of the 
programs, only the orbit-raising version is discussed here. The 
nomenclature used for the various parameters follows that used in the 
computer listings in which all the program variables use only capital 
letters. In the listings, multiple statements on a line are separated 
by semi-colons. 

The computer program uses specific impulse (XISP) as the inde- 
pendent variable, over the range from 3000 to 18,000 sec. For argon, 
the beam voltage VB is given in terms of the specific Impulse XISP as 

VB= [xiSP/199]^ . (C-1) 
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table C- 1. 






c 

c 

c 

c 


HtAi I .M LOST Muptl 

pi M / 

‘^ATA A/s 6 !, so., loS./ •' 

Ilf - 

Tj lilts 


c ifi.V^*.'' Jl-ASL^r,! 

lilts 

PMoSi:('i’;'i'=-5*il«E''"'‘'’'*'‘-‘''''“tFx = .Cuef(i, 

IFI ifui?!/p’c-'*^^ POLtFhs-* 1 yi^*?''f.i-OA + ,7i 
1 r » HiUlvvt.Fw. i) SOLPFtl 'i ‘l^Fi.Jx+ u-* 

of kx;r quaudi 

5 ,tJ^fptLy/GlSH * ^ 

X.v.f,XNg/l:APiHt.TA) 

»l'i: jf S"""”' ' f "=eff.s.,uff . . 

IhOs iMOuA^iOO,/ 

i'^L’f:‘;?r.*.'!Lnti 


r«c, .,„=.. 

? iK W is;;::::::'-'""'-- 


^i'^Fih=nMt ,xi I •»- ' 
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TABLE C-1 (Continued) 




C <B/IAHD>UHVHI UtOKjo fKl dO-jMU-7a lUTHJ' 


HAUt i:l 


C THRUSTLR CoSI USING dOi hUuE 
C ONLY pay fOrt tHKUSTtA- tlFfc^UStU . , , 

UdErv0r=u01 MF tf AC » .ul . I . ) UofcNO I sf ACT *^UT 
CUSlsMa.FjlAALOGlOC^) f CPCoa.»Cuai*U3tNJl **.60/,6o 
C PPU COSIS AY ASSUMED PPU tFF OF VSX 
PPUP^HsP ( J )/T()TEFF 

uul lsi4,iF.b*HPuP«R**.ia;CRP=UNU*NOP**,68/.6b 
C INCLuDt DOWN PhOPelLAM IN CPKUP AND CUPfiECT PAYLOAD 

CPHuPsCk'^’P + i'IPD'NN) *10. JC I Ai\,K=. 1*100. *(X‘'P+f«iPDl^iV);M0P=VSP-M^ Uv'<fM 
C UStD IhPUSIEh /.EIGHT 


b 

t) 

bl 

bS 

60 

b2 


.tO*lvtRSf’**.> ^ 
MUKP+USfU);LUPF=l .Eb 
• ■ *PTur 


II# rnruciicn ^ ^ 

Ui>FD=/MA3S*NMLD*FAC 1 * 11 , **<01) ? CK&“1= I . t 
CkSSsI .L 4AMAS5** ,y»CUKi=0 U. *(nOPV-MUKH 
C0WU=l .tb/lO.*(MOK\/-MUPP)**.5»CWrt=lUU,"r 

CllPl=l.t5*IOTkIPJCL=300.*lX..lO+ySLD) ^ 

COST isCUkl+CRI ♦CHP*CI ANin + LHSM+CKSO + CPo + COkO 
CUtEF=TuThlP/L0kV+P,7t-0 * IOjPIP 
CuSYsCOtFF*tn&T I ♦LOPFiCUPT+CLlCPROP; AVCCSY =COSY/PSP 
XsXiSP/iOOO.nsAVGColMKt.tl.lO) GO UJ b 
Ys2,b*AL<'UUi(.l*Y>; iF lY.GE.T. J i>U lu b 
CALL PLuUX, y, IPEU) 

IPEhsP 

OlsAflS(oO.-i.O;l.V=AHo(hO.-D) ;D3sABS(10y.-0) ;ua5*‘uS(lbO.-0) 
lFflil,Ll..lO) CALL syHHuLlXrY,,lr'll»0,,-J) 

s ynhul (x.y^aif '* 12 » 0 , » i ) 

6Y^eUL(X,Y,.l,'^3#0.f-lJ 
SYi' AUL(X, Y, r-1 ) 


, * # M 


Call 

CALL 

Call 


IF luP.Ll .. 

If (L5 .lt 

IKD0,Lr..<J4) 

CGNTIrtUL 

CALL PLLmO.»o.»-3) ;C>0 bl Ja=1,10 
YsFlOAT lJX)#Ys2.‘S*ALOGlO(Y) 

CALL SV r-duL LO , , Y » . I # 20 , 0 . , -1) 

DU 65 Ja'=(2u» 1U0« lU 
Y=.FL0AT(jX);Y=2.5*ALUGlyfY) 

Call oYnhuh 0, , t , . 1 ,2 a» 0, »-i ) 

DU oO JxscfOy » 1 000 r 1 00 , ^ 

Yaf LUAT IJX); Ys2,b*ALUL.lO(Y) 

CALL SY^oaL(0.,Yr.l»2‘4»0.,-i) 

Du b2 JXsl » b 
XsFLOAKJXl 

Call aYrv.riUL(XrO, » . 1 f 15»0.»-U 
Call PL0r(().,O.r-5)JCALL POFF • 
31{)P?tND 
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TABLE C- 2 . 


ri 


PAGk 1 


C 


C 

C 


GRhn f.Hf ^•-uUH^5iOSP'rL^^W''»•■•'^^^♦ 

tslAw; 5 V 5 X'u;r;;^v'a^^s^ 

Kir;' wi'i^PP /P.hBf . 


LV 1 fU NS I UN A I b ) . K o # » ) 'i 

nitA VCoEF Va.bBr . 

KliHf -'lii 


Upon PA ib ‘ ^ r ^ l,,> OLL 

!:iRU'iouo:?Sui=.^#«^p=-‘^ 

uu.o jrAat=i,»i„ 


no 0 ufASt=l»i 


uo b 

KjiiudVn 

E?‘i i-^o’% vr.Pi UPVCK I U 1) ; VCUEP APVLOtP u J 

?t^tUisb6;.N=;sniMfc 

F i" 1 ’» i = *1 i 1 1 -3 Si 5 i‘ C u V = A L n b 1 U « u t - 1 S * f L u X ) 

lCUKVt = S . r . ,ot j*y4H lixi 


lCUKVt = :>, .., ,v \ni ► f h S-. l‘n<J ‘•XXKLUXi .7i 
IF ( ir.uPvL .tU. » ) -., I /:''iAXxFLU.<4 

ItimiU'E-.ll;:!) ltlkr'F=-;l?v.»vFL'jx..‘.<'s 


l;’'('vr.runcKux.'/-) x'‘ 

t,U IN y\ 


|ir •■ ' , * oUH « \ i>' '.’’’F * -F • ' 

j;^^i'.’.&i(i'H.ivvi;utF«.iv,../ioao.)*-i. 

CUM n put ^ 

x.-AbS=.07b*y*M.5'3 

UP 1 fl Slip L Wu 1 SI , ^ 


bJ 


UP 7 ^ Slip L V /U I SI 

x^FiXr.u/FxPua lA) 

XNiK=Xf (J-Xl*iF 

. . ^ « j r 4 4 A 


jeNiKsXi^ , i ii 

HiFTs)iF'H*i5si'‘*«iSF»'.’bt6FStr. 

1 0 (P (0 ) ) . l . lbl 

fi:Et'’5lS5bM5k‘F2?.).)MC..FF=f.FF.SI 
•-■'MsKJ£T/TuUFF 
u=p.it r/PMU,' 


;fF A8ULEFF*.^b 

KilV/Ki't, .,, 

i^sii'ur75Su5bSo(iF.u)t<».-i 

UKk’oM-^Jl^xlNE FOK VAi<U.=U IHPUSttK LlFl, 

79 ji!Fi>=!;,';'^rES./l ,:f.LT=FAUOKilFa 9 CIOA.LI.I.) FACn',.= l. 


-9 TUTPlPf /LlJf ALTsFAI.IUKJ 

*V.iu ALTu«* I.i I uPPUaP I U i * U . FPOH) 


e 1 


11 


u 


I'h (au)wvIgI .X^'H.O -F xP C-f'El a J ) 

IKUFl.lT . /) Ml >U 3/ 

OOUNsXT IKL tuO 
57 U>1 Pll'= » li'ik tA 1 1 -l 
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TABLE C-2 (Continued) 




PftC’l 1:1 


,bfc/,brt 


Flip I i 1 ' n * I U T 1 H Ll - 5u 0 . * ^ -tut Uot U ) 

CosIlsCoPl I ANc.tUHi>'];tf(;PSbfCHV**L''AO 

>Vp=rr;l KPMt./Lv)Sl;^¥2^^1U 
xi^i*>p/joyy» V— 


b 

o 



AuSllbO.-O) 


I V r. 

CUNT t^Ut 
cui'i I i i'fOt 


1.0 b^ TXswn 

Cttut^^b vipUL ^ ^ ^ • 1 ( } J-' >r,J H 

SlOFMOi) 
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TABLE C-2 (Continued) 


C <BrtAKl)>0KVi^^S1 .f OhM «'0-JAl\i-7a PAGl Itl 


C IMRUSTtf^ COST USiOu ttO* KULF 

C ONLY pay hOh THraioltK l.lFt UStt) il 

uSfcf»OI = Nni;lKFAC I«Lr,l .) O6h(i0lsFAcr*NaT I 

CoSl=a^,Ei*ALUG|0(OJ ;LPI=2.*CUbr*USLf'.ur**.no/.t>b 
C PPU COSTS AT ASSuFitD PPU FFF OF 95)1 

PPUP*Nh=H(J)/T0rEFh f 

UNn=A,3Eb*PPUPwF**. 12;Chh=UF:iT AN0P**.6d/.6b I 

C iNCLUOe DU.-*N PNOPFLLaM iN CPRoP AwD CUKRECI PAYLOaO » 

CPROPs ( xMP + mPlv.'N ) *10. ;CTANK = . 1*100,; ''■•SPsMSP-M'Oim'J 
C USED ThPuSTErt ♦\F.1 GhT 

USfc U=XMA5S*iMMol)*FAC I * C 1 , + KO I ) ; CKS '.si .t.«*MhSW**.i f 

CKSS=l ,ta*MWSS**.9;COKliAO.*(nOKV-M(.)kP + UScO; ;C0PF = 1 .fcS I 

CuPi/sl .t5/ 1 0.* (MOhV-MuRP) **,3;CRA5 s100.*PToT 
CuPIsl .E3*ruTRlP;CL = 500.*(Xi-i0 + USFU) 

CuSl l=CURl+CRP*CTjvHK + LHSKFCRS6*CPvj + C0hD T 

(,uFFF=TuTkIP/LOhV+?. 7t-4* lOTKiP | 

COSr=(:nEFF*COSTi+CR r+CriPF+COPT+CL+CPRuP? Av(»CST=COST/MSP ^ 

XAF=lU,/.Xt'tO 

x=x iSP/iOifO. 

IFCJ.EO.U Ysrii; f*XXF ; iFtJ.FO.c) YsMKPP*XXF ’ ? 

TF(J.Ln,3) VeXX:;F*FU,<F; CF ( J.tO.a) Y=XXF*MRPT I 

IFCJ.fc'J.S; YaXKF »t.RSM; IKj.tO.oJ r = XXK*!/iRSS 
IF (J.to.7; .Y-ex?.;- M vVm; IF (J.Fu.a) YsXXF*MPIH,N 

ir(j.tu.9; ’YsiviSp**XXF » 

IFIY.LT.Oj GO TO S 
CALL PLUT(X, Y, JPFiO 
iPtNaP 

U ( J.OL . I ) to Tu S 


6 CoN I I (mUL 
V DO 52 JXsirb 
XsFLOATCJx) 

5<; CAIL SYMBUL(X,0,, .1, li,o.,-l) 
CALL PLuTl0.,0,#-i) ;C mLL HOFF 
STOP’LL NO 
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Ol = rtBSC30.-0);OdsABS(bO.-D);Oi = A83(lOu.-DJ ;tia=ABS(l50,-rj) 
IF {L1.LT,,10) CALL S YMHUL. ( X , V , , I , M 1 , 1) . » “U 
IF ( 02. L r . . 24 ) call . S Yr.4iOL ( X , Y , . I , , 0 . » - 1 ) 

IF (05. LT.. 24) CALL SYMBOL C X , Y 1 , Vi5, 0 1 J 
IF(UJ.LT..24) call SYMKCL(X,Y,.l,M4,o.,-t ) 
lPt.ua2 
LuNl IhUt' 


I 

1 
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\} 

Nf 


M 

w 

0 

il 


lAOLic 


C 

C 


c fRj eo.ju,-7.i ujsrm 

^ Ii **h.V’ (**?. cREftAOO«^\ 

'fP!; 

iNiiiAL cSfiiJiriols 

CALL POr F 

call^p6,V UZl 

E wf ^ |V|oJ. 5 S ! b 1 7 liiE ^ ^ ^ ^ ^ ^ ^ ^ j 

P = H(J)*1000. 

Ilpiiiii i1lf“iiiH:!Si!!:;i; 

Ih (VH.L I .vCWUXj bU Tu Pfl 

ci’i!‘,uj'*'‘'"‘'''''’^'''^f‘*<''“''‘i'«.)**-i.sj 

XMA.bS=,t)7b*i>A*j ^5^ 

H£ TAaDfcLV/GlS^ 

^:"-f=J<^'^U/ExPlBL TA3 
XmPsXmO-XMF 

StC=TiMt /100.*8.t«y£6 

‘-SiO/StC 

.V.piJ = P.JET/P.l;Ou 

It.UilMgij*90o./tU*l;) 

4T2l5!:”Ji5'-°'’'“‘'^''>*‘'-'' 

'*''l iSllEa^*' '■'“’* '"‘WAbLt tHRUSTt'R L]F£ 

IF(«5Ky;Gll*viy'4!l®?,'',*^''’'*‘'«Pf‘*'^UI(P«MRPTtft(lSM»HKS!>Aitl<^ 

w'^ 1 “5^ 1 7'"*^^ * ( 1 xP ( •f F ! A 1 ) 

,/ nO^:^. = XTlM^;GO ri 79 ^ ■ 

3/ mT«iPs 


PA lit I 


2\i 

i'i 






1# 




.95 


1 • ) PAC i LKw j , 
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TABLE C-4. 


<B.AK1»MaSS.^0k;6 h. 1 cO-JAf.;.7o 


3;i4PM 


c 

c 

c 


LrK)T 

A. ru.ciio.) 


pagl j 


DO S^jSt.tCJ.^j CPfcNiS 



CHuStf'* 


• • • •» ^ r V ( 

C X X= 1 /■ , 

i = J ^ Ifii^O 

. gisii-‘.oUj;!;'i^^'<'->'>'<'o«. 

j 3 l.'n';n ‘?f '*' ' ''"^''''‘^•9»‘l Vr/iOu(i. ) «»•! .5 I 
2[' rjy«-.ij7a*0**l ,313 
t'PJA = |jf-|.V/6iSH 
XivF = X(v|.j/t-,xP(f3tTA) 
x,.iH = X|‘lO-x^:F ' 

M=i,t'4 >• ■ x-HpiAblt fH.xUSrtH LIF^ 

f jm£/ fLi/^^^ j 

i tHiS (-xAtj-.ui L. l^.. . 


it* 1C 


wi-iry I r h-ui i i n»; T v c. •« i 

iiiK&i i;Sii3Fi«tf?f;&«K,aKf,,, 


>^’lJ.SH = XviL;-viOijil '-BUI - 1 -Uo .,<nr)U + • -uPX lOP H ■ + j,-, 


364 
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C <H»<AKf'>i“'A.Si3 Ok » fe 


TABLE C-4 (Continued) 
♦■fii «?0-JAM»7e iSlA'kH 


f’AGL IJI 


( I j.t, j 

rnCw - i • * . ,w LUuOsl .f b/ i U. * -JOH) ** i 

r A • V ; i * . 9 ; c T U = 1 0 0 u . * , lOjp JU K J * * . 

I** tC<)Pb*CtlLiJ 


+cjspifLOi>stc ro+c'/I n‘ 

v’-y^i ci5 '' V(.CS f itObT/ ( -lOiiPn Y,<S) 

V ^ i t ^ ^ y. T A y ^ i' V F ' A »■ M • •r M. . ) i» u 


IP 
IP 

IP (D3.I 


ij(^ TU *j 


lO b 


, \=f?*!?,''«Li^liI0(Y) ; if ( Y,Gt,'/IT 
^ I ^ y # I P'P.t^ ) 
iPFj; = CPt.M 

lall Jj vyhOu ( a » y , , 1 , -la, i> * , •!) 


r 


bl 




h'l 


be! 


ih (ua.L I 

CUN I I I'jijt, 

CU'^TMJc. 

»■ Kl II LOb I IC PtARKo 

^=fLUrtT ( JA) ;Ys?.5*ALUbIy( Y) 
Call bYf-.ytiLio. , f , . l,^a^o. ,-i ) 

^A=PLOA r cjx) 

CAI L b r .'iHuL ( X, 0 . , . 1 , 1 0 . » - 1 ) 

o I u r , t '’j I' 


10 
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®'P0OR 




TABLE C-5. 


C <8lf AKU>MAStS1 .KOf<? 1 t-Kf 2(l-jAN-ifO 


f^AGc 1 


OW OF<bIT COST PRLAKDOhN 

HtAL W^'Oi), ,IJUP»NuT,MjOP,i»(ObP,MUrt* It'i 

KtAL l»UoM,Mul i MuSS » MOpi'i » MuS 

Di MtfjS 1 U J A (») J , M ( (5 ) , McS ( J J , Mf y ( 1 0 1 ) 

Data r/<;0w50. , 100 poo, » 500, , 1000./ 

Data a/ 5,, lU.,Pl». rO. rO,/ 

.t I ID.i t I r ft .1 T Kill I .*1 ^ I i r tT K( n <j j t , i 


c 

c 


ON OUlil I COoT MUDtL iMrtlITtN 
MUDIFIFD UtC /7 FuK Xh, MC 
M 1 = 2 P; c* 50 ;Mi=i ;Ma=P 
iNiiJALl/: \i&nlu'\b CO.iSIAl,Tb 
YhS = lu. ;*nif = .t’;Knh>s.i* 

DtLtf = l t iJ‘L = jtb.* YHS 
CALL Pi.UTSCSO) 

XM(? = P5.t i; iCUh V!:.= i 
DO 6 JCASt=l r I ;CPt.lV = ? 
h fJCASt.t'j.P) lPLNsS 
DO b J=l,10 
LAlL PuFF 

NtXJ LlivL SuOiiS DlivN PLoMF^V SU y*jU Ca- 
PtAL/(5, + ) 7 ZZ;a« 1 I E(S,* J Z/ZrFACTuP 
LALL PiIn;PC1)=2u. 


Ht iITlL aaPO at mML lP /7 


NiJVriF.K rnt lUk'i/cS 


FLUX = l.ttb;XXFLDX = ALOblO(l.t:-i5'*FLliX) 

!{:( JCUPvt.tu.n SuLtFFr-. 1 9 1 P * X xF LU/ t . 7 i 

S0LtFF=-,l7/5*XxFLUx+ .90 
SULhFFs-. 129*xXFLJX + .99'i 


IF(lCU»?vF.fU.2) 

IKiCUi5VF..t^J.S) 

IPtN=3 
PKfJljsP( l)» 1 000, 

CXX=17. 

ryvv'=5 

lUJ.Lo.l) 1VVV=1 

i)D S 1 = 1 , 1000, WW 

X1SF=CXX<FL0AT(1-1 ) +10GG. 

DiSP = ,009ci*XIbP 
Vi>=3,PB7*X13P/l 99.»X lSP/199, 

If*iOiJ=K'^iUi>/ Oil 

U (VH.LT .SOUS. ; GO i 0 20 
D=40.*Sij«» ( IMOD/lb. J 
tiU TO 53 

20 D = '5U.*SOH) ( iMun/P.flb* l''B/l()oO. ) *★-! ,5 j 
35 CDMIlMJt 

X-sAbSs . 0 7tt^D* * 1 , 5b 

bLTA = (*FLV/&iSi- 

Xi'lF iX^.d/FAPlbL f A) 

xr'’,PrXivlO-XMF 

sec=TiMt:/ioo,*b.69eo 

PJF IrxMP^lilbPAGISPA.bbto/SEC 

t XPuDxs-.bOi + ALoGiO (P( 1 J ) + J , lb5 

ALPPU=lO.**t.XPONX 
L{-F=.b*4*Vr}/ I VH + PPU, ) 

AbSUMf kPo tFF=COi.Sl Af<T5.9b 

J,pTtFj:=tFF*S0LFFF*.9b;Plur=PJL T/TuTtFF 
• _NKaD = PvIET/Pf-iOD; TEQ = IMOD*oOO,/(D*Dj 
U'A=-1 .27*AL0G1 0( lE 0)+a.u;Ll = 10,*itTtM 
IA^F OUT NfcXr LINE FuP VAPlAdLh IHKUbTtR LIFE 
L r=l .LA 

^ FACTDK=TIM&/(LT/?a, 1 

IF (FAlTuH.LT. 1.) FALTOPsI, 

(JUT = NiVUii*f ACTuR*(i ,tK'uT J 
NOPsMF.Oti* U .+KOP) 

uo.t = .,Ei...,-«v. COST=',e.S.«J 
PPOPWK = P( 1 )/TOTl-,FF 

U'j 1 J = (4 , it b *HPljP»<P* * , 1 2 J C DP= JN i r * <90 P * *,htt/,6rt 

*XMAbb^iJOT ;(';(ipp = ,o0 i * ALr'Pu*H T i.iT ; MuPs X""P : Mof'I = , 1 *.'. uP 
^ • J + 1 •/''<DP-1 . ) ^'^USbsVUbS* 1 1 + , U 

MUSf ! = , 2 * (i''!0 I Ff'iOFF t MQP + I'iGP T + MO 3S J 




TABLE C-5 (Continued) 


C <HrtAKl»MASLSr.^OK^i <»0-jAN-7d 


iPtMOuP.SyilMu) 

^10s^=AMtl-.v^aop ^ 

KSlsHlKiS'-Si-SSi!'"”"'"™* 

ifu'^'sl V-xIfIcuv' M f '='<'F*cI)Sh 

f>o io s 

C-jMl ItgUl 
CU 'j T ii'.Ui 
KH^ii-lA j ( lOtlO.i) 


PAGt l:l 


®®Rj®7 AU PA'tK f 

paEQOK quality 

* 

'L/ < 
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TABLE C-6. 


C <8AAul»MikSI .KUKM F»?X ?0-JAN-7tt i:b6PM 


PA,Gh 1 


Oi>i ORBIT hASS BREAKOQUN 

«t AL NMUO, I hOU r L f » NOP # (MU T , MuOP , -lUBP , MIU' • 
heal MUBMf MuT.MiPP,tMOP,^nHr»MOSS# ^lOP.M, .^ 
DiMENSlON A(5),P(6),PESl?J,?YYll01'J 
DATA P/£»0 .,b 0. r 100,# ?y0*F560.» 1000./ 
^AtA A/l5,,10,f20, #0,f0,/ 

OHOIT Cost itOUfcL .-.WITIF.R BY BILL ,'iM< 


#If.O 

uS 


BILL .HARO Af hPL 1977 


C ON OHOIT Cost itOUfcL f.PlTlF.N BY 
C MODUItD OtC 77 FOk aE.EIC. 

• M1=22;N.2=0;M3=1;'-i = ^ 

INlTTALlZf VAiaUl'.-* l- ARAinfc' Tf. MS 
LifB,; YHS=i.o,;Ruu.B;i<ups,a 
DtL v = j 6, S; T lt<1t = jAS , * Y t«S 
CALL PLuTSfJOj 
XMUs?>3.c 5F lCIJKVfc. = :S 
OU o JCASk = l # i ;CFtf’=2 
Th ( JCASt .t‘U.2) tPci'.=5 
I < I 5 J ” 1 » 6 
LAlL Po*-''F 

RcAOCb;*) Z/Zfvn-. I lt(b# *i zZ/#F«C10R 
CrtLL MOt ;H(1)-2 u. 

fliix=i .(■ lb; xxflox=aliigio( 1 .t-ib*FLUx) 

IF { lCo>' ve,fc u. 1 ) SuLth Fs-. 19l2*XxFL0X + . 7i 
If ( lCuMvF,tM,2) Si.*LtKF=-, l7/b*XXFLUX + .<»u 
iFdCURVF.FQ.i) S.iLETF = -. I29*XXFLUX + .<i95 
IFF. J=i 

P.iWl!u = P( J ) » 1 UOu. 

CXXsW. 

IVVv=b 

IFCJ.t'j.l) IVWVsl 
Ou b 1=1 # 1000# I VVv 

XiSP = CXX*FLOAl (1-n MOOO, 
f;isp=.oo9tj*xisp 
VD = i.2P/*XlBP/lsl9.*XlvSP/l99. 

IFOO=PMuD/Vd 

lF(vP.L1.300b.) GU 10 20 
DiSO.*SURl(iMOl)/lb.) 

GU ro 3i- 
0=iO. *avJRl liHuD/2,eb* I Vb/1 000, ) **-l .b) 

XM'ASa = ,0 7B*D**l .3b 
rtfc f A=OE LV/L;iSM 
X^FrXMO/exPlRtTA) 

X-iP = XciU-XMF 

StC=rit-*L/100.*8.ftaF6 

PJET = XMPxGli.P*GiSP*.5L6/SEC 
exPuNX=-.b03*ALOGlO(P(l ))+l .1S3 
C PAlbk PrU UASb bY HtDuNUA.JCY FACTUR 
ALPPU=(l.*HrjPj * lO.A’JfxPUMX 
. tFF5,b<T*Vo/l\/t> + 22U, ; 

L AbSuMt PPu tFKsCOKSl AivTs.Vb 

. n i T L F F i t Fi* *4LYU E F F * . 9 b ; P f 0 r = P J t T / T U T £ F F 

NMnu-=>PJiiT/Pf<1Qu;iEQ=If^(iP*goO./(0*D) ’ 

lEu)+A.A;Lr = 10.** I EM 

FttC I (JM^ i^/ (L 1 /29 . ) 

1F(FUTTjP.LI. 1,) FAClUPsI. 
iMul ='‘JMnu*-F actor » ( 1 . +KijT j 
NuP=<J/iUl.*1 I . + HOP> 

ppuptvr< = p(l )/UiTtFP 

Uun=«. 5Eb*ppuPr.K A *, i2;cop=u:'in *riup»x-.63/.6a 
i;'ioT=P, »^.yA.S:>*i\:0 i ; ••f)PP = .uoi* AL pPo*piuT ; MUP=XyP;f'.(iPl = , 1 
N'U.,'b= 1 uO , A ( I , + 1 , t -txp I Ij ( ) + 1 y 0 . * t Nijl / N Jp- J , ) ; ■'^'•uSSsi .(JiiS* ( H , 
MObi 1= ,2 * ( >•'] I +I';ijpp + i ,L.P iiP i tMOBS) 

PU‘! = 0. *.00 I API OT ; lOOFsI'.uSy + .'iU i •('.‘‘.OPP + .V'i'P + yuPr 


2u 

33 


c 


AU 


1 ) 
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I 


TABLE C-6 (Continued) 


C <dviAwn>MA.S| .FOk; 1 ffti ?0-JAN-7a irSAF* -! 


i : I 


6 

5u 

77 


IF(!00(jP.r.I.XMo) Gu 10 S 
MuSPsAMO-iviOUP 
CUP1=AMK*(10. + . 1*100. J 

CUSo=l .t*l*^ioSa** ,‘#;lTus1000,*wOJP ’ 

^ 7^0 . * viuop ? t;o*'<= 1 00 . *p 1 01 ; tnpa= i . E 3 * r r ht 

+t.OP I lOPT +CoSi*i+L0SS*Cr0tC01 tCOW + COPa f GOOD 
i5}SP::XM0-MfJoP;AVGGS1=C0al/(,-1OSP*VPsJ 
X = XlbP/30 00.;xxHslO,/AfKi 
IFl.J.FO.l ) YsXXF *:HIT; IK J.tO,?) V = XXK*iviOPP 
Y = X XF *MUP ; 1 F ( J , tO, /4 ^ y = XXF*MOPr 

if V Y = XXF*^OSP 

IKv.bF.o.) bo 10 b; IK y.LT .0, ) GO lU S 
C«LL PLuKXKKPF;.) ^ ^ 

iFtN = f>; IF(J.Gt.P) Go If 5 

l>l=Hflb ( rpcf = ^^3(bO,"^0 ;O3sA0b( IU0,-D) ;i.'A = AfiS( t SO.-O) 

iKol.L ..24J Call oYF.BuL ( X , Y, . I ,r-H , 0 . , -IJ vj 

If 0P.U..^A CALL 3YMOL(X,ylKj.y2;o:|-{) 

IKoS,Li..20) Call oYnBuI ( x , Y , , I , ,^ 15 ! q . » * 1 J 

ipF.SiScPrN' ’ “'"•''’‘"-‘CO.I.Mii.oCil) 

CUFJ I I mi IL 
CuNl Jt'.l'L 

FljRfiAl ( lOtlO, i) 

f Ur *0=' /F lu.2» 'P=*,Fi0.7# ' 1 = ' » 1 7) 

Call PLjri0.rO,,-3);CALL PUFF 
vST<)p;ti.o 
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T1.1. «l«lon u d«lv»d from th. erproaolon 


where 


sp 


v/g 


(UlI 

^1 

L 8 

J 



2eVB 
If 


\ " 0.81 
^2 « 0.067 
Y “ 0.958 

M " 6.63 X 10”^® kg 


(C-2) 


(C-3) 


^'or Hg, Xe» or Kr propellants tK« 
c-l shooM U multlpuT^t 

uitipiied by the naas ratios listed below: 


! PrODPl1nn«- 


Hg 

n/M Argon 

5.022 

Xe 

3.287 

Kr 

2.098 


For example, at 4000 sec, VB - 404 V a . 

the same n^, 02 , and y) . After the nr ^ <assumlng 
the beam current per module I u cal^T the beam voltage, 
(beam) power as assumed module 


IjjQD " 


(C-4) 


pervoaoc. aqulvalent^r'tta ““-"InS an optica 

— rP propanant. It la Ll “ 

- fotal ratio of o.f. a Jrcnrlt ^ “T" t “““ '' “ 

from a SO-cm-dlamater throatar Or < -33 A can ba axtractad 

“ter. Or. In general, at R . 0.7. 

i^oD ■ vcoef/ vb y-'3 


\ioooy 



(C-5) 
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where the VCOEF for different propellants is listed below: 


Propellant 

VCOEF. A 

VCRIT. V 

Hg 

2.33 

3462 

Xe 

2.88 

3005 

Kr 

3.60 

2590 

Ar 

5.20 

2026 
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The parameter VCRIT Is the voltage at which the extracted beam current 
is the maximum possible for a given thruster diameter because of thermal 
limitations. For the 30-cm-diameter thruster* a thermal analysis 
predicted a maximum possible beam current of 15 A; in general 

^OD " 

The computer program solves for the thruster diameter at a given VB and 
^OD “^^”8 Eq. C-5 for VB < VCRIT and Eq. C-6 for VB> VCRIT. The 
parameters VCRIT and VCOEF for Ar'*' become 12.04 A and 1,158 V for R = 

0.4 and 3.57 A and 2,605 V for R « 0.9. 

After the thruster diameter la calculated, the thruster mass XMASS 
la calculated using the empirical for-iula 

XMASS = 0.078 (kg) , (C-7) 


and the PPU specific mass (kg/kW) is calculated from the empirical 
formula 


ALPPD - 10 [--“3 log (PMOD) + 1.153] _ 

where PMOD is in kW. 


CC-8) 


From the assumed velocity change DELV and initial mass XMO, the 
final mass XMF and propellant mass XMP are calculated using the rocket 
equation. The jet power PJET is then calculated from the propellant 
mass and trip time SEC as follows: 


PJET » 


XMP [ g * XISP) ^ 
2 (SEC) 


(C-9) 
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Ih. total input po«« mi is cslculstsd by dividing the Jit power by 
the product of the thruster, solsr cell, end the m efficiencies. PPU 
efficiency is assumed to be constant at Os 95. Thus, 


-1 


PTOT - PJET (EFF * SOLEFF * .95) » 


(C-10) 


where the thruster efficiency is given by 


0.84 VB 


EFF - (yg + 220) 


(C-11) 


E,.C-11 sssunes . dischsrge pouer of 200 eV/lon end e fixed loss of 
20 eV/lon. The 0.84 lector In Eg. C-11 results from using the sene ness 
efflelency end loss lectors as for the 30-en thruster et 2 A of Hg 
been current. The solar cell efficiency Is detemined from the total 

electron flux (input) and a degradation model. 

The number of operating modules NMOD is calculated from the jet 

power by dividing the latter by the assumed module beam power. 


NMOD - PJET/pjjQjj 


(C-12) 


Stand-by redundancy for both the thrusters (ROT) and power processors 
(ROP) is Included In the nodeling. Thus MOP. the total nutber of PPO 

modules, is given by 

HOP - IMOD [1 + EOP] . 

which asaunes that the life of the PPU Is much greater th«. that of the 
nlsslon. The nunber of thruster nodules. NOT. Is caloulated fron the 
thruster llfetlne LT and the total trip tine TOTRIP as 


not w HMOD /tOTOIP.^ ^ + gol) 


(C-14) 

For most cases, thruster llfetlne Is taken as a constant equal to lo'' 
hr. However, the conputer program has an optional empirical foMula 
that calculates thruster Ufetlne as a function of the bean current 
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density. In Eq. C-lA, the total trip time is obtained from the sum 
of the down time (DOWN) and the up trip time (TIME). This calculation 
is done iteratively since, the down time depends on the down propellant 
mass, MPDWN, which is not known until the vehicle mass, MORV, is cal- 
culated for the up trip. Thus, 

MPDWN ■ MORV [1 - J • 

In Eq. C-15, it is assumed that the same GISP, DELV, and PTOT are 
used for the down trip and for the up trip. 

After the above calculations are made, it is a straightforward 
algebraic calculation to calculate the other system masses and coats. 
Since these formulas are in the computer listings and have been 
discussed in the text of this report, they are not repeated here. A 
summary of the mass and cost nomenclature used for both the orbit- 
raising and on-orbit models is shown in Table C-7 . 


15) 
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Table C-7. Summavy ef Komenclature for Coat Hod^la 


M/^S DEFINITIONS 


Orbit 

Raising 

MRT 

On 

Orbit 

MOT 

Thrusters » Gimbals* etc* 

MRPP 

MOPP 

PPU 

MORP 

MOP 

Proppellant (up) 

MPOVIN 


Propellant (down) 

MRPT 

MQPT 

Tankage 

MRSM 

MOSM 

Structures and Mechanisms 

URSS 

MOSS 

Subsystems 

MRU 

MOW 

Power 

MORV 

MOOP 

Vehicle 

MSP 

MOSP 

Payload 

Orbit 

Raising 

CRT 

On 

Orbit 

COT 

COST DEFINITIONS 

Thruster Costs 


COP 

FPU Costs 

HM 

COPT 

Propellant and Tankage 

CPROP 


Propellant 

CTANK 

— 

Tank 

CRSM 

COSM 

Structures and Mechaalsms 

CRSS 

coss 

Subsystems 

CRU 

COW 

Power 

CORD 


DDT and E Cost 

COPT 

r3PS 

Time Dependent Cost 

CORI 

COI 

Integration and Testing 

COPF 


Costs not associated with electric ORV 

Cl 


Launch Cost 



CTO 

Earth to Orbit Cost 
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appendix d 

RELIABILIIY AND REDUNDANCY 

Calculations were made to determine if 
presently used In the e «. j redundancy factors 

practice, ™<, « aTslTa^^TlTbeC 

The anelyeic that follow, ,how, that the PPO ,„d th “"Pared. 

factor, ehould be approximately -O.IS to 0 20 to 

However, chaope, 1 „ redwndaacy factor, produce onlyTalu 

rha coat per Ulopram of payload calculated by the co,:tdel 

A. ANALYSIS 

- ~”=: “r.r.r.r r-x.™ r 

rate of a single unit la A (fan ^ T 

aloele unit operatlnp for t hour, wUl“e « 

R(t) - e 

This exponential law is derived on the 

random (i.e., chance failures) and th failures are 

the Polaaon proceaa. Mathematically thl,°meV'°tb°““ 

o failure, ocourrlup lu ,^::i" — 1' ^raTr,. 


P(n, T) - e 

n: 


(E-.2) 
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Figure D-1 shows eq. D-S plotted as a function of the number of 
operating units m for a redundancy of 0.15 so that the number of 
standbys was given by 


n •» 0,15 m (Integer part) . (0,-4) 

The sawtooth behavior In the computed reliability Is caused by the 
step function changes in the number of standby units. For example, 
as m varies from 14 to 19, the FPU reliability decreases because 
the number of standbys (two) remains constant while at m - 20 the 
standbys increase to three, producing a sudden Jump in reliability. 

The results shown in Figure D-1 were computed for t ■ 10^ hr 
and an assumed X of 10 ^ and 5 x 10 ^ for the thruster and PPU, 
respectively. We are also assuming the reliability of the switching 
necessary to turn on a standby is included in these failure rates. 

Figure E-1 shows that above approximately 200 units the total thruster 
reliability is greater than 0.99, while for the PPU this reliaMlity 
is achieved at about 40 units. 

Reference to Table D-1 or D-2 shows the number of operating 

units for a wide range of module power and I and indicates that the 

sp 

number of operating modules varies from 'v 2,000 to 25 depending 
on the module power 'and 1^^. The results shown in Figure D-1 imply 
a higher redundancy than 0.15 would be needed for the lower number of 
modules. However, comparing the data in Tables D~1 and D-2 shows that 
the cost/kg changes less than IX for more than a 100% change in PPU/ 
thruster redundancy. Thus, the use of a constant redundancy factor 
produces a total thruster and PPU reliability which varies depending 
on the module power level and beam current density. Henceforth, a 
PPU and thruster reliability of 0.2 will be assumed that will provide 
a worst case reliability of > 0.99 for nearly all cases. 
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Table D~2. Tabulation of Orbit Raising Cost Model Paraiaeters as a 

Function of I with Module Power as a Parameter • 

®P 

Same input data as Table I, except: 

Thruster redundancy = 0»2 




PPO redundancy ®= 0.2 


rs 

Y 

sec 

Number of 

Total 

Thruster 

Thruster 

P 

Payioad 

Payload 

Operating 

Number of 

Diameter, 

Life, 

total 

Mass 

10-3 X kg 

Cost, 

Units 

Thrusters 

cm 

hr 

MN 

$/kg 

2t}Kvi 









e Ji 

104 4 

182.6 

17396.1 

23.2 

467.4 

2 31.8 

6000. 

97C 

7530 

66.3 

3711,5 

28,7 

463.3 

291.2 

3000. 

1310 

30425 

32.3 

1240.3 

35.4 

40 3.2 

429.8 

10000. 

1650 

58956 

21.8 

806.1 

42.7 

321,1 

645.6 

12000. 

paOD= 50KW 

19 S0 

71118 

18,2 

806.1 

50.3 

229.0 

1020. 3 

4 000 . 

252 

417 

268.8 

17396.3 

23.2 

505.5 

194.6 

6000 . 

388 

3012 

104.8 

3711.5 

28.7 

505,6 

2 36.6 

3000. 

524 

12170 

51,1 

1240.3 

35.4 

453.9 

331.6 

Io000. 

660 

23586 

34.5 

806.1 

42.7 

381.6 

471,6 

12000» 

•’ye 

28447 

28.7 

806.1 

50.3 

300.0 

674.0 

Pt*:OD=i0O[\W 

4000. 

126 

208 

408.4 

17396,3 

2 3.2 

525.2 

178.1 

6000. 

194 

1506 

148.2 

3711.5 

26.7 

52'',! 

212.8 

8000. 

262 

6085 

72.2 

1240.3 

35,4 

479.4 

290.6 

lBb00. 

330 

11793 

46.7 

806.1 

4 2.7 

412,1 

402,2 

±2000. 

398 

14223 

40.6 
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50,3 
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555.6 
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4000. 

63 
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23.2 
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97 
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20 s. 6 

3711.5 

28.7 

542.3 
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131 

3042 

102.1 

1240.3 

35.4 

497.5 

264.6 
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165 

5896 

68.9 

806.1 

42.7 

4 33.6 

359.6 

12000. 

199 

7111 

57.4 

806.1 

50.3 

361.1 

486.5 

FWOL!=500KW 

400U. 

25 

41 

913.2 

17396.3 

23.2 

552.2 

158,0 

6000. 

38 

301 

3 31.4 

3711.5 

28.” 

556.0 

184.6 

8000. 

52 

121’’ 

161.4 

1240.3 

35.4 

SI 3. 6 

24 3.4 

10000. 

66 

2 353 

109,0 

806.1 

*42,'’ 

452.8 

325.7 

12000, 

79 

2844 

90.8 

806.1 

50.3 
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appendix e 

ON-ORBIT MISSION AV MODEL 


As part of the on-orblt mission cost modeling analysis. It was 
necessary to estimate AV requirements. Thrust la required for both 
station-keeping (orbit perturbations) and for attitude control. 
Obviously, the magnitude of these forces and torques Is highly’ 
dependent on satellite mass and shape. Based on a review of space 

power satellite (BPS) studies, the following equation for total Impulse 
per year (N-aec/yr) was developed: 


Impulse - 50 M_^ + 150 A + 6 x 10^ 
s a 

where M^ - satellite total mass, kg 

Ag = satellite area exposed to sunlight, 
K = eartti gravitational parameter 
= 4 X 10^ km^/sec^ 



Rg - satellite orbit radius. 


The first term represents north-south stationkeeping, the second 
erm accounts for solar pressure (orbit and attitude perturbations), 
and the third term accounts for gravity gradient torques. The sole 
purpose of this equation Is to estimate the magnitude of the dV 
requirements for a range of satellite sites and is not considered to be 
more than a rough estimate. The details Involved In obtaining this 
equation are discussed in this section. 


A. STATIONKEEPING PERTURBATIONS 

Of the many possible perturbations that might affact a satellite 
orbit, four were considered: (1) east-west forces, caused by 
"triaxlallty;" (2) north-south forces, caused by solar and lunar 
perturbations; (3) solar pressure; and (4) microwave emission from an 
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ant Hast-„nst forces and nlnrowave radiation „are naglactad aa 

being small compared with others. 

per yea? “"f ^ ® «0ulra»ent of ahout 50 „/aac 

year. Oalng cha relationship between average thrust and AV. 


T = M ^ 
s At 


and 


(E-2) 


(E-3) 


Impulse = T At 
Thus, the first term in Eq. E-1 is simply 50 M 

s 

that for ^ ^°>^ce on typical arrays is approximately 55% of 

for a perfect reflector. This translates into a pressure of 

P^p = (0.55) (2) (4.4 X 10-6) = 4.7 x lo'^. N/m^. (e_ 4 ) 

The force associated with solar pressure is 


Fgp = 4.7 X 10-6 A , N , 


(E-5) 


and the annual impulse needed to balance this force 


is approximately 


^sp ~ ^s* N.sec/yr. 


(E-6) 


B. ATTITUDE CONTROL PERTURBATIONS 


The two walor attitude p.rturbatlone coneid.red were gravity 
gra lente and 8»lar-preeeure/CG-mlsallgnwent. Gravity gradient tor 
ate l„portant If the wowente of inertia of the eateuL“ e 

on eat mte L°fT"'"r" 

aatelllte! ^n-tia for large-area 

(a.g., SPS) are significantly different. In addition when 

arge waseea and large dlwenelone are Involved, the gravity gradient 
torques are high. srauient 
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Clearly, without a given satellite design, gravity gradient 

torques can only be grossly estimated. The order of magnitude of this 
torque was assumed to be 


where 


Torque = M A 
1^ s o s 


(E-7) 


The angular 
mately 


^ - satellite orbit rotational frequency 

1/2 



7.3 X 10 rad/sec at synchronous altitude, 
momentum per year associated with this torque is approxl- 


Momentum (GG) = 3 x 10^ 



N»m»sec/yr . 


(S-8) 


The impulse required to compensate this momentum can be estimated by 

assuming propulsion forces are applied with a lever arm (LA) length 
equal to 


LA = 1/2A , 

s * 


(E-9) 


such that 


Impulse (GG) = 6 x 10^ 



(E-10) 


Solar prasaura torqua results from a misalignment of the tenter 
of pressure and the satellite center of mass. If the misalignment Is 

« of a typical satellite dimension (l.e., 0.01 A . the momentum 
imparted would be 


Momentum (sp) = l.SA^^^ , N-m-sec/yr . (E- 

Since misalignments probably could be kept to a few percent of typical 
dimensions, this perturbation would be small compared with others. 
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C. AV ESTIMATE 

An estimate of the propulsion AV requirement 

fi 7 9 

Eq. E-1. Using an M =10 kg and an A = 10 

s s 

requirement per year at synchronous altitude is 

9 

Impulse = 1.65 X 10 , N*sec/yr , 


or 


AV = 1.65 X 10 , m/sec/yr. 


can be obtained from 
the total impulse 

(E“12) 


(E-13) 


This AV value was used in the on-orbit cost model calculations. 



